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ABSTRACT 


Lower  confidence  limit  estimation  procedures  for  the  reliability  of  several 
systems  are  developed  and  their  accuracies  evaluated  using  computer  simulation. 
The  procedures  use  test  data  on  components  of  the  system  which  can  have  failure 
times  with  either  exponential  or  Weibull  distributions  or  both.  Testing  scenarios  for 
the  components  can  be  truncated  by  number  of  failures  or  by  planned  test  times. 

Although  the  evaluation  effort  was  focussed  on  series  systems  in  this  thesis,  the 
procedures  readily  apply  to  other  systems  as  described  in  the  thesis.  The  evaluations 
demonstrate  the  procedures  to  be  quite  accurate  when  sufficient  component  testing 
is  performed. 

Two  FORTRAN  computer  programs  were  written  to  perform  the  evaluation. 
They  are  annotated  in  Users’  Guides  and  can  be  used  to  determine  the  accuracy  of 
these  approximate  lower  confidence  limit  procedures  for  a  given  specific  system  and 
associated  set  of  input  parameters. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  the  computer  programs  developed  in  this 
research  may  not  have  been  exercised  for  all  cases  of  interest.  While  every  effort  has 
been  made,  within  the  time  available,  to  ensure  that  the  programs  are  free  of 
computational  and  logic  errors,  they  cannot  be  considered  validated.  Any  application 
of  these  programs  without  additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


This  thesis  develops  approximate  lower  confidence  interval  procedures  for  the 
reliability  of  complex  systems  using  test  data  on  components  of  the  system.  The 
accuracies  of  these  procedures  are  also  assessed  using  computer  simulation.  The 
procedures  can  be  used  for  any  complex  system  whose  reliability  does  not  decrease 
when  the  reliability  of  any  one  of  the  components  is  increased. 

The  failure  times  of  the  continuously  operating  components  are  assumed  to 
have  either  an  exponential  or  Weibull  distribution.  Parameters  in  both  distributions 
are  assumed  to  be  unknown.  The  Weibull  distribution  is  used  to  model  the  lifetime 
probability  distribution  of  electronic  components  with  non-constant  failure  rate 
functions.  It  is  also  used  to  model  the  lifetime  probability  distributions  of  mechanical 
devices,  since  their  failure  rate  functions  are  usually  increasing  with  operating  time. 

Lower  confidence  limit  estimation  procedures  for  system  reliability  are  needed 
during  the  development  phase  of  systems  to  provide  indications  of  a  contractor’s 
ability  to  meet  a  stated  system  reliability  goal  as  development  progresses  and  the 
results  of  test  progr^^’c  become  available.  "Hiese  procedures  are  also  needed  to 
assess  the  reliability  of  systems  that  have  been  operating  in  the  field  for  some  time 
and  have  accumulated  histories  of  failure  data  and  unique  configurations  of  modified 
or  repaired  components. 

Few  textbooks  on  reliability  treat  the  problem  of  system  reliability  interval 
estimation.  Those  that  do  usually  limit  the  discussion  to  series  or  parallel  systems. 
Moreover,  the  procedures  they  present  are  not  adaptable  to  other  more  complex 
systems.  Mann,  Shafer,  and  Singpurwalla  [Ref.l  pp  487-524]  provide  one  of  the  better 
treatments  of  a  variety  of  these  methods  in  Chapter  10  of  their  book.  This  chapter 
provides  an  excellent  summative  discussion  of  the  many  procedures  that  were 
developed  from  1954  to  1974.  However,  none  of  the  procedures  reviewed  in  their 
book  can  accomodate  the  use  of  test  data  from  a  mix  of  components  with  both 
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exponential  and  Weibull  failure  time  distributions.  The  procedures  presented  in  this 
thesis  does  accomodate  this  type  of  system  with  a  mixture  of  different  component 
types.  In  addition,  the  procedures  presented  in  this  thesis  can  accomodate  an 
additional  mix  of  components  for  which  only  attribute  data  has  been  collected. 

Procedures  developed  in  this  thesis  are  extensions  of  a  procedure  developed  by 
Mylire,  Sanders  and  Rosenfeld  [Ref.2].  In  their  paper,  they  assume  the  failure  times 
of  continuously  operating  components  have  exponential  distributions  with  associated 
failure  rates,  X,.  The  test  data  on  the  remaining  components,  the  number  of  observed 
failures/]  in  n,  tests,  are  assumed  to  have  Poisson  distributions  with  associated  means 
ni.q,.  They  assume  the  ratios  X/Xj ,  q/qj  and  X/qj  are  known  and  develop  confidence 
interval  estimation  procedures  for  system  reliability  that  use  this  information.  They 
also  show  that  the  accuracy  of  their  procedure  is  not  very  sensitive  to  moderate 
inaccuracies  of  these  ratios.  This  suggests  that  it  might  be  possible  to  estimate  the 
ratios  from  the  data  as  part  of  the  interval  estimation  process  and  not  suffer 
significant  loss  of  accuracy  in  the  interval  estimates.  Estimating  these  ratios  is  part 
of  the  procedures  developed  in  this  thesis. 

This  thesis  also  provides  an  annotated  computer  program  that  can  be  used  to 
assess  the  accuracies  of  the  lower  confidence  limit  procedures  when  applied  to  any 
specific  system.  Sufficient  annotations  are  provided  throughout  the  program  in 
Appendix  C.  This  program  provides  the  user  with  a  means  for  verifying  the  accuracy 
of  these  proposed  lower  confidence  limit  procedures  for  his  specific  system  and 
testing  program,  that  is,  sample  sizes  and  type  of  truncation.  This  capability  will 
allow  the  user  to  answer  many  "what  if  type  of  questions. 
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II.  THEORY 


A.  Interval  Estimation  Procedure  for  Exponential  Failure  Times 

A  system  is  defined  to  be  quasi-coherent  if  an  increase  in  reliability  of  any  one 
of  its  components  does  not  cause  a  decrease  in  system  reliability.  The  components 
of  a  quasi-coherent  system  do  not  need  to  be  statistically  independent.  However, 
throughout  this  thesis,  it  is  assumed  that  all  components  are  statistically  independent. 

Suppose  a  quasi-coherent  system  has  k  components  and  the  distribution  of  the 
failure  time  of  component  i  is  exponential  with  failure  rate  Then  the  system 
reliability  can  be  written  as  a  function  of  ,i  =  1,  2,  ...,  k  as  follows: 

^.(0  =g(  Aj.X,,. . t,)  ...(2.1) 


where  r,  is  the  operating  time  for  component  i.  Let  m  be  any  one  of  the  k 

components  and  r,  =  ,  for  /  =  1,  2 . k.  Then  equation  (2.1)  may  be  viewed 

as 


^.(0  =  g(  Kr  tj.t^ . t,  )  ...  (2.2) 

If  the  r/s  are  known  and  were  an  upper  100(l-a)%  confidence  limit  for 

the  corresponding  lower  confidence  limit  for  RJt)  would  be: 

^j(0l(o)  ~  si  )  "■ 

Specifically,  if  we  have  a  series  system  of  independent  components,  so  that 


Rsit)  =  exp(  ^ 

1-1 

=  expi  ) 


...  (2.4) 


1-1 


then. 


^sif)L(a)  =  exp( 


1=1 


...  (2.5) 


If  AJ,  items  of  component  /  are  tested  until  /  failures  occur,  7,  denotes  the  total  test 

k 

time  accumulated  on  all  the  //,  items,  and  F  =  E-^i  expression 

(-1 

2\  hr  T 

m  t  i 
i*l 


has  a  Chi-square  distribution  with  2F  degrees  of  freedom.  See  Bain  and  Engelhardt 
[Ref.3].  The  corresponding  100(l-a)%  upper  confidence  limit  for  is 


m,U(a) 


A  tt,2f 


1-1 


...  (2.6) 


where  is  the  100(l-a)th  percentile  point  of  a  Chi-square  distribution  with  2F 
degrees  of  freedom. 

If  the  testing  on  component  i  is  terminated  when  a  total  test  time  of  T,  has 
been  accumulated  by  all  n,  items,  then  the  equation  for  becomes 


m.U{a) 


X~tt.2(UD 

2E'-,r, 

»-l 


...  (2.7) 


In  this  case  /,  is  random  and  so  is  F. 

If  testing  on  each  of  the  n,  items  of  component  i  are  tested  until  a  planned  test 
time  or  failure,  and  failed  items  are  replaced  immediately,  then  equation  (2.7)  will 
be  the  exact  expression  for  .  If  failures  are  not  replaced,  then  equation  (2.7) 

is  approximate.  See  Lee,  Bain  and  Englehardt  [Ref.3  pp  486-495].  Department  of 
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Defense  document  NAVSEA  OD29304B  "Reliability  and  Availability  Evaluation 
Program  Manual"  [Ref.4  p  5-42]  provides  nearly  exact  procedures  for  when 

testing  is  terminated  by  planned  test  time  for  each  item  tested  and  failures  are  not 
replaced. 

The  values  of  the  r/s  are  assumed  to  be  unknown  in  this  thesis.  When  testing 
is  terminated  by  the  number  of  failures,  a  nearly  unbiased  estimator  for  r-  is 


where  A,  =(/,  -!)/ 7’.  and  the  index  m  denotes  the  component  with  largest  value 

of  A .  .  The  ratio,  (f^l )/Ti ,  is  an  unbiased  estimator  for  A,  (see  Appendix  A).  If  1  /  A^ 
were  unbiased  for  1/A^  then  f.  would  be  an  unbiased  estimator  for  r,.  Replacing 

with  in  equation  (2.8)  will  yield  an  unbiased  estimator  f-  for  r,. 

Multiplying  by  this  constant  /„/(/„, -7)  is  nullified  by  a  cancellation  with  the  same 
constant  in  the  final  equation  for  the  system  reliability  lower  confidence  limit,  so 
equation  (2.8)  is  used  to  estimate  r,.  Using  estimator  r.  for  r,  ,  equation  (2.6) 

becomes 


X'«.2r 


I«1 


...  (2.9) 


It  is  important  to  note  that  the  index  m  denotes  the  component  for  which  X- 

=  (fi-l)/Ti  is  the  largest  among  all  the  components  in  the  system.  The  corresponding 
equation  for  the  100(l-a)%  lower  confidence  limit  on  the  reliability  of  a  series 
system  is 
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...  (2.10) 


* 

/=1 

The  corresponding  lower  confidence  limit  for  the  reliability  of  any  quasi- 
coherent  system  is  given  by  equation  (2.3)  with  r,  replaced  by  r-  . 

In  this  thesis,  equation  (2.8)  with  i.  =  fJT-  will  also  be  used  to  estimate  r, 

under  exponential  assumptions  when  testing  is  terminated  after  an  accumulated  test 
time  is  achieved  (tnincated),  and  when  at  least  two  components  have  at  least  one 
failed  test  item.  This  is  done  because  failures  will  not  be  replaced  in  the  time 
truncated  test  plans  that  are  simulated  in  this  thesis.  In  this  type  of  testing  both  fj 
and  7,  are  random.  Under  time  truncation,  it  is  possible  that  no  failures  will  occur 
on  any  component  tested  in  which  case  equation  (2.8)  is  undefined.  Also,  if  only  one 
component  has  one  failure  and  the  remaining  components  have  zero  failure,  equation 
(2.8)  would  be  zero  for  all  i  except  the  case  when  i  =  m. 

All  of  the  confidence  limit  procedures  in  this  thesis  have  a  common  special 
method  for  computing  the  lower  confidence  limit  of  system  reliability  in  the  two 
cases  of  zero  or  one  failure.  This  feature  amounts  to  a  modification  to  equations 
(2.3)  and  (2.10). 

When  either  zero  failures  or  one  failure  have  occured  among  all  components, 
the  test  data  is  examined  for  each  component  to  determine  the  total  number,  A, ,  of 
equivalent  component  mission  tests  (/  =  1,  2,  ...,  k).  These  A/s  and  the  system 
configuration  are  analyzed  to  determine  the  equivalent  number  of  mission  tests.  A, 
for  the  system  that  would  have  occurred  if  A; ,  A, ,  ...,  of  these  k  components 
were  assembled  into  systems. 

For  a  series  system,  this  A  will  be  equivalent  to  min  {  Nj ,  N2,  N,^  }.  The 

100(l-a)%  lower  confidence  limit  of  system  reliability  if  zero  failures  occurred  is  then 
computed  directly  as  follows: 
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...  (2.11) 


If  exactly  one  failure  occurred  among  all  k  components, 
solution  for  p  in  the  equation 

+  Np^'^{\-p)  =  a  —  (2-12) 

These  confidence  limit  equations  are  the  standard  binomial  lower  confidence  limit 
equations.  Equations  (2.11)  and  (2.12)  are  part  of  the  set  of  equations  used  to 
compute  ^/Ol(b)  truncated  interval  estimation  procedures  in 

this  thesis. 

It  is  important  to  remember  that  the  symbol  7,  in  equations  (2.7),  (2.8)  and 
(2.9)  denote  total  accumulated  test  time  for  component  i  ;  that  is 

T  ~  T  "■  (2-13) 

/'I 


where  Tjj  is  the  test  time  accumulated  on  the  jih  test  item  of  component  i  and  n,  is 
the  number  of  test  items  of  component  i  being  tested. 


B.  Interval  Estimation  Procedure  for  Weibull  Failure  Times 

Consider  a  series  system  with  k  components.  Let  the  time  to  failure,  Xj  ,  of 
component  i  have  a  Weibull  distribution  with  density 


/,(',)  =  )  ,  (,  >  0 


...  (2.14) 


Then 


/?,((,)  =  exp(  )  ,  (,  >  0  -  (2-15) 
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and 


=  n  exp{-xf‘f/‘} 

i-l 

.exp{-S  -(.2.16) 

I-l 

=  exp{-A;52r.f/'} 

I-l 

where  A;*  =  A;^'  ,  A„*  is  any  one  of  the  \  ’ ,  i  =  1,  2,  k  ,  and  r,  =  the 

ySj’s  are  known,  then  A'/'  will  have  a  constant  failure  rate  A/'  and  the  procedures 
described  in  Section  A  can  be  used  to  obtain  ^^(0/,(b)  "'^th  replaced  by  in 

equation  (2.13). 

Suppose  0,  is  unknown  and  ,  AQ^,^ ,  ...,  are  the  ordered  failure  times 
under  either  type  of  truncated  testing  for  component  i  in  the  system.  Solutions 
and  X.  for  0i  and  A,  in  the  two  equations  given  in  equation  (2.17)  are  the 


maximum  likelihood  estimates  for  and  A,.  See  Mann  and  others  [Ref.l  pp  189-191]. 
These  equations  are  used  for  both  types  of  test  truncation.  If  for  component  i, 
testing  is  terminated  on  the  /;'*  failure,  then  t^^  =  A^^  in  equation  (2.17).  The 
solution,  3,  ,  is  a  biased  estimator  for  Bain  [Ref.5  pp  220]  provides  a  table  of 

constants  BinJ  which  depends  on  number  of  test  items  n,  such  that  =  3,^(«,) 


is  a  nearly  unbiased  estimator  for 


M _ 

*  («,-/, )'il' 


fi 


j  =  7E 

Pi  Ji  ;>1 
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and 


...  (2.17b) 


If  the  testing  for  component  i  is  terminated  at  failure  or  at  a  given  time  for 
each  of  the  /i,  items  on  test,  then  4  =  in  equations  (2.17a)  and  (2.17b). 

Now,  let 

Tij  =  ,  with  i  =  1,  2 . k  (2.18) 

j  =  1,  2,  ...,  n. 


In  this  thesis,  the  distribution  of  is  approximated  by  the  exponential  distribution 
^  « 

with  failure  rate  ‘  s  and  procedures  similar  to  those  in  Section  A  are  used 


to  obtain  the  lower  confidence  limit  on  system  reliability.  Define 


(•  _  fi 

'  ~  T: 


...  (2.19) 


where  7,.  =  J]  7^  ,  i  =  1,  2,  ...,  k  .  Let  X’  =  max^. 
Note  that  this  defines  the  index  m.  Define  r.  by 


■  I  ■ 


...  (2.20) 


for  both  types  of  test  truncation  plans. 

Then  an  approximate  100(l-a)%  upper  confidence  limit  for  is  given  by 
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...  (2.21) 


m,U(_a) 


2E  f,T, 


where 


F'  = 


k 


Hfi 


k 


1  *  E/, 

i-l 


if  test  till  fi  **  failure 
for  all  components 

if  test  till  specified  time 
’  for  all  components 


...  (2.22) 


The  corresponding  approximate  100(l-a)%  lower  confidence  limit  Rs(0L(tt)  for  the 
reliability  R/t)  of  a  series  system  is  given  by 


=  exp{  -il 


m,U(a.) 


EV/ 


...  (2.23) 


I-l 


when  at  least  two  components  have  at  least  one  failure.  Equations  (2.11)  and  (2.12) 
also  apply  here  when  the  total  failures  over  all  components  is  either  zero  or  one. 

The  accuracies  of  these  approximate  confidence  interval  procedures  were 
evaluated  by  using  computer  simulations  which  are  described  in  the  next  chapter. 
During  this  evaluation  process,  the  degrees  of  freedom  in  the  expressions  X^a,2F‘ 


in  equation  (2.21),  X^^2F  in  equation  (2.9)  and  in  equation  (2.7)  were 

increased  and  decreased  from  the  defined  values  of  F‘  and  F  given  by  these 
equations.  The  purpose  of  these  modifications  was  to  find  more  accurate  lower 
confidence  limit  procedures.  The  specific  increases  and  decreases  are  described  in 
Chapter  III.  The  results  show  that  for  some  cases  the  procedures  with  modified 
degrees  of  freedom  are  more  accurate. 
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III.  COMPUTER  SIMULATION 


A.  Test  Plan  1  :  Testing  n,  Until  ^  Failures  (RETPl) 

RETPl  is  a  program  written  in  FORTRAN,  on  the  Amdahl  mainframe 
computer,  which  performs  the  computer  simulation  of  the  random  failure  times  of 
the  different  types  of  components  in  the  system.  A  documentation  of  this  program 
and  its  associated  subroutines  is  included  in  Appendix  B. 

The  program  accepts  input  parameters  via  an  input  file  INl.DAT.  For  each 
replication,  it  generates  the  failure  times  for  all  the  component  items  included  in  the 
test  plan  using  a  uniform  random  number  generating  subroutine  LRNDPC.  A  quick 
evaluation  of  LRNDPC  (see  Appendk  D)  by  plotting  U(n  +  1)  vs  U(n)  illustrates  the 
uniformity  of  the  routine.  The  program  determines  the  total  test  time  accumulated 
for  each  component  in  the  system  and  computes  the  estimates  of  the  key  parameters 
and  the  consequent  lower  confidence  limit  for  system  reliability  for  that  replication. 
The  process  is  repeated  1000  times.  When  all  replications  are  done,  the  routine 
EVAL  processes  the  lower  confidence  limit  estimates  from  all  1000  replications  and 
determines  the  two  measures  of  accuracy  for  the  run,  namely  RSLOW  and  LEVEL. 

RSLOW  is  the  100(1-Q!)  percentile  of  the  ordered  set  of  lower  confidence  limits 
from  the  1000  replications  computed  in  a  run.  The  true  reliability  of  the  system  is 
RS.  The  closer  RSLOW  is  to  RS,  the  greater  the  accuracy  of  the  procedure  under 
evaluation  in  the  run.  If  the  procedure  is  exact,  RSLOW  will  be  equivalent  to  RS. 
To  be  conservative,  RSLOW  should  always  be  lower  than  RS. 

LEVEL  measures  the  proportion  of  1000  lower  confidence  limits,  from  a  run 
with  1000  replications,  which  are  lower  than  the  true  system  reliability  RS.  The 
closer  LEVEL  is  to  the  specified  confidence  level  for  the  procedure,  l-a,  the  better 
the  procedure.  Values  of  LEVEL  greater  than  (l-a)  reflect  an  under-estimation  of 
RS  which  is  conservative.  Values  of  LEVEL  lesser  than  l-a  signal  an  over¬ 
estimation  of  RS  which  may  be  undesirable. 
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Simulation  runs  are  performed  using  RETPl  for  all  combinations  of  failure 
time  distributions  and  levels  of  key  input  parameters  listed  below. 

(a)  System.  -  8  Exponential  components  in  Series  (Case  1) 

-  8  Weibull  components  in  Series  (Case  2) 

-  4  Exp  and  4  Wei  (Mixed)  components  in  Series  (Case  3) 

(b)  True  System  Reliability  (RS). 

-  Hi  (greater  than  0.9)  (Type  A) 

-  Lo  (greater  than  0.8)  (Type  B) 

(c)  Level  of  Significance  (a). 

-  0.1 
-  0.2 

(d)  Degrees  of  Freedom  for  statistic  (DF)  as  a  function  of  the  total  number  of 
failed  test  components  (NFC)  and  total  number  of  system  components  (NCOMP). 

-  DF  =  2  •  NFC 

-  DF  =  2  •  (  NFC  +  NCOMP  ) 

-  DF  =  2  •  (  NFC  -  NCOMP  ) 

-  DF  =  2  •  NFC  -  NCOMP 

e)  Test  Plan  for  each  component. 

-  Test  5  until  5  failures 

-  Test  15  until  15  failures 

-  Test  15  until  11  failures 

-  Test  15  until  7  failures 

-  Test  15  until  3  failures 


For  the  8  exponential  components  in  Case  1,  the  mission  time  for  each  of  the 
component  is  chosen  to  be  10  hrs.  The  program  will  accomodate  different 
component  mission  times.  The  chosen  values  of  the  scale  parameters,  X;  ,  were 
different  depending  on  whether  the  system  is  highly  reliable  (Type  A)  or  one  with 
a  lower  reliability  (Type  B).  The  ratios  between  the  largest  and  the  smallest  failure 
rate  was  chosen  to  be  8  and  4.5  respectively  for  Type  A  and  Type  B  systems. 
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For  the  8  Weibull  components  in  Case  2,  the  mission  time  for  each  of  the 
components  was  chosen  to  be  10  hrs.  The  chosen  values  of  the  scale  parameters,  Xj, 
were  different  depending  on  whether  the  system  is  highly  reliable  (Type  A)  or  one 
with  a  lower  reliability  (Type  B).  The  ratio  between  the  largest  and  smallest  failure 
rate  was  chosen  to  be  8  for  both  system  types.  The  shape  parameter  is  chosen  to  be 
2  for  all  cases.  The  program  will  accomodate  any  value  greater  than  zero  for  the 
shape  parameter. 

A  mixture  of  exponential  and  Weibull  components  with  those  parameters 
described  in  the  last  two  paragraphs  is  chosen  for  the  Type  A  and  Type  B  systems 
of  Case  3. 

Each  simulation  run  of  1000  replication  results  in  an  output  file  OUTl.DAT. 
raw  output  from  all  the  RETPl  runs  are  summarized  in  tabular  form  and  placed 
in  Appendix  E.  Each  table  corresponds  to  a  specific  run  case  and  system  type 
combination. 


B.  Test  Plan  2  ;  Testing  for  a  Specified  Planned  Test  Time  (RETP2) 

RETP2  is  another  program  written  in  FORTRAN,  on  the  Amdahl  mainframe 
computer,  which  performs  the  computer  simulation  of  the  random  failure  times  of 
the  different  types  of  components  in  the  system.  A  documentation  of  this  program 
and  its  associated  subroutines  is  included  in  Appendix  C. 

The  structure  of  this  program  is  quite  similar  to  that  of  RETPl  described  in 
Section  A  of  this  chapter.  The  program  accepts  input  parameters  via  an  input  file 
IN2.DAT.  For  each  replication,  it  generates  the  failure  times  for  all  the  component 
items  included  in  the  test  plan  using  LRNDPC.  The  program  then  determines  the 
number  of  failed  test  components  for  each  component  in  the  system  and  computes 
the  estimates  of  the  key  parameters  and  the  consequent  lower  confidence  limit  for 
system  reliability  for  that  replication.  The  process  is  repeated  1000  times.  When  all 
replications  are  done,  the  routine  EVAL  processes  the  lower  confidence  limit 
estimates  from  all  the  1000  replications  and  determines  the  two  measures  of  accuracy 
for  the  run,  namely  RSLOW  and  LEVEL.  The  definitions  of  these  two  measures 
were  discussed  in  the  Section  A. 

Simulation  runs  are  performed  using  RETP2  for  all  combinations  of  failure 
time  distributions  and  levels  of  key  input  parameters  listed  below. 

(a)  System.  -  8  Exponential  components  in  Series  (Case  4) 

*  8  WeibuII  components  in  Series  (Case  5) 

-  4  Exp  and  4  Wei  (Mixed)  components  in  Series  (Case  6) 

(b)  True  System  Reliability  (RS). 

-  Hi  (greater  than  0.9)  (Type  A) 

-  Lo  (greater  than  0.8)  (Type  B) 

(c)  Level  of  Significance  (a). 

-  0.1 
-  0.2 
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(d)  Degrees  of  Freedom  for  statistic  (DF)  as  a  function  of  the  total  number  of 
failed  test  components  (NFC)  and  total  number  of  system  components  (NCOMP). 

-  DF  =  2  •  (  1  +  NFC  ) 

-  DF  =  1.3  •  2  •  (  1  +  NFC  ) 

(e)  10  values  of  A:  =  0.25,  0.5,  1,  2,  3,  4,  5,  10,  20,  30 

where  K  is  a  factor  such  that  the  expected  number  of  failures  for  an  exponential 
component  during  the  specified  total  test  time  for  that  component  is  0.6  times  K. 
This  K  factor  determines  the  bounds  of  the  expected  total  number  of  failed  test  items, 
E[NFC].  The  accuracy  of  the  lower  confidence  limit  procedures  are  highly  correlated 
with  E[NFC]. 

For  the  8  exponential  components  in  Case  4,  the  mission  time  for  each  of  the 
component  is  chosen  to  be  5  hrs.  The  program  can  accomodate  different  component 
mission  times.  The  chosen  values  of  the  scale  parameters,  Xj,  were  different 
depending  on  whether  the  resultant  system  is  highly  reliable  (Type  A)  or  one  with 
a  lower  reliability  (Type  B).  The  failure  rate  was  chosen  to  be  0.001  and  0.005 
failures/hr  for  all  the  components  respectively.  Total  test  time  to  be  accumulated 
by  each  component  is  computed  according  to  the  following  method.  For  each 
component  i ,  /,  represents  the  amount  of  operating  time  required  to  result  in  a  40% 
survival  probability,  that  is,  an  expected  failure  of  0.6  component  with  an  exponential 
failure  time  distribution.  The  computation  for  t,  is  as  follows; 

^iUi)  =  exp(-Xr.) 

=  0.4 

<,  =  -(l)ln(0.4) 

r, ,  the  total  amount  of  test  time  to  be  accumulated  for  component  i  would  be 
K  times  f,  which  will  result  in  an  expected  number  of  0.6  times  K  failed  items  for  this 
component.  E[NFC]  will  then  be  8  times  that  number,  since  there  are  8  such 
components  in  the  system. 
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For  the  8  Weibull  components  in  Case  5,  the  mission  time  for  each  of  the 
component  is  chosen  to  be  15  hrs.  The  chosen  values  of  the  scale  parameters,  \  , 
were  0.005  failures/hr  for  the  Type  A  system  and  0.01  failures/hr  for  a  Type  B 
system.  The  shape  parameter  is  chosen  to  be  2.  For  each  Weibull  component,  a 
maximum  of  20  test  items  are  tested  to  failure.  Estimation  of  E[NFC]  and  thus  the 
total  test  time  to  be  accumulated  for  each  component  is  based  on  the  exponential 
failure  time  model  described  in  the  earlier  paragraphs. 

A  mixture  of  exponential  and  Weibull  components  with  those  parameters 
described  in  the  last  two  paragraphs  is  chosen  for  the  Type  A  and  Type  B  systems 
of  Case  6. 

Each  simulation  run  of  1000  replication  results  in  an  output  file  OUT2.DAT. 
The  raw  output  from  all  the  RETP2  runs  are  summarized  in  tabular  form  and  placed 
in  Appendbc  F.  Each  table  corresponds  to  a  specific  run  case  and  system  type 
combination. 
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IV.  RESULTS  AND  DISCUSSION 

The  results  of  the  simulation  runs  are  summarized  and  discussed  in  this  section. 
Tables  lA,  IB,  2A,  2B,  3A  and  3B  in  Appendix  E,  and  Tables  4A,  4B,  5A,  5B,  6A 
and  6B  in  Appendix  F  present  the  accuracy  results  in  tabular  form  for  all  run  cases 
that  were  simulated:  a  few  of  these  tables  appear  in  this  section  to  facilitate 
discussion  of  the  results. 

A.  Test  Plan  1  ;  Testing  Until /j  Failures  (RETPl) 

Table  lA  displays  the  simulation  results  for  Case  1  Type  A.  In  this  case,  the 
system  is  comprised  of  8  components  in  series.  The  failure  time  of  each  component 
has  an  exponential  distribution.  The  failure  rates  of  the  8  components  range  from 
0.0002  failures/hour  to  0.0016  failures/hour.  The  mission  time  of  the  system  is  10 
hours  and  the  mission  operating  time  of  each  component  was  also  set  to  10  hours. 
The  component  mission  times  do  not  need  to  be  equal  to  the  system  mission  time 
for  the  procedures  evaluated  in  this  thesis.  This  was  discussed  in  Chapter  II  and  is 
allowed  for  in  all  of  the  lower  confidence  limit  equations.  System  reliability,  RS,  in 
Table  lA  is  0.931  .  Throughout  this  case,  all  of  these  parameters  remain  fixed. 

In  simulation  number  1  (S/N:  1)  in  Table  lA,  five  items  for  each  of  the  8 
components  in  the  system  are  tested  until  they  fail.  Thus  the  number  of  failed 
components  (NFC)  is  40.  One  set  of  this  40  failure  times  is  randomly  generated  for 
each  simulation  run.  For  this  set  of  data,  four  90%  lower  confidence  limits  and  four 
80%  lower  confidence  limits  are  computed.  The  four  limits  correspond  to  the  values 
assigned  to  the  degrees  of  freedom  parameter  F  in  the  symbol  X^b.2F  which  is  a  factor 
in  the  upper  confidence  limit  equation  for  .  The  four  different  methods  for 

computing  this  degree  of  freedom  appear  in  the  "Deg  of  Freedom"  column.  NCOMP 
denotes  the  number  of  components  in  the  system.  Thus  for  each  simulation  run, 
eight  lower  confidence  limits  are  computed.  After  1000  replications  of  this 
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Table  lA  ;  8  Exp  in  Series,  RS  =  0.931  (Hi) 

min  A  =  0.0002  f/hr,  max  X  =  0.0016  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

Measures  of  Accuracy 

a 

RSLOW 

LEVEL 

1 

Test  5  until 

2*NFC 

0.1 

0.919 

0.982 

5  failed. 

(80) 

0.2 

0.919 

0.960 

NFC  =  40 

2»(NFC  + 
NCOMP) 

(96) 

0.1 

0.906 

1.000 

0.2 

0.905 

0.999 

2*NFC- 

NCOMP 

(72) 

0.1 

0.927 

0.949 

0.2 

0.927 

0.880 

2*(NFC- 

NCOMP) 

(64) 

0.1 

0.934 

0.821 

0.2 

0.934 

0.702 

2 

Test  15  until 

2*NFC 

0.1 

0.928 

0.955 

15  failed. 

(240) 

0.2 

0.927 

0.908 

NFC  =120 

2*(NFC+ 

NCOMP) 

(256) 

0.1 

0.923 

0.990 

0.2 

0.923 

0.975 

2*NFC- 

NCOMP 

(232) 

0.1 

0.930 

0.916 

0.2 

0.930 

0.833 

2*(NFC- 

NCOMP) 

(224) 

0.1 

0.932 

0.844 

0.2 

0.932 

0.747 

3 

Test  15  until 

2*NFC 

0.1 

0.927 

0.955 

11  failed. 

(176) 

0.2 

0.926 

0.916 

NFC  =  88 

2*(NFC+ 

NCOMP) 

(192) 

0.1 

0.921 

0.996 

0.2 

0.920 

0.988 

2»NFC- 

NCOMP 

(168) 

0.1 

0.930 

0.916 

0.2 

0.929 

0.843 

0.1 

0.933 

0.843 

0.2 

0.932 

0.735 
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Table  lA  : 


8  Exp  in  Series,  RS  =  0.931  (Hi)  (Cent...) 

min  X  =  0.0002  f/hr,  max  X  =  0.0016  f/hr,  UT  =  10  hrs 


S/N 

Test 

Deg  of 

Measures  of  Accuracy  | 

Plan 

Freedom 

Q 

RSLOW 

LEVEL 

4 

Test  15  until 

2*NFC 

0.1 

0.924 

0.970 

7  failed. 

(112) 

0.2 

0.923 

0.931 

NFC=56 

2*(NFC+ 

NCOMP) 

(128) 

0.1 

0.915 

0.998 

0.2 

0.913 

0.994 

2*NFC- 

NCOMP 

(104) 

0.1 

0.929 

0.919 

0.2 

0.928 

0.853 

2*(NFC- 

NCOMP) 

(96) 

0.1 

0.934 

0.835 

0.2 

0.933 

0.720 

5 

Test  15  until 

2*NFC 

0.1 

0.915 

0.986 

3  failed. 

(48) 

0.2 

0.912 

0.975 

NFC  =  24 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.891 

1.000 

0.2 

0.888 

1.000 

2*NFC- 

NCOMP 

(40) 

0.1 

0.927 

0.944 

0.2 

0.926 

0.860 

■nSlH 

0.1 

0.939 

0.753 

0.2 

0.939 

0.634 

simulation  are  run,  the  2  measures  of  accuracy  RSLOW  and  LEVEL  are  computed. 
The  lower  confidence  limit  procedures  are  exact  if  RSLOW  =  RS  in  which  case 
LEVEL  =  1-a. 

Table  1 A  displays  the  accuracy  results  for  5  different  sampling  plans  which  are 
described  in  S/N:  1,  2,  3,  4  and  5. 

A  comparison  of  the  four  values  of  RSLOW  for  each  of  these  five  sampling 
plans  reveal  that  the  lower  confidence  limit  procedure  with  degrees  of  freedom  equal 
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to  2*NFC-NCOMP  is  the  most  accurate  lower  confidence  limit  procedure.  In  S/N 
1,  for  example,  the  RSLOW  value  of  0.927  is  the  largest  such  value  below  the  RS 
value  of  0.931.  Values  of  RSLOW  above  RS  are  optimistic  and  not  as  desirable  as 
values  of  RSLOW  which  are  equi-distant  below  RS. 

The  values  of  RSLOW  and  LEVEL  are  based  on  1000  replications  and  their 
accuracy  merit  should  roughly  be  measured  to  the  nearest  one  hundredth.  That  is 
we  should  round  0.927  to  0.93  and  compare  it  with  RS  =  0.93.  It  is  evident  that  the 
lower  confidence  limit  procedure  with  degrees  of  freedom  equal  to  2*NFC-NCOMP 
is  very  accurate  for  all  cases  simulated. 

Table  2A  displays  the  accuracy  results  of  Case  2  for  Type  A  systems.  In  this 
case,  the  8  components  cormected  in  series  have  the  shape  parameter  p  -  2  and  the 
scale  parameters  varying  between  0.001  and  0.008  failures/hr.  Mission  time  is  10 
hours  and  each  component  has  this  same  mission  or  utilization  time  (UT). 
Inspection  of  Table  2A  reveals  the  following: 

(1)  More  than  5  items  of  each  component  should  be  tested  until  failure  for 
any  of  these  procedures  to  be  reasonably  accurate. 

(2)  If  15  items  of  each  component  are  tested  until  all  fail,  then  these 
procedures  will  be  reasonably  accurate  when  the  degrees  of  freedom  is  either 
2*NFC-NCOMP  or  2*(NFC-NCOMP). 

(3)  The  procedures  are  reasonably  accurate  for  80%  confidence  level  when 
the  truncation  is  not  below  7  out  of  15  items. 

(4)  The  procedures  are  slightly  conservative  at  the  90%  confidence  level  when 
the  truncation  is  7  out  of  15  items  or  11  out  of  15  items. 

There  are  numerous  ways  to  modify  these  lower  confidence  limit  procedures 
to  effect  improvements  in  their  accuracy.  One  avenue  is  to  modify  the  estimate  for 
the  shape  parameter,  /3.  Some  very  recent  work  in  the  literature  provides  a  method 
for  estimating  P  that  differs  greatly  from  the  maximum  likelihood  estimator  (MLE) 
and  does  not  require  computer  iteration. 
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Table  2A  : 


8  Wei  in  Series,  RS  =  0.980  (Hi) 

min  X  =  0.001  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

1 

Test  5  until 

5  failed. 

NFC=40 

2*NFC 

(80) 

0.1 

0.947 

0.992 

0.2 

0.930 

0.989 

2*(NFC+ 

NCOMP) 

(96) 

0.1 

0.937 

0.994 

0.2 

0.918 

0.993 

2*NFC- 

NCOMP 

(72) 

0.1 

0.951 

0.989 

0.2 

0.937 

0.986 

2*(NFC- 

NCOMP) 

(64) 

0.1 

0.956 

0.985 

0.2 

0.943 

0.981 

2 

Test  15  until 

15  failed. 

NFC =120 

2*NFC 

(240) 

0.1 

0.978 

0.918 

0.2 

0.974 

0.913 

2*(NFC+ 

NCOMP) 

(256) 

0.1 

0.977 

0.931 

0.2 

0.972 

0.924 

2*NFC- 

NCOMP 

(232) 

0.1 

0.979 

0.914 

0.2 

0.975 

0.901 

2*(NFC^ 

NCOMP) 

(224) 

0.1 

0.980 

0.904 

0.2 

0.975 

0.889 

3 

Test  15  until 

11  failed. 

NFC=88 

2*NFC 

(176) 

0.1 

0.982 

0.876 

0.2 

0.977 

0.860 

2*(NFC  + 
NCOMP) 
(192) 

0.1 

0.980 

0.894 

0.2 

0.975 

0.882 

2*NFC- 

NCOMP 

(168) 

0.1 

0.983 

0.861 

0.2 

0.978 

0.839 

HRESaH 

0.1 

0.983 

0.840 

0.2 

0.979 

0.819 
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Table  2A  : 


8  Wei  in  Series,  RS  =  0.980  (Hi)  (Cent...) 

min  X  =  0.001  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC=56 

2*NFC 

(112) 

0.1 

0.987 

0.800 

0.2 

0.981 

0.779 

2*(NFC+ 

NCOMP) 

(128) 

0.1 

0.985 

0.839 

0.2 

0.978 

0.824 

2*NFC- 

NCOMP 

(104) 

0.1 

0.988 

0.776 

0.2 

0.982 

0.753 

2*(NFC- 

NCOMP) 

(96) 

0.1 

0.989 

0.746 

0.2 

0.983 

0.732 

5 

Test  15  until 

3  failed. 

NFC  =  24 

2*NFC 

(48) 

0.1 

0.994 

0.621 

0.2 

0.991 

0.584 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.993 

0.705 

0.2 

0.988 

0.685 

2*NFC- 

NCOMP 

(40) 

0.1 

0.995 

0.548 

0.2 

0.992 

0.514 

0.1 

0.9% 

0.468 

0.2 

0.993 

0.417 

Table  3A  displays  the  results  of  Case  3  for  Type  A  systems.  In  this  case  8 
components  are  connected  in  series.  Four  of  them  have  failure  times  with 
exponential  distributions  and  the  remaining  four  have  failure  times  with  Weibull 
distributions  each  with  shape  parameter,  equal  to  2. 
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Table  3A  :  4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.980  (Hi) 

min  X  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


Test 

Deg  of 

S/N 

Plan 

Freedom 

Measures  of  Accuracy 


Test  5  until 
5  failed. 

NFC=40 


Test  15  until 
15  failed. 

NFC  =120 


Test  15  until 
11  failed. 

NFC=88 


2*NFC 

(80) 


2*(NFC+ 

NCOMP) 

(96) 


2*NFC- 

NCOMP 

(72) 


2*(NFC- 

NCOMP) 

(64) 


2*NFC 

(240) 


2*(NFC+ 

NCOMP) 

(256) 


2*NFC- 

NCOMP 

(232) 


2*(NFC- 

NCOMP) 

(224) 


2*NFC 

(176) 


2*(NFC  + 
NCOMP) 
(192) 


2*NFC- 

NCOMP 

(168) 


a 

RSLOW 

LEVEL 

0.1 

0.979 

0.942 

0.2 

0.978 

0.905 

0.1 

0.975 

0.987 

0.2 

0.974 

0.976 

0.1 

0.981 

0.881 

0.2 

0.980 

0.805 

0.1 

0.983 

0.771 

0.2 

0.982 

0.684 

0.1 

0.981 

0.863 

0.2 

0.980 

0.800 

0.1 

0.979 

0.941 

0.2 

0.979 

0.898 

0.1 

0.981 

0.881 

0.2 

0.980 

0.805 

0.1 

0.982 

0.725 

0.2 

0.981 

0.631 

0.1 

0.981 

0.864 

0.2 

0.980 

0.801 

0.1 

0.979 

0.951 

0.2 

0.978 

0.907 

0.1 

0.982 

0.802 

0.2 

0.981 

0.698 

0.1 

0.982 

0.702 

0.2 

0.982 

0.591 
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Table  3A  : 


4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.980  (Hi)  (Com...) 
min  \  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy  || 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC=56 

2*NFC 

(112) 

0.1 

0.981 

0.865 

0.2 

0.980 

0.787 

2*(NFC+ 

NCOMP) 

(128) 

0.1 

0.978 

0.952 

0.2 

0.978 

0.920 

2*NFC- 

NCOMP 

(104) 

0.1 

0.982 

0.769 

0.2 

0.982 

0.676 

2*(NFC- 

NCOMP) 

(96) 

0.1 

0.983 

0.644 

0.2 

0.983 

0.523 

5 

Test  15  until 

3  failed. 

NFC =24 

2*NFC 

(48) 

0.1 

0.982 

0.843 

0.2 

0.981 

0.762 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.976 

0.970 

0.2 

0.975 

0.941 

2*NFC- 

NCOMP 

(40) 

0.1 

0.984 

0.684 

0.2 

0.984 

0.580 

0.1 

0.987 

0.459 

0.2 

0.987 

0.356 

Inspection  of  Table  3A  reveals  the  following: 

(1)  The  two  procedures  corresponding  to  degrees  of  freedom  equal  to  2*NFC 
and  2*NFC-NCOMP  are  reasonably  accurate  for  all  5  simulation  cases. 

(2)  The  procedures  appears  to  be  nearly  equally  accurate  for  both  80%  and 
90%  confidence  levels. 


24 


B.  Test  Plan  2  :  Testing  for  a  Specified  Planned  Test  Time  (RETP2) 

In  the  simulations  for  test  plan  2,  components  were  tested  until  failure  or  until 
some  plaimed  test  time  scenario.  Failed  items  were  not  replaced.  Components 
whose  failure  times  had  exponential  distributions  were  tested  until  a  pre-determined 
total  test  time  was  accumulated  for  their  type  of  component.  Components  whose 
failure  times  had  Weibull  distributions  were  tested  until  failure  or  a  pre-determined 
planned  test  time  for  that  test  item.  The  latter  truncation  plan  is  needed  for  Weibull- 
type  items  in  order  to  use  the  maximum  likelihood  estimates  [as  in  equation  (2.17)] 
to  solve  for  ^  . 

Inspection  of  Table  4A  reveals  that  the  lower  confidence  limit  procedure  for 
degrees  of  freedom  equal  to  2*(1+NFC)  is  quite  accurate  when  enough  testing  is 
done  to  make  the  expected  number  of  failures,  E[NFC],  greater  than  or  equal  to  4.8. 
This  testing  constraint  is  well  within  the  domain  of  constraints  set  on  testing  in 
development  programs  for  major  systems  within  the  Department  of  Defense. 

Examination  of  Table  5  A  reveals  that  the  lower  confidence  limit  procedure  for 
degrees  of  freedom  equal  to  2*(1  +  NFC)  is  moderately  accurate  when  enough  testing 
is  done  to  make  E[NFC]  greater  than  or  equal  to  9.  The  accuracy  diminishes  slightly 
as  E[NFC]  increases.  This  could  be  corrected  by  decreasing  the  degrees  of  freedom 
slightly  to  make  RSLOW  slightly  larger. 

The  results  displayed  in  Table  6A  show  that  the  lower  confidence  limit 
procedure  for  degrees  of  freedom  equal  to  2*(1  +  NFC)  is  quite  accurate  when 
enough  testing  is  done  to  make  E[NFC]  greater  than  or  equal  to  9.6. 
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Table  4A  : 


8  Exp  in  Series,  RS  =  0.961  (Hi) 
X  =  0.001  f/hr,  UT  =  5  hrs 


S/N 

Degrees 

of 

Freedom 

K  /  E[NFC] 
(TT) 

a 

Measures  of  Accuracy  |j 

RSLOW 

LEVEL 

1 

2*(1+NFC) 

0.25  /  1.2 
(225) 

0.1 

0.950 

0.851 

0.2 

0.935 

0.851 

0.5  /  2.4 
(450) 

0.1 

0.957 

0.857 

0.2 

0.954 

0.857 

1.0  /  4.8 
(900) 

0.1 

0.957 

0.941 

0.2 

0.957 

0.850 

2.0  /  9.6 
(1800) 

0.1 

0.958 

0.916 

0.2 

0.960 

0.850 

3.0  /  14.4 
(2700) 

0.1 

0.959 

0.916 

0.2 

0.959 

0.809 

4.0  /  19.2 
(3600) 

0.1 

0.959 

0.937 

0.2 

0.960 

0.843 

5.0  /  24 
(4500) 

0.1 

0.960 

0.926 

0.2 

0.960 

0.814 

10.0  /  48 
(9000) 

0.1 

0.960 

0.924 

0.2 

0.960 

0.809 

20.0  /% 
(18000) 

0.1 

0.960 

0.914 

02 

0.961 

0.820 

30.0  /  144 
(27000) 

0.1 

0.%1 

0.906 

0.2 

0.%1 

0.804 
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Table  5A  : 


8  Wei  in  Series,  RS  =  0.956  (Hi)  (*) 
A  =  0.005  f/hr,  UT  =  15  hrs 


S/N 

Degrees 

of 

Freedom 

K  /  E(NFC1 
(TT) 

a 

Measures  of  Accuracy  | 

RSLOW 

LEVEL 

1 

2*(1+NFC) 

0.25  /  1.2 
(45) 

0.1 

1.000 

0.186 

0.2 

1.000 

0.158 

0.5  /  2.4 
(90) 

0.1 

0.986 

0.501 

0.2 

0.979 

0.458 

1.0  /  4.8 
(180) 

0.1 

0.%7 

0.767  • 

0.2 

0.960 

0.732 

2.0  /  9.6 
(360) 

0.1 

0.957 

0.879 

0.2 

0.952 

0.854 

3.0  /  14.4 
(540) 

0.1 

0.952 

0.934 

0.2 

0.946 

0.922 

4.0  /  19.2 
(720) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

10.0  /  48 
(1800) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

20.0/% 

(3600) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

30.0  /  144 
(5400) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

(•)  20  test  items  for  each  Weibull  component. 


Table  6A  :  4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.958  (Hi)  (•) 

A(exp)  =  0.001  f/hr,  UT(exp)  =  5  hrs 
X(wei)  =  0.005  f/hr,  UT(wei)  =  15  hrs 


S/N 

Degrees 

of 

Freedom 

Measures  of  Accuracy  | 

a 

RSLOW 

LEVEL 

1 

2*(1+NFC) 

0.25  /  1.2 
(225) 

(45) 

0.1 

1.000 

0.620 

0.2 

0.995 

0.451 

0.5  /  2.4 
(450) 

(90) 

0.1 

0.982 

0.623 

0.2 

0.975 

0.573 

1.0  /  4.8 
(900) 

(180) 

0.1 

0.971 

0.736 

0.2 

0.965 

0.684 

2.0  /  9.6 
(1800) 

(360) 

0.1 

0.964 

0.803 

0.2 

0.960 

0.765 

3.0  /  14.4 
(2700) 

(540) 

0.1 

0.960 

0.874 

0.2 

0.956 

0.841 

4.0  /  19.2 
(3600) 

(720) 

0.1 

0.960 

0.873 

0.2 

0.957 

0.839 

5.0  /  24 
(4500) 

(900) 

0.1 

0.959 

0.887 

0.2 

0.956 

0.861 

10.0  /  48 
(9000) 

(1800) 

0.1 

0.959 

0.891 

0.2 

0.956 

0.862 

20.0  /% 
(18000) 
(3600) 

0.1 

0.959 

0.892 

0.2 

0.955 

0.867 

0.1 

0.960 

mm 

0.2 

0.956 

0.862  1 
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The  accuracy  results  of  simulations  performed  in  this  thesis  cannot  be  extended 
to  systems  that  differ  significantly  from  those  simulated  here.  However,  it  can  be 
said  that  the  procedures  which  are  accurate  for  series  systems  are  also  usually 
accurate  for  1-out-of-it  parallel  systems  because  system  reliability,  can  be  written 
in  terms  of  component  reliabilities,  as 

R,  =  1  -  n  (1  -  «,) 

i-1 

k 

Thus,  upper  confidence  limits  on  J][  (1  -  Ri)  will  yield  a  lower  confidence 

1-1 

limit  on  R^.  The  accuracy  of  the  upper  confidence  interval  procedures  for 

k 

( 1  -  ,  a  series-type  problem,  should  be  nearly  the  same  as  those  obtained 

i=l 

in  this  thesis  because  equations  like  (2.9)  would  be  replaced  with  equations  for  the 
lower  confidence  limit  on  and  would  look  like 


X 


m,L(a) 


X  l-.2f 

i-1 


...  (4.1) 


If  the  degrees  of  freedom  parameter  F  is  large,  as  it  is  in  the  cases  simulated  in  this 
thesis,  the  associated  Chi-Square  distribution  with  2F  degrees  of  freedom  is  nearly 
symmetric  about  its  mean  so  the  lower  tail  has  nearly  the  same  shape  as  the  upper 
tail.  This  characteristic  should  yield  very  similar  accuracies  for  as  was 

obtained  for  X^  in  this  thesis.  These  accuracy  comparisons  translate  directly 

to  similar  comparisons  about  the  accuracies  of  the  associated  lower  confidence  limit 
procedures  for  the  system  reliability  of  series  and  parallel  systems. 
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Table  7  :  Summary  of  Procedure  Accuracy  by  Simulation  Category  (RETPl) 


S/N 

Key  Parameters  and 
their  Levels 

Observation  and  Discussion 
based  on  RSLOW  and  LEVEL 

1 

Run  Cases 

■  All  Exponential  (1) 

■  All  Weibull  (2) 

■  Mixed  (3) 

■  For  all  exponential  systems,  accurate 
procedures  were  developed  for  component 
sample  sizes  >  5. 

■  For  all  Weibull  systems,  component 
sample  sizes  should  be  >  15  with 
truncation  at  r  >  7  failures. 

2 

System  Reliabilities 

■  Hi  (>  0.9)  (Type  A) 

■  Lo  (>  0.8)  (Type  B) 

Accurate  procedures  were  developed  for 
both  cases  if  sample  sizes  are  adequate. 

3 

Levels  of  Significance 

■  a  =  0.1 

■  a  =  0.2 

Accuracy  varied  slightly  depending  on 
system  type  and  test  plan,  but  accurate 
procedures  exist  for  both  levels. 

4 

Degrees  of  Freedom 
.  DF  =  2*NFC 

■  DF  =  2*(NFC+NCOMP) 

■  DF  =  2*(NFC-NCOMP) 

■  DF  =  2*NFC-NCOMP 

Greatest  accuracy  as  follows; 

■  All  exponential  :  2*NFC-NCOMP 

■  All  Weibull  :  2*NFC-NCOMP 

or  2*(NFC-NCOMP) 

5 

Test  Plan 

■  Test  5  until  5  failures 

■  Test  15  until  15  failures 

■  Test  15  until  11  failures 

■  Test  15  until  7  failures 

■  Test  15  until  3  failures 

Accurate  procedures  existed  for  all  run 
cases  for  all  system  types  except  for  the  all 
Weibull  system  where  number  of  test 
items  n  >  15  and  number  of  failures  r 
should  be  >  7. 
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Table  8  :  Summary  of  Procedure  Accuracy  by  Simulation  Category  (RETP2) 


S/N 

Key  Parameters  and 
their  Levels 

Observation  and  Discussion 
based  on  RSLOW  and  LEVEL 

1 

Run  Cases 

■  All  Exponential  (4) 

■  All  Weibull  (5) 

■  Mixed  (6) 

Procedures  were  accurate  for  DF  = 

2*(1  +  NFC)  when  enough  testing  was 
done  to  make  the  expected  number  of 
failed  components  e[nFC]  >  9. 

2 

System  Reliabilities 

■  Hi  ( >  0.9)  (Type  A) 

■  Lo  (>  0.8)  (Type  B) 

Same  as  S/N(l). 

3 

Levels  of  Significance 

■  a  =  0.1 

■  tt  =  0.2 

Same  as  S/N(l). 

4 

Degrees  of  Freedom 

■  DF  =  2*(1  +  NFC) 

■  DF  =  1.3*2*(1  +  NFC) 

DF  =  2*(1  +  NFC)  was  the  most  accurate 
procedure. 

5 

Test  Plan 
■  K  factors  of 

0.25,  0.5,  1,  2,  3, 

4,  5,  10,  20,  30 

K  should  be  chosen  so  that  E(NFC]  >  9. 

Tables  7  and  8  provide  cursory  summaries  of  some  constraints  needed  to  assure 
the  existence  of  one  or  more  accurate  lower  confidence  limit  procedures  among  the 
procedures  that  were  evaluated.  The  simulation  scenarios  are  divided  into  five 
categories  for  this  summarization. 
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V.  APPLICATION  EXAMPLES 

Based  on  the  procedures  evaluated  by  the  RETPl  and  RETP2  runs,  four 
different  test  plans  and  failure  time  data  were  constructed  to  illustrate  the  use  of  the 
procedures  in  providing  a  lower  lOO(l-a)  %  confidence  limit  for  the  system  reliability 
of  a  series  system  with  different  types  of  components. 


CASE  1  :  8  Exponential  Components  In  Series 

-  TEST  PLAN  1  -  Test  15  until  7  falls  for  each  component 

I .  Raw  Data 


Comp 

i 

T(l) 

T(2) 

Ordered 

T(3) 

Failure 

T(4) 

Times  (h) 
T(5) 

T(6) 

T(7) 

1 

300.0 

*00.0 

500.0 

600.0 

700.0 

800.0 

900.0 

2 

350.0 

*50.0 

550.0 

650.0 

750.0 

850.0 

950.0 

3 

AOO.O 

500.0 

600.0 

700.0 

800.0 

900.0 

1000.0 

k 

*50.0 

550.0 

650.0 

750.0 

850.0 

950.0 

1050.0 

5 

500.0 

600.0 

700.0 

800.0 

900.0 

1000.0 

1100.0 

6 

550.0 

650.0 

750.0 

850.0 

950.0 

1050.0 

1150.0 

7 

600.0 

700.0 

800.0 

900  0 

1000.0 

1100.0 

1200.0 

8 

650.0 

750.0 

850.0 

950.0 

1050.0 

1150.0 

1250.0 

II .  Data  Sunxnary 


Comp  ER(1)* 


i 

UT(1) 

NC(l) 

KF(1) 

TT(1) 

ELM(i) 

ER(i) 

TT(i) 

1 

5.0 

15 

7 

11*00.0 

0.00053 

1.00000 

11*00.0 

2 

5.0 

15 

7 

12150.0 

0.000*9 

0.93827 

11*00.0 

3 

5.0 

15 

7 

12900.0 

0.000*7 

0.88372 

11*00.0 

k 

5.0 

15 

7 

13650.0 

0.000** 

0.83516 

11*00.0 

5 

5.0 

15 

7 

1**00. 0 

0.000*2 

0.79167 

11*00.0 

6 

5.0 

15 

7 

15150.0 

0.000*0 

0.752*8 

11*00.0 

7 

5.0 

15 

7 

15900.0 

0.00038 

0.71698 

11*00.0 

8 

5.0 

15 

7 

16650.0 

0.00036 

0.68*68 

11*00.0 

III.  Estimation  Procedure  for  RSLOW 


Parameter 

df 

Value 

ALPHA 

0.1 

NFC 

56 

CHISQD 

LMU 

112  131.56  (df  -  2  •  NFC  ,  CHISCJD  from  tables) 

0.00072 

RSLOW 

0.976*7 
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CASE  2  :  8  Weibull  Components  in  Series 

-  TESX  plan  1  -  Test  15  until  7  fails  for  each  component 


I .  Raw  Data 


Comp 

i 

T(l) 

T(2) 

Ordered 

T(3) 

Failure 

T(l) 

Times  (h) 
T(5) 

T(6) 

T(7) 

1 

10.0 

20.0 

30.0 

10.0 

50.0 

60.0 

70.0 

2 

20.0 

30.0 

10.0 

50.0 

60.0 

70.0 

80.0 

3 

30.0 

10.0 

50.0 

60.0 

70.0 

80.0 

90.0 

4 

10.0 

50.0 

60.0 

70.0 

80.0 

90.0 

100.0 

5 

50.0 

60.0 

70.0 

80.0 

90.0 

100.0 

110.0 

6 

60.0 

70.0 

80.0 

90.0 

100.0 

110.0 

120.0 

7 

70.0 

80.0 

90.0 

100.0 

110.0 

120.0 

130.0 

8 

80.0 

90.0 

100.0 

110.0 

120.0 

130.0 

110.0 

II.  Data 

Comp 

i 

Summary 

UT(i) 

NC<1) 

NF{i) 

TT(i) 

ELM(i) 

ER(i) 

ER(i)* 

TT(i) 

1 

5.0 

15 

7 

5.3E+01 

1.32E-01 

1.00000 

5.3E+01 

2 

5.0 

15 

7 

7 . 2E+01 

9.79E-05 

0.71106 

5.3E+01 

3 

5.0 

15 

7 

9.3E+01 

7.51E-05 

0.57328 

5.3E+01 

4 

5.0 

15 

7 

1.2E+05 

5.98E-05 

0.15131 

5.3E+01 

5 

5.0 

15 

7 

l.lE+05 

1.85E-05 

0.36812 

5.3E+01 

6 

5.0 

15 

7 

1.7E+05 

l.OlE-05 

0.30152 

5.3E+01 

7 

5.0 

15 

7 

2.1E+05 

3.37E-05 

0.25577 

5.3E+01 

8 

5.0 

15 

7 

2.1E+05 

2.87E*05 

0.21777 

5.3E+01 

III.  Estimation  Procedure  for  RSLOW 


Parameter 

df 

Value 

BETA 

2.0 

ALPHA 

0.1 

NFC 

56 

CHISQD  112  131.56  (df  -  2  *  NFC  ,  CHISQD  from  tables) 

LMU  1.55E-01 

RSLOW  0.98«97 


IV.  Workarea 


Cocnp 

1 

Ordered  Failure  Times  (h)  raised  to  the  power  of  BETA 
T'd)  T'(2)  T’(3)  T'(l)  T'(5)  T’(6)  T'(7) 

1 

l.OE+02 

l.OE+02 

9.0E+02 

1 . 6E+03 

2 . 5E+03 

3.6E+03 

1.9E+03 

2 

l.OE+02 

9.0E+02 

1.6E+03 

2 . 5E+03 

3.6E+03 

1.9E+03 

6.1E+03 

3 

9.0E+02 

1.6E+03 

2.5E+03 

3.6E+03 

1 . 9E+03 

6.1E+03 

8.1E+03 

4 

1.6E+03 

2 . 5E+03 

3.6E+03 

1.9E+03 

6.1E+03 

e.lE+03 

l.OE+01 

5 

2.5E+03 

3.6E+03 

1 . 9E+03 

6.1E+03 

8.1E-t03 

l.OE+01 

1 . 2E+01 

6 

3.6E+03 

1.9E+03 

6.1E+03 

8.1E+03 

1 . OE+01 

1.2E+01 

l.lE+01 

7 

1.9E+03 

6.1E+03 

8. lE+03 

1 . OE+01 

1.2E+01 

l.lE+01 

1 . 7E+01 

8 

6.1E+03 

8. lE+03 

l.OE+01 

1 . 2E+01 

l.lE+01 

1.7E+01 

2.0E+01 
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CASE  3  ;  8  Exponential  Components  in  Series 
-  TEST  PLAN  2  “  Test  until  TT(i)  for  each  component 

I .  Raw  Data 


No  data  except  noting  the  number  of  failures  for  each  component  (NF(i)) 

II.  Data  Summary 


Coop 

1 

UT<i) 

NF(1) 

TT(i) 

ELM(i) 

ER(i) 

ER(i)* 

TT(i) 

ER(i)* 

UT(i) 

1 

5.0 

6 

2000.0 

0.00300 

0.85714 

1714.3 

4.3 

2 

5.0 

6 

2000.0 

0.00300 

0.85714 

1714.3 

4.3 

3 

5.0 

5 

2000.0 

0.00250 

0.71429 

1428.6 

3.6 

4 

5.0 

7 

2000.0 

0.00350 

1.00000 

2000.0 

5.0 

5 

5.0 

5 

2000.0 

0 . 00250 

0.71429 

1428.6 

3.6 

6 

5.0 

5 

2000.0 

0.00250 

0.71429 

1428.6 

3.6 

7 

5.0 

4 

2000.0 

0.00200 

0.57143 

1142.9 

2.9 

8 

5.0 

7 

2000.0 

0.00350 

1.00000 

2000.0 

5.0 

III.  Estimation  Procedure  for  RSLOU 


Parameter  df  Value 


ALPHA  0 . 1 

NFC  45 

CHISQD  92  109.76  (df 

LMU  0.00427 


2  *  (  1  +  NFC  )  ,  CHISQD  from  tables) 


RSLOU 


0.87180 


CASE  4  :  8  Welbull  Components  in  Series 

-  TEST  PLAN  2  -  Test  20  until  TT(i)  for  each  component 


I .  Raw  Data 


Comp 

i 

TT(i) 

T(l) 

Ordered 

T(2) 

Failure 

T(3) 

Times  (h) 
T{4) 

T(5) 

T(6) 

1 

60.0 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

2 

60.0 

15.0 

25.0 

35.0 

45.0 

55.0 

3 

60.0 

20.0 

30.0 

40.0 

50.0 

60.0 

4 

60.0 

25.0 

35.0 

45.0 

55.0 

5 

60.0 

30.0 

40.0 

50.0 

60.0 

6 

60.0 

35.0 

45.0 

55.0 

7 

60.0 

40.0 

50.0 

60.0 

8 

60.0 

45.0 

55.0 

II.  Data  Suimary 


Comp 

i 

UT(i) 

NC(i) 

NF(i> 

ET(i) 

ELM{i) 

ER(i) 

ER(i>* 

ET(i> 

1 

5.0 

20 

6 

6 . OE+04 

l.OlE-04 

1.00000 

6. OE+04 

2 

5.0 

20 

5 

6.  lE-t-04 

8.18E-05 

0.81118 

5.0E+04 

3 

5.0 

20 

5 

6 . 3E+04 

7.94E-05 

0.78704 

5.0E+04 

4 

5.0 

20 

4 

6 . 5E+04 

6.20E-05 

0.61499 

4 . OE+04 

5 

5.0 

20 

4 

6.6E+04 

6.04E-05 

0.59919 

4 . OE+04 

6 

5.0 

20 

3 

6.7E+04 

4.45E-05 

0.44090 

3 . OE+04 

7 

5.0 

20 

3 

6.9E+04 

4.35E-05 

0.43179 

3 . OE+04 

8 

5.0 

20 

2 

7.0E+04 

2.86E-05 

0.28394 

2.0E+04 

III.  Estimation  Procedure  for  RSLOW 


Parameter 

df 

Value 

BETA 

2.0 

ALPHA 

0.1 

NFC 

32 

CHISQD  66  81.09  (df  -  2  *  (  1  +  NFC  )  ,  CHISQD  from  tables) 

LMU  1.28E-04 

RSLOW  0.9842S 

IV.  Workarea 


Comp  Ordered  Failure  Times  (h)  raised  to  the  power  of  BETA 


i 

TT'(l) 

T'(l) 

T'(2) 

T’(3) 

T’(4) 

T'(5) 

T'(6) 

1 

3.6E+03 

l.OE+02 

4 . OE+02 

9 . OE+02 

1.6E+03 

2.5E+03 

3 . 6E+03 

2 

3 . 6E+03 

2.3E+02 

6 . 3E+02 

1 . 2E+03 

2.0E+03 

3 . OE+03 

3 

3.6E+03 

4 . OE+02 

9.0E+02 

1.6E+03 

2 . 5E+03 

3.6E+03 

4 

3.6E+03 

6 . 3E+02 

1.2E+03 

2 . OE+03 

3.0E+03 

5 

3.6E+03 

9 . OE+02 

1.6E+03 

2 . 5E+03 

3.6E+03 

6 

3.6E+03 

1.2E+03 

2.0E+03 

3.0E+03 

7 

3.6E+03 

1.6E+03 

2.5E+03 

3.6E+03 

6 

3.6E+03 

2.0E+03 

3.0E+03 
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VI.  CONCLUSION 


Some  of  the  lower  confidence  limit  procedures  developed  and  evaluated  in  this 
thesis  are  reasonably  accurate  for  the  series  systems  simulated  for  test  plans  with 
sample  sizes  and  truncation  scenarios  that  are  usually  experienced  in  DoD  aquisition 
programs.  The  accuracy  of  these  methods  can  be  varied  by  modifying  the  degrees 
of  freedom  parameter,  F,  in  the  x^«,f  term  in  equation  (2,9),  namely: 


_ 

m.t/(o)  - 1 - 

i-1 


The  computer  program  can  be  modified  with  modest  effort  to  accomodate  specific 
complex  coherent  systems  so  long  as  the  components  have  failure  time  distributions 
that  are  exponential  or  Weibull.  This  means  that  the  computer  program  provided  in 
this  thesis  can  be  used  to  develop  a  reasonably  accurate  lower  confidence  limit  for 
the  system  reliability  of  a  specific  complex  quasi-coherent  system  with  independent 
components.  This  can  be  done  by  choosing  the  failure  distribution  and  associated 
parameters  of  the  components,  the  corresponding  test  plan  parameters  and  the 
desired  level  of  confidence.  The  simulation  can  then  be  run  for  this  set  of 
parameters  for  various  equations  for  the  degrees  of  freedom  parameter,  F,  to 
determine  an  equation  for  F  that  yields  a  lower  confidence  limit  with  a  satisfactory 
degree  of  accuracy. 

When  testing  is  truncated  on  the  number  of  failures,  a  reasonably  accurate 
procedure  existed  for  all  3  cases  of  systems  (all  exponential,  all  Weibull  and  mbced) 
that  were  simulated  when  (a)  the  sample  size  of  the  components  was  10  or  larger, 
and,  (b)  the  ratio  of  the  number  of  failures  to  the  sample  size  was  at  least  0.5.  When 
testing  was  truncated  by  planned  test  time,  reasonably  accurate  procedures  were 
found  for  cases  where  the  expected  number  of  failures  was  at  least  7. 
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VII.  RECOMMENDATIONS 

The  computer  program  developed  in  this  thesis  facilitates  the  development  of 
lower  confidence  limit  procedures  for  explicit  quasi-coherent  systems.  Systems  other 
than  series  systems  with  large  numbers  of  components  (eg.  30)  should  be  simulated 
to  test  the  versatility  of  the  general  lower  confidence  limit  methods  used  here. 

Modified  estimates  for  the  0  (shape)  parameter  in  the  Weibull  failure  time 
distribution  and  the  parameter  r  =  should  be  explored  in  an  attempt  to  find 
more  accurate  procedures. 
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APPENDIX  A  :  Derivation  of  Formula  Used 


Suppose  T  has  an  exponential  distribution  with  failure  rate  Suppose  '^(ly  '^(2y 
are  the  first  r  ordered  statistics  in  a  random  sample  of  size  n  from  this 
exponential  distribution.  Let  5  be  defined  by 

=  E  ^(0  (« 

j-i 

It  is  well  known  that  2\S  has  a  Chi-Square  distribution  with  2r  degrees  of  freedom 
(See  Ref.3  p  488).  The  maximum  likelihood  estimator  for  X  is  given  by 

i  =  L 
S 

We  seek  an  unbiased  estimator  for  X  .  Suppose  X  has  a  Chi-Square  distribution 
with  2r  degrees  of  freedom,  then 

and  the  integral  of  this  function  from  zero  to  infinity  equals  1  (See  Ref.3). 


«• 

£(— )  =  f — i — x^'"^^’*exp(-— )4ic 

J2T(r)  2' 

=  r - 1 - A:^'^"’^‘^exp(-— )tic 

2T(r)]z-^T{r-l)  2' 


2T{r) 

1 

2(r-l) 
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Then 


=  £(1)  .  r2XE(^)  = 


rX 

Fi) 


Therefore 


E(!^X)  =  £(^)  =  X 

r  o 
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APPENDIX  B  ;  Users’  Guide  for  RETPl 


Reliability  Estimation  Test  Plan  1  (RETPl). 
by  YEE,  Kah-Chee 
July  91 


1.  Brief  Description. 

RETPl  is  a  computer  program  written  in  FORTRAN  that  runs  on  the  Amdahl 
mainframe  at  NPGS.  It  allows  the  user  to  simulate  exponential  and  Weibull  failure 
times  of  component  items  being  tested  to  evaluate  the  accuracy  of  a  confidence  limit 
estimation  procedure  based  on  Type  II  data  censoring  (that  is,  testing  n^  items  of 
component  i  until  fj  of  them  failed). 

2.  Program  Input.  (INl.DATA) 

The  input  of  the  program  are  specified  to  the  program  via  an  input  file  called 
INl.DAT.  A  sample  input  file  is  shown  below. 

This  file  contains  the  inputs  required  by  the  RETPl  model. 

Update  only  the  numerical  values  between  dotted  lines  as  appropriate. 

Do  not  delete  any  of  the  comment  lines.  (INl.DAT)  14  May  91 


Value 

Type 

Units 

Description 

Variable 

16807.0 

Real 

• 

initial  random  seed 

ISEED 

8 

Int 

- 

total  #  of  components  in  system 

NCOMP 

4 

Int 

- 

#  of  Exponential  components 

NEXP 

4 

Int 

- 

y/  of  WEIbull  components 

NWEI 

0 

Int 

- 

#  of  GEOmetric  components 

NGEO 

0.01 

Real 

- 

tolerance  for  MLE 

TOL 

0.20 

REAL 

- 

DESIRED  SIGNIFICANCE  LEVEL 

ALPHA 

1000 

Int 

- 

y/  of  replications  desired 

NREP 

3 

Int 

test  case  number 

1  -  all  EXP 

2  -  all  WEI 

3  -  EXP  +  WEI 

4  -  EXP  +  WEI  +  GEO 

TCN 

8 

Int 

- 

number  of  cut  sets 

NCS 
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TEST  PLAN  :  Testing  NC(I)  items  of  component  i 

until  NF(I)  of  them  fails.  (Use  REAL  numbers  ONLY!!!) 


Comp 

Comp 

Comp  Parameters 

Util 

Test  Plan  Inputs 

Number 

Type 

Scale 

Shape 

Time/Cycle 

y/  Comp 

y/  Failed 

I 

TY(I) 

PARMl(I) 

PARM2(I) 

UT(I)  UC(I) 

NC(I) 

NF(I) 

Int 

Int 

Real 

Real 

(hrs)  Int 

Int 

Int 

1.0 

1.0 

0.0020 

1.0 

10.0 

15.0 

3.0 

2.0 

1.0 

0 . OOAO 

1.0 

10.0 

15.0 

3.0 

3.0 

1.0 

0.0060 

1.0 

10.0 

15.0 

3.0 

A.O 

1.0 

0.0080 

1.0 

10.0 

15.0 

3.0 

5.0 

2.0 

0.0020 

2.0 

10.0 

15.0 

3.0 

6.0 

2.0 

O.OOAO 

2.0 

10.0 

15.0 

3.0 

7.0 

2.0 

0.0060 

2.0 

10.0 

15.0 

3.0 

8.0 

2.0 

0.0080 

2.0 

10.0 

15.0 

3.0 

Note  : 

TY(I) 

-1  EXPONENTIAL  P(surv)  -  exp( -PARM1)*T) 

TY(I) 

-2  WEIBULL  P(surv)  -  exp( - (PARM1*T)**PARM2) 

TY(I) 

-3  GEOMETRIC  P(surv)  -  PARM1**N 

SYSTEM 

CONFIGURATION  : 

Identification  of  CUT  SETs 

-  min  groups 

of  components 

that  have 

to  fail 

for  the  system  to  fail. 

Cut  Set 

y/  in  Set 

List  of  Components  in  Cutset 

J 

C0MP(J,1) 

C0MP(J,2)  . 

.  .  up  to  C0MP(J,1)  components 

1 

1 

10  0  0 

0  0  0  0  0 

2 

1 

2  0  0  0 

0  0  0  0  0 

3 

1 

3  0  0  0 

0  0  0  0  0 

A 

1 

A  0  0  0 

0  0  0  0  0 

5 

1 

5  0  0  0 

0  0  0  0  0 

6 

1 

6  0  0  0 

0  0  0  0  0 

7 

1 

7  0  0  0 

0  0  0  0  0 

8 

1 

8  0  0  0 

0  0  0  0  0 

3.  Program  Flow  and  Lx)gic.  (NAMEl.DEF,  PARMl.DEF  and  RETPl.FOR) 

Input  parameters  are  first  read  in  by  the  program  by  calling  the  INPUT 
subroutine.  The  program  then  evoke  the  SIM  subroutine  which  generates  the 
random  failure  times  and  compute  the  key  statistics  required  in  the  procedure.  The 
next  subroutine  EVAL  determines  the  measures  of  accuracy  for  run.  REPORT  is 
the  subroutine  which  generates  the  output  file  for  the  run  OUTl.DAT. 


The  variables  in  the  program  RETPl.FOR  are  described  in  the  file 
NAMEl.DEF  as  listed  below. 
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c . 

C  This  file  contains  the  declaration  for  input  and  output  variables 
C  used  in  the  the  RETPl  model.  (NAMEl.DEF)  14  May  91 

C . 

C  Input  Variables. 

C  . 

C  ISEED  -  initial  random  seed  selected. 

C  SEED  -  current  random  seed. 

C  RS  -  true  overall  series  system  reliability. 

C  ALPHA  -  level  of  significance  desired. 

C  NREP  -  number  of  replications  desired  for  the  simulation. 

C  TPN  -  test  plan  number  (1). 

C  TCN  -  test  case  number  (1,  2,  3  or  4) . 

C  NCOMP  -  total  number  of  components  in  the  system. 

C  NEXP  “  number  of  components  with  EXP  failure  times. 

C  NWEI  -  number  of  components  with  WEI  failure  times . 

C  NGEO  -  number  of  components  with  GEO  failure  times. 

C  TOL  —  desired  tolerance  for  MLE  of  WEI  shape  parameter. 


c 

Distribution: 

Exponential 

WEIbull 

GEOmetric 

c 

TY(I)  -  type: 

1 

2 

3 

c 

PARMd.I)  : 

Scale(l/hr) 

Scale(l/hr) 

Prob 

G 

PARM(2.I) ; 

- 

Shape 

- 

C  UT(I)  —  utilization  time  (hrs)  for  component  i  (EXP  and  WEI). 

C  UC(I)  -  utilization  cycles  for  component  i  (GEO  only). 

C  NC(I)  -  number  of  test  samples  (sample  size)  for  emponent  i. 

C  NF(I)  -  desired  number  of  failures  in  test  for  component  i. 

C  NCS  -  number  of  cut- sets  for  the  system. 

C  COMP(J,K)  -  kth  parameter  of  cut-set  j  (first  being  the  no.  of 
C  components  belonging  to  the  cut-set). 

C 

C  Assumed  Variables. 

C  . 

C  MAXCOMP  “  maximum  number  of  components  allowed  in  the  system. 

G  MAXREP  -  maximum  number  of  replications  permitted. 

C  MAXCUT  -  maximum  number  of  cut-sets. 

C 

C  Program  and  Output  Variables. 

C  . 

C  RS  -  true  overall  system  reliability. 

C  TT(I)  -  total  accumulated  failure  time  (hr)  for  component  i 
C  (EXP  and  WEI  only) . 

C  TC(I)  -  total  accumulated  cycles  to  failure  (incl.  failure  cycle) 

C  for  component  i  (GEO  only) . 

C  EBETA(I)  -  estimate  for  shape  parameter  of  component  i  (if  Weibull) . 

G  RELl(J)  -  actual  reliability  for  cut-set  j. 

C  REL2(J)  -  computed  reliability  for  cut-set  j  for  current  replication. 
G  ELM(I)  -  estimated  component  failure  rate  (1/hrs)  for  component  i. 

C  ELMAX(M)  -  max  estimated  component  failure  rate  for  rep  m  (1/hrs). 

C  ER(I)  -  ratio  of  estimated  failure  rate  to  LMAX. 

C  NFC(M)  -  total  number  of  failed  test  components. 
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C  LMU(M)  -  upper  confidence  bound  for  failure  rate  (1/hrs), 

C  RSL(M)  -  lower  confidence  limit  estimated  for  system  reliability 
C  ...  for  the  mth  replication. 

C  ORSL(M)  -  ordered  RSL(M)  (ascending) . 

C  RSLOW  -  (l-ALPHA)xlOO  percentile  of  set  of  RSL(M) . 

C  LEVEL  -  achieved  confidence  level,  ie.  proportion  of  RSL(M)  that 
C  are  lesser  than  RS  (conservative  estimate) . 

C 

C--  END  OF  NAMEl.DEF . - . --- 

C 


Together  with  the  main  program  in  RETPl.FOR  are  the  other  subroutines 
needed  in  the  simulation.  The  declaration  of  variables  is  done  in  the  file 
PARMl.DEF.  Relevant  descriptions  are  included  as  comment  lines  in  the  source 
code  to  help  explain  the  program  segments.  A  listing  of  PARMl.DEF  and 
RETPl.FOR  is  given  below. 

c . 

C  This  file  contains  the  declaration  for  input  and  output  variables 
C  used  in  the  the  RETPl  model.  (PARMl.DEF)  14  May  91 

C . 

INTEGER  MAXCOMP,  MAXREP 

PARAMETER (  MAXCOMP  -  100  ,  MAXREP  -  1000  ,  MAXCUT  -  20  ) 

REAL*8  I SEED,  SEED 

INTEGER  NREP,  TCN,  NCOMP,  NEXP,  NWEI .  NGEO,  NCS , 

*  NC(MAXCOMP),  NF(MAXCOMP),  TY ( MAXCOMP ) ,  NFC (MAXREP) , 

*  UC(MAXCOMP),  TC(MAXCOMP),  COMP (MAXCUT .MAXCOMP) 

REAL*8  RS,  ALPHA,  UT ( MAXCOMP ) ,  TT (MAXCOMP ) , 

*  PARM (2, MAXCOMP ) ,  ELM ( MAXCOMP ) ,  ER (MAXCOMP ) , 

*  LMU(MAXREP),  RS L( MAXREP ) ,  ORSL(MAXREP) , 

*  ELMAX ( MAXREP ) ,  RSLOW,  LEVEL,  TOL,  EBETA (MAXCOMP ) , 

*  RELl (MAXCUT) ,  REL2 (MAXCUT) 

C 

COMMON/BLOCKl/ISEED,  SEED,  NREP.  TCN,  NCOMP,  NC,  NF,  NEXP,  NWEI , 

*  NGEO,  NCS,  TY,  NFC,  UC.  TC,  COMP 
COMMON/BLOCK2/RS ,  ALPHA.  UT,  TT.  PARM,  ELM,  ER,  LMU, 

*  RSL,  ORSL,  ELMAX.  RSLOW,  LEVEL,  TOL.  EBETA, 

*  RELl,  REL2 
C 

C--  END  OF  PARMl.DEF  . 

C 

C . 

C  This  file  contains  the  main  program  and  the  subroutines 
C  for  the  Reliability  Estimation  Test  Plan  1  (RETPl)  model. 

C  (RETPl.FOR)  -  runs  on  a  IBM  PC  Compatible. 

C 

C  MAINFRAME  VERSION. . . 

C 

C  Test  Plan  1  :  Testing  NC(I)  items  for  component  i 
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Until  NF(I)  of  them  fails. 


C 
C 

C  by  Yee  Kah-Chee  SMC  2802. 

C  14  May  91. 

C . . 

C  1.  Main  Program  (RETPl) . 

C . 

PROGRAM  RETPl 

Include  the  declaration  files. 

INCLUDE  'NAMEl  DEF' 

INCLUDE  'PARMl  DEF' 

Read  in  input  data. 

CALL  INPUT 

Activate  simulation. 

CALL  SIM 

Process  and  evaluate  output  data. 
CALL  EVAL 

Generate  simulation  report. 

CALL  REPORT 

STOP 
END 


2.  Input  Initialisation  Subroutine  (INPUT). 

SUBROUTINE  INPUT 

This  subroutine  reads  in  the  inputs  for  the  RETPl  model. 

Include  the  declaration  file. 

INCLUDE  'PARMl  DEF' 

INTEGER  I,  J,  K,  DUM2(11) 

REAL*8  DUM1(7) 

Read  data  from  'INI. DAT'  designated  as  logic  unit  1. 
OPEN(UNIT-l) 
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READd.lO) 

10  FORMAT( IX. /////) 

READ(1,*)  ISEED 
READd,*)  NCOMP 
READd,*)  NEXP 
READd,*)  NWEI 
READd,*)  NGEO 
READ(1,*)  TOL 
READ(1,*)  ALPHA 
READd,*)  NREP 
READd,*)  TCN 
READd, 20) 

20  FORMAT (IX,///) 

READd,*)  NCS 
C 

READd, 30) 

30  FORMAT (IX, ////////) 

C 

DO  50  K  -  1,  NCOMP 
READ(1,*)  DUMl 
C 

I  -  NINT(DUM1(1)) 

TY(I)  -  NINT(DUM1(2)) 

C 

IF  (TY(I).EQd)  THEN 
PARMd.I)  -  DUM1(3) 
PARM(2,I)  -  DUM1(4) 
UT(I)  -  DUM1(5) 

NC(I)  -  NINT(DUM1(6)) 
NF(I)  -  NINT(DUM1(7)) 
EBETA(I)  -  0 

ELSEIF  (TY(I).EQ.2)  THEN 
PARMd.I)  -  DUM1(3) 
PARM(2,I)  -  DUM1(4) 
UT(I)  -  DUM1(5) 

NC(I)  -  NINT(DUM1(6)) 
NF(I)  -  NINT(DUM1(7)) 
ELSEIF  (TYd).EQ.3)  THEN 
PARMd.I)  -  DUM1(3) 
PARM(2,I)  -  DUM1(4) 
UC(I)  -  NINT(DUM1(5)) 
UT(I)  -  DUM1(5) 

NC(I)  -  NINT(DUM1(6)) 
NF(I)  -  NINT(DUM1(7)) 
END  IF 

50  CONTINUE 
C 

READ(1,60) 

60  FORMAT (IX, ///////////) 

DO  80  I  -  1,  NCS 
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READd,*)  DUM2 
J  =  DUM2(1) 

C0MP(J,1)  -  DUM2(2) 

DO  70  K  -  1,  C0MP(J,1) 
C0MP(J,K+1)  -  DUM2(K+2) 
70  CONTINUE 
80  CONTINUE 

CLOSE (UNIT-1) 

RETURN 

END 


3.  Subroutine  (SIM) 


SUBROUTINE  SIM 

This  subroutine  simulates  NREP  possible  outcomes  of  the  test  plan 
desired  in  order  to  obtain  the  raw  estimates  of  LMU(M)  and  RSL(M) 
for  each  of  the  replication. 

Include  the  declaration  file 
and  declare  local  variables. 

INCLUDE  'PARMl  DEF' 

INTEGER  I,  J,  K,  M,  ISUM,  KEY 

REAL* 8  UNI 

REAL* 8  SUM,  PROD, 

*  FT(MAXCOMP),  OFT(MAXCOMP) 

SEED  -  ISEED 

Compute  overall  true  system  reliability  RS. 


RS  -  1.0 
DO  30  J  -  1,  NCS 
PROD  -  1.0 

DO  20  I  -  1,  C0MP(J,1) 

K  -  C0MP(J,I+1) 

PROD  -  PR0D*(  1  -  SURV(TY(K) ,PARM(1,K) ,PARM(2,K) ,UT(K))  ) 
20  CONTINUE 

RELl(J)  -  1.0  -  PROD 
RS  -  RS  *  RELl(J) 

30  CONTINUE 

Start  of  Simulation. 

C  (Initialize  replication  counter  M) . 

C 

M  -  1 

DO  WHILE  (M.LE.NREP) 

C 

C  PRINT  35,  M 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


35  FORMAT (IX, 'Replication  ',14) 

Test  Plan  ;  Sample  and  determine  unknown  TT(I) 

.  with  known  NC(I)  until  NF(I)  fails. 

Generate  NC(I)  failure  times,  put  them  in  ascending  order 
with  the  smallest  failure  time  on  the  top  of  the  list. 

DO  70  1  -  1,  NCOMP 

DO  40  K  -  1,  NC(I) 

CALL  LRNDPC( SEED, UNI, 1) 

IF(TY(I) .EQ.l)  THEN 

FT(K)  -  -L0G(UNI)/PARM(1,I) 

ELSEIF(TY(I) .EQ.2)  THEN 

FT(K)  -  (1.0/PARM(1,I))*(-LOG(UNI))**(1.0/PARM(2,I)) 
ELSEIF(TY(I) .EQ.3)  THEN 
FT(K)  =  1.0 

DO  WHILE  (UNI.LT.PARMd,!)) 

FT(K)  -  FT(K)  +  1.0 
CALL  LRNDPC( SEED, UNI, 1) 

ENDDO 

ENDIF 

40  CONTINUE 


Bubble  Sort  the  failure  times  in  ascending  order. 


CALL  BUBBLE(NC(I) ,FT,0FT) 


Compute  the  total  time  accumulated  in  the  test  and  the  estimate 
for  the  failure  rate  of  the  component  as  in  the  procedure. 


IF  (TY(I).NE.2)  THEN 
SUM  -  0.0 

DO  50  K  -  1,  NF(I) 

SUM  “  SUM  +  OFT(K) 

50  CONTINUE 

TT(I)  -  FL0AT(NC(I)-NF(I))*0FT(NF(I))  +  SUM 
IF  (TY(I) .EQ.l)  THEN 

ELM(I)  -  FL0AT(NF(I)-1)/TT(I) 

ELSEIF  (TY(I).EQ.3)  THEN 

ELM(I)  -  FL0AT(NF(I)-1)/TT(I) 

ENDIF 

ELSEIF  (TY(I).EQ.2)  THEN 
PRINT  55,  M,  I 

55  FORMATdX, 'REP  -  ',13,  '  Comp  =  ',13,/) 

C 

CALL  MLESHAPE(0FT,NC(I)  ,NFd ) , TOL,  1 . DO , EBETAd ) ) 
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C 

C 


EBETA(I)  -  BN(NC(I))*EBETA(I) 

SUM  -  0.0 

DO  60  K  -  1,  NF(I) 

SUM  -  SUM  +  OFT(K)**EBETA(I) 

CONTINUE 

TT(I)  -  FLOAT(NC(I)-NF(I))*OFT(NF(I))**EBETA(I)  +  SUM 
ELM(I)  -  FL0AT(NF(I))/TT(I) 

ENDIF 

CONTINUE 


Determine  the  total  number  of  failed  test  items. 
ISUM  -  0 

DO  80  I  -  1,  NCOMP 
ISUM  -  ISUM  +  NF(I) 

80  CONTINUE 

NFC(M)  -  ISUM 


Determine  the  maximum  failure  rate  estimate 
and  identify  that  component. 

ELMAX(M)  -  0.0 
KEY  -  0 

DO  90  I  -  1,  NCOMP 

IF  (ELM(I)  .GT.ElilAX<M))  THEN 
ELMAX(M)  =  ELM(I) 

KEY  -  I 
ENDIF 

90  CONTINUE 


Compute  the  ratios  of  the  failure  rate  estimates  to  their  maximum. 

DO  100  I  -  1,  NCOMP 

ER(I)  -  ELM(I)  /  ELMAX(M) 

100  CONTINUE 


Determination  of  UiU(M) 

C 

SUM  -  0.0 

DO  110  I  -  1,  NCOMP  ' 

SUM  -  SUM  +  (ER(I)*TT(I)) 

110  CONTINUE 
C 

LMU(M)  -  CHISQD(1-ALPHA,2*(NFC(M) -NCOMP) )/(2*SUM) 

C 

C  Compute  estimate  of  overall  reliability  RSL(M)  for  the  system. 
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RSL(M)  -  1.0 
DO  130  J  -  1,  NCS 
PROD  =1.0 

DO  120  I  =  1,  COMP(J.l) 

K  -  C0MP(J.I+1) 

IF  (TY(K).EQ.l)  THEN 

PROD  -  PR0D*(1  -  SURV(TY(K) ,LMU(M)*ER(K) ,EBETA(K) ,UT(K))) 
ELSEIF  (TY(K).EQ.2)  THEN 

PROD  -  PR0D*(1  -  SURV(TY(K) , (LMU(M)*ER(K) )**(!. /EBETA(K)) , 
*  EBETA(K) ,UT(K))) 

ELSEIF  (TY(K).EQ.3)  THEN 

PROD  -  PROD*(l  -  SURV(TY(K) ,1.D0-LMU(M)*ER(K) ,0.D0,UT(K))) 
ENDIF 

120  CONTINUE 

REL2(J)  -  1.0  -  PROD 
RSL(M)  -  RSL(M)  *  REL2(J) 

130  CONTINUE 


Increment  replication  counter. 
M  =  M  +  1 


ENDDO 

END 


4.  Subroutine  (EVAL) . 


SUBROUTINE  EVAL 

This  subroutine  calls  BUBBLE  to  sort  the  array  RSL(NREP)  in 
ascending  order  to  get  an  ordered  array  ORSL(NREP) .  It  also 
determine  the  estimate  for  RSLOW  at  the  specified  significance 
level  ALPHA  and  the  value  of  LEVEL  in  which  ORSL(LEVEL)  is  closest 
to  the  true  reliability  RS . 

Include  the  declaration  files 
and  declare  the  local  variables. 

INCLUDE  'PARMl  DEF' 

INTEGER  INDEX 
REAL*8  DIFF 

Order  the  array  RSL(NREP)  in  ascending  order. 

DO  10  M  =  1,  NREP 
ORSL(M)  -  RSL(M) 

10  CONTINUE 
C 

C  Bubble  Sort.  Sink  the  larger  of  the  pair. 
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C 

CALL  BUBBLE (NREP.RSL.ORSL) 

C 

C  Determine  the  (1-ALPHA)%  lower  confidence  bound  for  the  system 
C  reliability. 

C 

RSLOW  -  0RSL(NINT(NREP*(1 -ALPHA))) 

C 

C  Finding  the  %  confidence  level  for  the  true  reliability  RS. 

C  (ie.  the  proportion  of  RSL(M)  lesser  than  RS) 

C 

DIFF  -  1.0 
INDEX  -  0 

DO  200  M  -  1,  NREP 

IF  (ABS(ORSL(M)-RS) .LT.DIFF)  THEN 
DIFF  -  ABS(ORSL(M)-RS) 

INDEX  -  M 
ENDIF 

200  CONTINUE 
C 

LEVEL  -  FLOAT (INDEX) /NREP 
Record  evaluated  parameters  in  RAW1.DAT  (unit  2). 


OPEN (UNIT-2) 

WRITE(2,300) 

300  FORMAT (IX,  '  M  LMU(M)  ELMAX(M)  RSL(M) ’ , 

*  '  ORSL(M)  NFC(M)') 

DO  500  M  -  1,  NREP 

WRITE(2,400)  M,LMU(M) ,ELMAX(M) ,RSL(M) ,0RSL(M) ,NFC(M) 
400  F0RMAT(1X.I6,2F12.7,2F12.7,I10) 

500  CONTINUE 

CLOSE (UNIT-2) 

RETURN 

END 


5.  Report  Generation  Subroutine  (REPORT). 


SUBROUTINE  REPORT 
C 

C  This  subroutine  record  the  simulation  results  into  the  ’ 0UT1.DAT' 
C  file  as  logic  unit  3. 

C 

C  Include  the  declaration  files 
C  and  declare  local  variables. 

C 

INCLUDE  'PARMl  DEF' 

INTEGER  I,  J,  K,  DUM(IO) 

C 
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Write  to  output  file  '0UT1.DAT'  designated  as  logic  unit  3. 

0PEN(UNIT=3) 

C 

WRITE(3,10) 

WRITE(3,20)  NREP 
WRITE(3,25) 

WRITE(3.26) 

WR1TE(3.30) 

WRITE(3,40) 

WRITE ( 3 , 50 )  ISEED , NCOMP , ALPHA , TOL . NCS , TCN 
C 

WRITE(3,60) 

DO  200  I  -  1,  NCOMP 

WRITE(3.70)  I,TY(I) ,PARM(1,I),PARM(2,I) ,UT(I) ,NC(I) ,NF(I) 

200  CONTINUE 
WRITE(3.80) 

WRITE(3,90) 

DO  300  I  -  1,  NCOMP 

WRITE(3,100)  I,NF(I) ,TT(I) .ELM(I) ,ER(I) ,EBETA(I) 

300  CONTINUE 

WRITE(3,110) 

WRITE(3,120) 

DO  500  J  -  1,  NCS 
DO  400  K  -  1,  10 
DUM(K)  -  COMP(J,K) 

400  CONTINUE 

WRITE(3,130)  J,DUM,REL1(J) .REL2(J) 

500  CONTINUE 

WRITE(3,1A0) 

WRITE(3,150)  RS,ELMAX(NREP) ,LMU(NREP) ,RSL0W, LEVEL 
C 

10  FORMAT (IX, 'OUT 1 .DAT  :  Output  File  of  the  RETPl  simulation') 

20  FORMAT ( IX, '  after  ',15,'  replications',/) 

25  FORMAT (IX, 'COMMENTS  :  8  COMPONENTS  IN  SERIES  ') 

26  FORMAT(lX,'  DF  -  2  *  (NFC  -  NCOMP)  ' ,/) 

30  FORMAT ( IX, ' Input  Parameters:',/) 

40  FORMATdX,'  ISEED  NCOMP  ALPHA  TOL  NCS  TCN',/) 

50  FORMAT(1X,F10.1,I8,F8.4,F8.5,2I6,/) 

60  FORMATdX,'  ITY(I)  PARMl(I)  PARM2(I)  UT(I)  NC(I)  NF(I)',/) 
70  FORMAT(lX,I3,I6,2F9.5,F8.2,2I8) 

80  FORMAT (IX,/, 'Output  Parameters  for  the  LAST  Replication:',/) 

90  FORMATdX,'  INF(I)  TT(I)  ELM(I)  ER(I)', 

*  '  EBETA(I)',/) 

100  FORMAT(lX,I3,I6,E16.7,2F14.7,F14.7) 

110  FORMATdX,/, 'Cut-Set  Data:  '  ,/) 

120  F0RMAT(1X,'  J  NUM  Component  List 

*  '  RELl  REL2(M)',/) 

130  FORMATdX, 13, 15, 913, 2F12. 9) 

140  FORMATdX,/,'  RS  ELMAX(M)  LMU(M)  '  , 
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* 


RSLOW 


LEVEL' ,/) 


150  F0RMAT(1X,5F12.7,/) 

C 

CLOSE (UNIT-3) 

C 

RETURN 

END 

C . 

C  This  portion  of  the  file  contains  functions  and  subroutines 
C  used  in  the  RETPl  model . 

C  -  14  May  91 

C  -  by  Yee  Kah-Chee  SMC  2802 

C . 

C  A.  Random  Number  Generating  Subroutine  (LRNDPC) . 

C  (Courtesy  of  Mr.  David  Lim  Hung-Heng) 

C . 

SUBROUTINE  LRNDPC  (DSEED.U.N) 

INTEGER  N,  I 

REAL*8  U(N) 

REAL*8  D31M1,  DSEED,  D31 

D31M1-2**31  -  1 
D31  -2**31 

DATA  D31M1/2147483647.D0/ 

DATA  D31  /2147483648.D0/ 

DO  5  I-1,N 

DSEED  -  DMOD(950706376.D0*DSEED.D31Ml) 

DSEED  -  DMOD(16807.D0*DSEED,D31Ml) 

5  U(I)  -  DSEED  /  D31 
RETURN 
END 


B.  Survivability  Function. 


FUNCTION  SURV( TYPE, PARI, PAR2, UTIL) 

C 

C  This  function  returns  the  survival  probability  of  the  component  of 
C  different  types  (TYPE)  with  scale  (PARI)  and  shape  (PAR2)  parameters 
C  given  the  specified  utilization  times  or  cycles  (UTIL). 

C 

INTEGER  TYPE,  N 
REAL*8  PARI,  PAR2 ,  UTIL 
C 

IF  (TYPE.EQ.l)  THEN 

SURV  -  EXP(-(PAR1*UTIL)) 

ELSEIF  (TYPE.EQ.2)  THEN 

SURV  -  EXP(-((PAR1*UTIL)**PAR2)) 

ELSE 

N  -  NINT(UTIL) 

SURV  -  PAR1**N 
END  IF 
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C 


END 


C.  Bubble  Sort  Routine  in  ASCENDING  Order. 


SUBROUTINE  BUBBLE(N,LIST,OLIST) 

This  subroutine  performs  a  bubble  sort  in  increasing  order  (ie.  sink 
the  greater  numeral)  for  the  first  N  terms  in  an  array  LIST  and 
returns  the  result  in  OLIST. 

LOGICAL  DONE 
INTEGER  N,  K,  PAIR 
REAL*8  LIST(*),  OLIST(*) 

Sink  the  larger  of  the  pair. 

DO  50  K  -  1.  N 

OLIST(K)  -  LIST(K) 

50  CONTINUE 

PAIR  -  N  -  1 
DONE  -  .FALSE. 

DO  WHILE  ( .NOT. DONE) 

DONE  -  .TRUE. 

DO  100  K  -  1,  PAIR 

IF  (OLIST(K) .GT.0LIST(K+1))  THEN 
TEMP  -  OLIST(K) 

OLIST(K)  -  0LIST(K+1) 

0LIST(K+1)  -  TEMP 
DONE  -  .FALSE. 

ENDIF 

100  CONTINUE 

PAIR  -  PAIR  -  1 
ENDDO 
END 

. . 

C  D.  Unbiasing  Factor  for  Biased  MLE  for  Weibull  Shape  Parameter. 

C . - . - . 

FUNCTION  BN(I) 

C 

C  This  functon  returns  the  value  of  the  unbiased  factor  for  the  biased 
C  maximum  likelihood  estimate  of  the  shape  parameter  of  a  Weibull 
C  distribution  with  a  sample  size  of  N. 

C 

INTEGER  I 
IF  (I.LE.5)  THEN 

BN  -  (I*0.699)/(5.0) 

ELSEIF  (I.EQ.6)  THEN 
BN  -  0.752 

ELSEIF  (I.EQ.7)  THEN 
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BN  -  0.786 

ELSEIF  (I.EQ.8)  THEN 
BN  -  0.82 

ELSEIF  (I.EQ.9)  THEN 
BN  -  0.8395 
ELSEIF  (I.EQ.IO)  THEN 
BN  -  0.859 

ELSEIF  (I.EQ.ll)  THEN 
BN  -  0.871 

ELSEIF  (I.EQ.12)  THEN 
BN  -  0.883 

ELSEIF  (I.EQ.13)  THEN 
BN  -  0.892 

ELSEIF  (I.EQ.14)  THEN 
BN  -  0.901 

ELSEIF  (I.EQ.15)  THEN 
BN  -  0.9075 
ELSEIF  (I.EQ.16)  THEN 
BN  -  0.914 

ELSEIF  (I.EQ.17)  THEN 
BN  -  0.9185 
ELSEIF  (I.EQ.18)  THEN 
BN  -  0.923 

ELSEIF  (I.EQ.19)  THEN 
BN  -  0.927 

ELSEIF  (I.EQ.20)  THEN 
BN  -  0.931 

ELSEIF  (I.LE.25)  THEN 

BN  -  0.931+(I-20)*0. 014/5.0 
ELSEIF  (I.LE.30)  THEN 

BN  -  0.945+(I-25)*0. 01/5.0 
ELSEIF  (I.LE.40)  THEN 

BN  -  0.955+(I-30)*0. 011/10.0 
ELSEIF  (I.LE.60)  THEN 

BN  =  0.966+(I-40)*0. 012/20.0 
ELSEIF  (I.LE.80)  THEN 

BN  -  0.978+(I-60)*0. 006/20.0 
ELSEIF  (I.LE.lOO)  THEN 

BN  -  0.984+(I-80)*0. 003/20.0 
ELSEIF  (I.LE.120)  THEN 

BN  -  0.987+(I-100)*0. 003/20.0 
ELSE 

BN  -  1.0 
ENDIF 
RETURN 
END 

C . 

C  E.  Biased  MLE  of  Weibull  Shape  Parameter. 

C . . 

SUBROUTINE  MLESHAPE (T , N , R , DEL, B , BNEW) 
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c 

c 

c 

c 

c 

c 


c 


c 

c 


c 


c 


c 


c 

c 

c 

c 

c 

c 

c 


This  subroutine  returns  a  biased  estimator  (BNEW)  for  a  Weibull 
shape  parameter  using  the  Newton-Raphson' s  Method  of  Successive 
Approximation.  The  data  parameters  consist  of  an  ascending  ordered 
list  of  failure  times  (T)  ,  S€imple  size  (N) ,  number  of  failed  samples 
(R) ,  tolerance  for  convergence  (DEL)  and  an  initial  estimate  of  the 
shape  parameter  (B) . 

LOGICAL  DONE 

INTEGER  N,  R,  I 

REAL*8  GFUNCT,  GPRIME,  B,  BOLD,  BNEW,  T(*) ,  DEL, 

*  TERMl,  TERM2,  TERMS,  SUMl,  SUM2,  SUMS,  SUM4,  STEP 

BNEW  -  B 

DONE  -  .FALSE. 


DO  WHILE  (.NOT. DONE) 

DONE  -  .TRUE. 

TERMl  -  FLOAT(N-R)*(T(R)**BNEW) 

TERM2  -  FLOAT(N-R)*(T(R)**BNEW)*LOG(T(R)) 

TERMS  -  FLOAT(N-R)*(T(R)**BNEW)*LOG(T(R))*LOG(T(R)) 

SUMl  -  0.0 

SUM2  -  0.0 

SUMS  -  0.0 

SUM4  -  0.0 


DO  50  I  -  1,  R 

SUMl  -  SUMl  +  T(I)**BNEW 

SUM2  -  SUM2  +  (T(1)**BNEW)*L0G(T(I)) 

SUMS  -  SUMS  +  (T(I)**BNEW)*L0G(T(I))*L0G(T(I)) 

SUM4  -  SUM4  +  L0G(T(I)) 

50  CONTINUE 


* 


GFUNCT  =  (SUM2+TERM2)/(SUM1+TERM1) 
-  (1.0/FLOAT(R))*SUM4 


(1.0/BNEW) 


GPRIME  -  (1.0/(SUM1+TERM1)**2)*(  (SUM1+TERM1)*(SUMS+TERMS) 

*  -  (SUM2+TERM2)**2  ) 

*  +  (1.0/BNEW**2) 


PRINT  60,  GFUNCT, GPRIME, BNEW 

60  F0RMAT(1X,  'GFUNCT  -',F8.S,'  GPRIME  -',F8.3,'  BNEW  -',F8.S) 


Control  magnitude  of  the  marching  step  towards  convergence 
as  no  more  than  0.1. 


IF  ( (GFUNCT. LT.O)  .AND.  (GPRIME. GT . 0) )  THEN 
STEP  ■=■  VMAX(- .IDO,  (GFUNCT/GPRIME)) 

ELSEIF  ((GFUNCT. GT.O)  .AND.  (GPRIME . LT . 0) )  THEN 
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STEP  -  VMAX(-.1D0.(GFUNCT/GPRIME)) 

ELSE 

STEP  -  VMIN( .IDO, (GFUNCT/GPRIME)) 

ENDIF 

BOLD  -  BNEW 
BNEW  -  BNEW  -  STEP 

Check  for  convergence  of  the  MLE  for  the  shape  parameter  B. 

IF  ( ABS( BOLD- BNEW ) .GT. DEL)  THEN 
DONE  -  .FALSE. 

ENDIF 

Avoid  overflow  error  due  to  large  MLE  value  caused  by  small 
GPRIME  (slope)  as  GFUNCT  approaches  to  near  zero. 

Stop  when  magnitude  of  BNEW  exceeds  7. 

IF  (BNEW. GT. 7.0)  THEN 
BNEW  -  BOLD 
DONE  -  .TRUE. 

ENDIF 

ENDDO 

RETURN 

END 


E.  Chi-Square  Quantile  Function. 


FUNCTION  CHISQD(P,N) 

C 

C  Modified  version  of  Algorithm  451  from  Comunications  of  the  ACM 
C  Aug  1977  Vol.16  No. 8  . 

C 

C  This  function  evaluates  the  quantile  at  the  probability  level  P 
C  (left  tail  area)  for  the  Chi-square  distribution  with 
C  N  degrees  of  freedom. 

C 

REAL* 8  P 
REAL  X 
INTEGER  IF 

DIMENSION  C(21) ,  A(19) 


DATA 

C/  1.565326E-3, 

* 

1.060438E-3, 

* 

-6.950356E-3, 

* 

-1.323293E-2. 

* 

2.277679E-2, 

* 

-8.986007E-3, 

* 

-1.513904E-2, 

* 

2.530010E-3, 
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* 

1.450117E-3. 

* 

5.169654E-3. 

* 

.1.153761E-2, 

* 

1.128186E-2. 

* 

2.607083E-2. 

* 

■0.2237368, 

* 

9.780499E-5, 

* 

-8.426812E-4. 

* 

3.125580E-3, 

★ 

-8.553069E-3, 

* 

1.348028E-4, 

* 

0.4713941, 

* 

1.0000886  / 

DATA  A/ 

1.264616E-2, 

* 

-1.425296E-2. 

* 

1.400483E-2, 

* 

-5.886090E-3, 

* 

-1.091214E-2, 

* 

-2.304527E-2, 

* 

3.135411E-3, 

* 

-2.728484E-4, 

* 

-9.699681E-3, 

* 

1.316872E-2, 

* 

2.618914E-2, 

* 

-0.2222222. 

* 

5.406674E-5, 

* 

3.483789E-5, 

* 

.7.274761E-4, 

* 

3.292181E-3, 

* 

.8.729713E-3, 

* 

0.4714045, 

* 

1.  / 

IF  (N-2)  10,  20,  30 
10  CALL  XFR0MP( .5*(1. -P) ,X,IF) 


20 


30 


-  -2.*L0G(1. -P) 


CHISQD  -  X 

CHISQD  -  CHISQD*CHISQD 
RETURN 
CHISQD 
RETURN 
F  -  N 
FI  -  l./F 

CALL  XFROMP(P,X,IF) 

T  -  X 

F2  -  SQRT(F1)*T 

IF  (N.GE. (2+INT(4,*ABS(T))))  GO  TO  40 

CHTSOD  -  (((((((C(1)*F2+C(2))*F2+C(3))*F2+C(4))*F2 

k  !c(5))*F2+C(6))*F2+C(7))*Fl+((((((C(8)+C(9)*F2)*F2 

t  +C(10))*F2+C(11))*F2+C(12))*F2+C(13))*F2+C(14)))*F1+ 

*  (((((C(15)*F2+C(16))*F2+C(17))*F2+C(18))*F2 

^  +C(19))*F2+C(20))*F2+C(21) 
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GO  TO  50 

40  CHISQD-(((A(1)+A(2)*F2)*F1+(((A(3)+A(4)*F2)*F2 

*  +A(5))*F2+A(6)))*F1+(((((A(7)+A(8)*F2)*F2+A(9))*F2 

*  +A(10) )*F2+A(11) )*F2+A(12) ) )*Fl+( ( ( ( (A(13)*F2 

*  +A(14))*F2+A(15))*F2+A(16))*F2+A(17))*F2*F2 

*  +A(18))*F2+A(19) 

50  CHISQD  -  CHISQD*CHISQD*CHISQD*F 
RETURN 
END 


F.  Standard  Normal  Variate  Computation  Subroutine. 


C 

C 

C 


SUBROUTINE  XFROMP(P,X,IFAULT) 

Algorithm  AS  24  J .R.STAT.SOC.  C.  (1969)  Vol.18.  No. 3. 

This  subroutine  computes  the  standard  normal  deviate  X  for 
the  specified  left  tail  area  P. 

REAL*8  P 
DIMENSION  A(5) 

DIMENSION  CONNOR  (17),  HSTNGS(6) 

DATA  CONNOR/  8 . 0327350124E-17 , 


* 

1.4483264644E-15, 

* 

2.4668270103E-14, 

* 

3.9554295164E-13, 

★ 

5.9477940136E-12, 

* 

8.3507027951E-11, 

* 

1.0892221037E-9, 

* 

1.3122532964E-8, 

* 

1.4503852223E-7, 

* 

1.4589169001E-6, 

* 

1.3227513228E-5, 

* 

1.0683760684E-4, 

* 

7.5757575758E-4, 

* 

4.6296296296E-3, 

* 

2.3809523810E-2, 

* 

0.1. 

* 

0.33333333333  / 

DATA  RTHFPI  / 
DATA  RRT2PI  / 
DATA  TERMIN  / 
DATA  HSTNGS  / 


* 

* 

* 


1.2533141373  / 

0.3989422804  / 

l.OE-11  / 

2.515517, 

0.802853, 

0.010328, 

1.432788, 
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C 


* 

* 


0.189269. 
0.001308  / 


I FAULT  -  1 

IF  ((P.LE.0.0) .OR. (P.GE.1.0))  GO  TO  100 
I FAULT  -  0 
C 

Get  first  approximation  XO  to  deviate  by  Hastings'  formula 
B  -  P 

IF(B.GT.0.5)  B  -  1.0  -  B 

F  -  -  LOG(B) 

E  -  SQRT(F+F) 

XO  -  -E  +  ((HSTNGS(3)*E+HSTNGS(2))*E+HSTNGS(1))/ 

*  ( ( ( HSTNGS ( 6 ) *E+HSTNGS ( 5 ) ) *E+HSTNGS ( 4 ) ) *E+1 . 0 ) 

IF  (XO.LT.O.O)  GO  TO  1 
XO  -  0.0 
PO  -  0.5 
XI  -  -RTHFPI 
GO  TO  7 

Find  the  area  PO  corresponding  to  XO 

1  Y  -  X0**2 

IF  (X0.LE.-1.9)  GO  TO  3 
Y  -  -0.5*Y 

(1)  series  approximation 

PO  -  CONNOR(l) 

DO  2  L-2,17 

2  PO  -  PO*Y  +  CONNOR (L) 

PO  _  (PO*Y+1.0)*XO 
XI  -  -(PO+RTHFPI)*EXP(-Y) 

PO  =  PO*RRT2PI  +0.5 
GO  TO  7 

(2)  continued  fraction  approximation 
C 

3  Z  -  1.0 A 
A(2)  -  1.0 
A(3)  -  1.0 
A(4)  -  Z  +  1.0 
A(5)  -  1.0 

W  -  2.0 
C 

4  DO  6  L-1,3.2 
DO  5  J-1,2 

K  =  L  +  J 
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KA  =  7  -  K 


C 

5  A(K)  -  A(KA)  +  A(K)*W*Z 
C 

6  W  -  W  +  1.0 
APPRXU  -  A(2)/A(3) 

APPRXL  -  A(5)/A(4) 

C  -  APPRXU  -  APPRXL 
IF  (C.GE.TERMIN)  GO  TO  4 
XI  -  APPRXL/XO 
PO  -  -X1*RRT2PI*EXP(-0.5*Y) 

C 

C  Get  accurate  value  of  deviate  by  Taylor  Series 
C  (XI,  X2 ,  X3  are  derivatives  for  the  Taylor  Series 
C 

7  D  -  F  +  LOG(PO) 

X2  -  X0*X1*X1  -XI 

X3  -  Xl**3  +  2,0*XO*X1*X2  -X2 
X  -  ((X3*D/3.0+X2)*D/2.0+Xl)*D  +  XO 
IF  (P.LE.0.5)  GO  TO  100 
X  -  -X 
100  RETURN 
END 


G.  Maximum  Function. 


FUNCTION  VMAX(X,Y) 

REAL*8  X,  Y 
IF(X.GT.Y)  THEN 
VMAX  -  X 
ELSE 

VMAX  -  Y 
END  IF 
RETURN 
END 

C . 

C  H.  Minimum  Function. 

C . 

FUNCTION  VMIN(X,Y) 

REAL*8  X,Y 
IF(X.LT.Y)  THEN 
VMIN  -  X 
ELSE 

VMIN  -  Y 
END  IF 
RETURN 
END 

4.  Program  Output.  (OUTl.DAT) 
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The  result  for  the  simulation  run  based  on  the  input  parameters  specified  in 
INl.DAT  are  computed  and  written  to  the  file  OUTl.DAT.  A  sample  of  this  file  is 
as  follows. 

OUTl.DAT  :  Output  File  of  the  RETPl  simulation 
after  1000  replications 

COMMENTS  :  8  COMPONENTS  IN  SERIES 

DF  -  2  *  (NFC  -  NCOMP) 


Input  Parameters: 


I  SEED 

NCOMP 

ALPHA 

TOL  NCS 

TCN 

16807.0 

8 

0.2000  0. 

01000  8 

3 

I 

TY(I)  PARMl(I) 

PARM2(I) 

UT(I)  NC(I) 

NF(I) 

1 

1  0.00200 

1.00000 

10.00 

15 

3 

2 

1  0.00400 

1.00000 

10.00 

15 

3 

3 

1  0.00600 

1.00000 

10.00 

15 

3 

4 

1  0.00800 

1.00000 

10.00 

15 

3 

5 

2  0.00200 

2.00000 

10.00 

15 

3 

6 

2  0.00400 

2.00000 

10.00 

15 

3 

7 

2  0.00600 

2.00000 

10.00 

15 

3 

8 

2  0.00800 

2.00000 

10.00 

15 

3 

Output  Parameters  for  the  LAST 

Replication: 

I 

NF(I) 

TT(I) 

ELM(I) 

ER(I) 

EBETA(I) 

1 

3 

0.1969091E+04 

0.0010157 

0.1703872 

0.0000000 

2 

3 

0. 3355078E+03 

0.0059611 

1.0000000 

0.0000000 

3 

3 

0.5999981E+03 

0.0033333 

0.5591815 

0.0000000 

4 

3 

0.6423594E+03 

0.0031135 

0.5223055 

0.0000000 

5 

3 

0.2563836E+16 

0.0000000 

0.0000000 

6.3524983 

6 

3 

0.8186629E+07 

0.0000004 

0.0000615 

2.9571536 

7 

3 

0. 5675037E+09 

0.0000000 

0.0000009 

4.2418658 

8 

3 

0. 5675519E+05 

0.0000529 

0.0088672 

2.3471863 

Cut-Set  Data: 

J  NUM  Component  List  RELl  REL2(M) 

1  1100000000  0.980198622  0.990279973 

2  1200000000  0.960789382  0.944286704 

3  1300000000  0.941764534  0.968452990 

4  1400000000  0.923116326  0.970502377 

5  1500000000  0.999600053  0.999999940 

6  1600000000  0.998401225  0.999680758 
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7  1700000000  0.996406436  0.999911249 

8  1800000000  0.993620396  0.988757312 


RS 

ELMAX(M) 

LMU(M) 

RSLOW 

LEVEL 

0.8089645 

0.0059611 

0.0057325 

0.8691733 

0.3380000 
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APPENDIX  C  :  Users’  Guide  for  RETP2 


Reliability  Estimation  Test  Plan  2  (RETP2). 
by  YEE,  Kah-Chee 
July  91 


1.  Brief  Description. 

RETP2  is  a  computer  program  written  in  FORTRAN  that  runs  on  the  Amdahl 
mainframe  at  NPGS.  It  allows  the  user  to  simulate  exponential  and  Weibull  failure 
times  of  component  items  being  tested  to  evaluate  the  accuracy  of  a  confidence  limit 
estimation  procedure  based  on  Type  I  data  censoring  (that  is,  testing  items  of 
component  i  until  a  specified  total  testing  time  is  achieved  for  each  of  them), 

2.  Program  Input.  (IN2.DATA) 

The  input  of  the  program  are  specified  to  the  program  via  an  input  file  called 
IN2.DAT.  A  sample  input  file  is  shown  below. 

This  file  contains  the  inputs  required  by  the  RETP2  model. 

Update  only  the  numerical  values  between  dotted  lines  as  appropriate. 

Do  not  delete  any  of  the  comment  lines.  (IN2.DAT)  20  Jun  91 


Value 

Type 

Units 

Description 

Variable 

16807.0 

Real 

- 

initial  random  seed 

ISEED 

8 

Int 

- 

total  of  components  in  system 

NCOMP 

4 

INT 

- 

//  OF  EXPONENTIAL  COMPONENTS 

NEXP 

4 

INT 

- 

#  OF  WEIBULL  COMPONENTS 

NWEI 

0 

INT 

- 

#  OF  GEOMETRIC  COMPONENTS 

NGEO 

0.01 

Real 

- 

tolerance  for  MLE 

TOL 

0.20 

REAL 

- 

DESIRED  SIGNIFICANCE  LEVEL 

ALPHA 

1000 

INT 

- 

#  OF  REPLICATIONS  DESIRED 

NREP 

3 

INT 

TEST  CASE  NUMBER 

1  -  all  EXP 

2  -  all  WEI 

3  -  EXP  +  WEI 

4  -  EXP  +  WEI  +  CYC 

TCN 

8 

Int 

- 

number  of  cut  sets 

NCS 
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C  TEST  PLAN  :  Testing  until  TT(I)  (total  test  time)  is  accumulated. 
C  (Use  REAL  numbers  ONLY! I ! ) 


Comp 

Number 

I 

Int 

Comp 

Type 

TY(I) 

Int 

Comp  Parameters 
Scale  Shape 

PARMl(I)  PARM2(I) 

Real  Real 

Util 

Time/Cycle 

UT(I) 

(hrs) (eyes) 

Test  Plan  Inputs 
Accumulated 

TT(I)  NC(I) 

(hrs) (eyes)  Int 

1.0 

1.0 

0.005 

1.0 

5.0 

5400.0 

20.0 

2.0 

1.0 

0.005 

1.0 

5.0 

5400 . 0 

20.0 

3.0 

1.0 

0.005 

1.0 

5.0 

5400.0 

20.0 

4.0 

1.0 

0.005 

1.0 

5.0 

5400.0 

20.0 

5.0 

2.0 

0.010 

2.0 

15.0 

2700.0 

20.0 

6.0 

2.0 

0.010 

2.0 

15.0 

2700.0 

20.0 

7.0 

2.0 

0.010 

2.0 

15.0 

2700.0 

20.0 

8.0 

2.0 

0.010 

2.0 

15.0 

2700.0 

20.0 

Note  ; 

TY(I) 

-1  EXPONENTIAL  P(surv) 

-  exp( 

-PARM2)*T) 

TY(I) 

-2  WEIBULL  P(surv) 

-  exp( 

- ( PARM1*T ) **PARM2 ) 

TY(I) 

-3  GEOMETRIC  P(surv) 

-  PARM1**T 

SYSTEM 

CONFIGURATION  : 

Identification 

of  CUT 

SETs 

-  min  groups 

of 

components  that  have  to  fail 

for  the  system  to  fail 

Cut  Set 

#  in  Set 

List  of  Components  in 

i  Cutset 

J 

C0MP(J,1) 

C0MP(J,2)  . 

up  to  C0MP(J,1)  components 

1 

1 

10  0  0 

0 

0  0 

0 

0  0 

2 

1 

2  0  0  0 

0 

0  0 

0 

0  0 

3 

1 

3  0  0  0 

0 

0  0 

0 

0  0 

4 

1 

4  0  0  0 

0 

0  0 

0 

0  0 

5 

1 

5  0  0  0 

0 

0  0 

0 

0  0 

6 

1 

6  0  0  0 

0 

0  0 

0 

0  0 

7 

1 

7  0  0  0 

0 

0  0 

0 

0  0 

8 

1 

8  0  0  0 

0 

0  0 

0 

0  0 

3.  Program  Flow  and  Logic.  (NAME2.DEF,  PARM2.DEF  and  RETP2.FOR) 

Input  parameters  are  first  read  in  by  the  program  by  calling  the  INPUT 
subroutine.  The  program  then  evoke  the  SIM  subroutine  which  generates  the 
random  failure  times  and  compute  the  key  statistics  required  in  the  procedure.  The 
next  subroutine  EVAL  determines  the  measures  of  accuracy  for  run.  REPORT  is 
the  subroutine  which  generates  the  output  file  for  the  run  OUT2.DAT. 


The  variables  in  the  program  RETP2.FOR  are  described  in  the  file 
NAME2.DEF  as  listed  below. 
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C . . . 

C  This  file  contains  the  declaration  for  input  and  output  variables 
C  used  in  the  the  RETP2  model.  (NAME2.DEF)  20  Jun  91 

C . - . - . 

C  Input  Variables. 

C  . 

C  ISEED  “  initial  random  seed  selected. 

C  SEED  -  current  random  seed. 

C  RS  “  true  overall  series  system  reliability. 

C  ALPHA  -  level  of  significance  desired. 

C  NREP  —  number  of  replications  desired  for  the  simulation. 

C  TPN  -  test  plan  number  (1). 

C  TCN  -  test  case  number  (1,  2,  3  or  4) . 

C  NCOMP  -  total  number  of  components  in  the  system. 

C  NEXP  -  number  of  components  with  EXP  failure  times. 

C  NWEl  -  number  of  components  with  WEI  failure  times. 

C  NGEO  -  number  of  components  with  GEO  failure  times. 

C  TOL  -  desired  tolerance  for  MLE  of  WEI  shape  parameter. 

C  Distribution:  Exponential  WEIbull  GEOmetric 

C  TY(I)  -  type:  1  23 

C  PARM(1,I):  Scale(l/hr)  Scale(l/hr)  Prob 

C  PARM(2,I):  -  Shape 

C  UT(I)  -  utilization  time  (hrs)  for  component  i  (EXP  and  WEI). 

C  UC(I)  -  utilization  cycles  for  component  i  (GEO  only). 

C  NC(1)  -  number  of  test  samples  (sample  size)  for  cmponent  i. 

C  NF(I)  -  desired  number  of  failures  in  test  for  component  i. 

C  NCS  -  number  of  cut -sets  for  the  system. 

C  C0MP(J,K)  -  kth  parameter  of  cut-set  j  (first  being  the  no.  of 
C  components  belonging  to  the  cut-set). 

C 

C  Assumed  Variables. 

C  . 

C  MAXCOMP  -  maximum  number  of  components  allowed  in  the  system. 

C  MAXREP  -  maximum  number  of  replications  permitted. 

C  MAXCUT  -  maximum  number  of  cut-sets. 

C 

C  Program  and  Output  Variables. 

C  . 

C  RS  -  true  overall  system  reliability. 

C  TT(I)  -  total  accumulated  failure  time  (hr)  for  component  i 
C  (EXP  and  WEI  only). 

C  TC(I)  -  total  accumulated  cycles  to  failure  (incl.  failure  cycle) 

C  for  component  i  (GEO  only) . 

C  EBETA(I)  -  estimate  for  shape  parameter  of  component  i  (if  Weibull) . 

C  RELl(J)  -  actual  reliability  for  cut-set  j. 

C  REL2(J)  -  computed  reliability  for  cut-set  j  for  current  replication. 
C  ELM(I)  -  estimated  component  failure  rate  (1/hrs)  for  component  i. 

C  ELMAX(M)  -  max  estimated  component  failure  rate  for  rep  m  (1/hrs) . 

C  ER(I)  “  ratio  of  estimated  failure  rate  to  IJIAX. 

C  ET(I)  -  same  as  TT(I)  except  that  these  are  for  Weibull  components. 
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C  NFC(M)  —  total  number  of  failed  test  components. 

C  LMU(M)  -  upper  confidence  bound  for  failure  rate  (1/hrs), 

C  RSL(M)  -  lower  confidence  limit  estimated  for  system  reliability 

C  ...  for  the  mth  replication. 

C  ORSL(M)  =  ordered  RSL(M)  (ascending). 

C  RSLOW  -  (l-ALPHA)xlOO  percentile  of  set  of  RSL(M) . 

C  LEVEL  -  achieved  confidence  level,  ie.  proportion  of  RSL(M)  that 
C  are  lesser  than  RS  (conservative  estimate) . 

C 

C--  END  OF  NAME2.DEF . - . - . 

C 


Together  with  the  main  program  in  RETP2.FOR  are  the  other  subroutines 
needed  in  the  simulation.  The  declaration  of  variables  is  done  in  the  file 
PARM2.DEF.  Relevant  descriptions  are  included  as  comment  lines  in  the  source 
code  to  help  explain  the  program  segments.  A  listing  of  PARM2.DEF  and 
RETP2.FOR  is  given  below. 

c . . . - . 

C  This  file  contains  the  declaration  for  input  and  output  variables 
C  used  in  the  the  RETP2  model.  (PARM2.DEF)  20  Jun  91 

C . 

INTEGER  MAXCOMP,  MAXREP 

PARAMETER (  MAXCOMP  -  100  ,  MAXREP  -  1000  ,  MAXCUT  -  20  ) 

REAL*8  ISEED,  SEED 

INTEGER  NREP,  TCN,  NCOMP,  NEXP,  NWEI ,  NGEO,  NCS , 

*  NC (MAXCOMP),  NF(MAXCOMP),  TY ( MAXCOMP ) ,  NFC (MAXREP) , 

*  UC(MAXCOMP),  TC(MAXCOMP),  COMP (MAXCUT, MAXCOMP) 

REAL*8  RS,  ALPHA.  UT ( MAXCOMP ) ,  TT ( MAXCOMP ) , 

*  PARM (2. MAXCOMP ) ,  ELM ( MAXCOMP ) ,  ER ( MAXCOMP ) ,  ET (MAXCOMP ) , 

*  LMU(MAXREP),  RSL(MAXREP) ,  ORSL ( MAXREP ) , 

*  ELMAX ( MAXREP ) ,  RSLOW,  LEVEL,  TOL,  EBETA (MAXCOMP) , 

*  RELl(MAXCUT) ,  REL2(MAXCUT) 

C 

COMMON/BLOCKl/ISEED.  SEED.  NREP,  TCN,  NCOMP,  NC ,  NF,  NEXP,  NWEI , 

*  NGEO,  NCS,  TY,  NFC,  UC,  TC ,  COMP 
COMMON/BLOCK2/RS,  ALPHA,  UT,  TT,  PARM,  ELM,  ER,  ET,  LMU, 

*  RSL,  ORSL,  ELMAX,  RSLOW,  LEVEL,  TOL,  EBETA, 

*  RELl ,  REL2 
C 

C--  END  OF  PARM2.DEF . - . 

C 

C . - . 

C  This  file  contains  the  main  program  and  the  subroutines 
C  for  the  Reliability  Estimation  Test  Plan  2  (RETP2)  model. 

C  (RETP2.F0R)  -  runs  on  a  IBM  PC  Compatible. 

C 

C  IBM  Mainframe  Version. 

C 
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Test  Plan  2  ;  Testing  until  accumulated  time  or  cycles  is  achieved 
.  for  component  i. 

by  Yee  Kah-Chee  SMC  2802. 

20  Jun  91. 


1.  Main  Program  (RETP2) . 


PROGRAM  RETP2 

Include  the  declaration  files. 

INCLUDE  'NAME2  DEF' 

INCLUDE  'PARM2  DEF' 

Read  in  input  data. 

CALL  INPUT 

Activate  simulation. 

CALL  SIM 

Process  and  evaluate  output  data. 
CALL  EVAL 

Generate  simulation  report. 

CALL  REPORT 

STOP 

END 


2.  Input  Initialisation  Subroutine  (INPUT). 


SUBROUTINE  INPUT 

This  subroutine  reads  in  the  inputs  for  the  RETP2  model. 
C 

C  Include  the  declaration  file. 

C 

INCLUDE  'PARM2  DEF' 

C 

INTEGER  I,  J,  K,  DUM2(11) 

REAL*8  DUM1(7) 

C 

C  Read  data  from  'IN2.DAT'  designated  as  logic  unit  1. 

C 


0PEN(UNIT=1) 
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READ(1,10) 

10  FORMAT (IX. /////) 

READ(1,*)  ISEED 
READ(1,*)  NCOMP 
READ(1,*)  NEXP 
READ(1,*)  NWEI 
READ(1,*)  NGEO 
READ(1,*)  TOL 
READ(1.*)  ALPHA 
READd,*)  NREP 
READ(1,*)  TCN 
READ(1,20) 

20  FORMAT (IX,///) 

READ(1,*)  NCS 
C 

READ(1,30) 

30  FORMAT (IX, ////////) 

C 

DO  50  K  -  1,  NCOMP 
READ(1,*)  DUMl 
C 

I  -  NINT(DUM1(1)) 

TY(I)  -  NINT(DUM1(2)) 

C 

IF  (TY(I).EQ.l)  THEN 
PARM(1,I)  -  DUM1(3) 
PARM(2,I)  -  DUM1(4) 
UT(I)  -  DUM1(5) 

TT(I)  -  DUM1(6) 

NC(I)  -  NINT(DUM1(7)) 
EBETA(I)  -  0 

ELSEIF  (TY(I).EQ.2)  THEN 
PARM(1,I)  -  DUM1(3) 
PARM(2,I)  -  DUM1(4) 
UT(I)  -  DUM1(5) 

TT(I)  -  DUM1(6) 

NC(I)  -  NINT(DUM1(7)) 
ELSEIF  (TY(I).EQ.3)  THEN 
PARM(1,I)  -  DUMl (3) 
PARM(2,I)  -  DUM1(4) 
UC(I)  -  NINT(DUM1(5)) 
UT(I)  -  DUM1(5) 

TT(I)  -  DUM1(6) 

NC(I)  -  NINT(DUM1(7)) 
ENDIF 

50  CONTINUE 
C 

KCAUd  ,  DU  ) 

60  FORMAT (IX, ///////////) 
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DO  80  I  -  1,  NCS 
READd,*)  DUM2 
J  -  DUM2(1) 

C0MP(J,1)  -  DUM2(2) 

DO  70  K  -  1,  C0MP(J,1) 
C0MP(J,K+1)  -  DUM2(K+2) 
70  CONTINUE 
80  CONTINUE 

CLOSE(UNIT-l) 

RETURN 

END 


3.  Subroutine  (SIM) 


SUBROUTINE  SIM 

This  subroutine  simulates  NREP  possible  outcomes  of  the  test  plan 
desired  in  order  to  obtain  the  raw  estimates  of  LMU(M)  and  RSL(M) 
for  each  of  the  replication. 

Include  the  declaration  file 
and  declare  local  variables. 

INCLUDE  'PARM2  DEF' 

INTEGER  I,  J,  K,  M,  ISUM,  KEY,  ICOUNT 
INTEGER  NCYC(MAXCOMP) .  NSYS 
REAL* 8  UNI 
REAL*8  SUM,  PROD, 

*  FT(MAXCOMP,200) ,  OFT(MAXCOMP, 200) , 

*  DFT(200),  DOFT(200) 

LOGICAL  FLAG 

SEED  -  ISEED 

Compute  overall  true  system  reliability  RS . 

RS  -  1.0 
DO  30  J  -  1,  NCS 
PROD  -  1.0 

DO  20  I  -  1,  C0MP(J,1) 

K  -  COMP(J,I+l) 

PROD  -  PROD*(  1  -  SURV(TY(K) ,PARM(1,K) ,PARM(2,K) ,UT(K))  ) 
20  CONTINUE 

RELl(J)  -  1.0  -  PROD 
RS  -  RS  *  RELl(J) 

30  CONTINUE 
C 

C  Start  of  Simulation. 

C  (Initialize  replication  counter  M) . 

C 
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M  =  1 

DO  WHILE  (M.LE.NREP) 

C 

C  Test  Plan  ;  Sample  and  determine  unknovm  NF(I) 

C  .  (with  known  NC(I),  Weibull  case)  until  TT(I)  is  reached. 

C 

C  Generate  NC(I)  failure  times,  put  them  in  ascending  order 
C  with  the  smallest  failure  time  on  the  top  of  the  list. 

C 

DO  60  I  -  1,  NCOMP 
C 

C  Exponential  components.  Test  each  component  until  it  fails,  check 
C  to  see  if  TT(1)  is  exceeded,  if  not,  carry  on  testing. 

C 

IF(TY(I) .EQ.l)  THEN 
SUM  -  0.0 
L  -  1 

DO  WHILE  (SUM.LE.TT(I)) 

CALL  LRNDPC (SEED, UNI, 1) 

FT(I,L)  -  -L0G(UNI)/PARM(1,I) 

SUM  -  SUM  +  FT(I,L) 

L  -  L  +  1 
ENDDO 

NF(I)  -  L  -  2 
C 

C  Weibull  components.  Generate  NC(I)  failure  times,  put  them  in 
C  an  ascending  order  with  the  smallest  failure  time  on  top  and 
C  determine  NF(I). 

C 

ELSEIF(TY(I) .EQ.2)  THEN 
DO  40  K  -  1,  NC(I) 

CALL  LRNDPC (SEED, UNI, 1) 

FT(I,K)  -  (1.0/PARM(1,I))*(-LOG(UNI))**(1.0/PARM(2,I)) 
DFT(K)  -  FT(I,K) 

40  CONTINUE 

CALL  BUBBLE (NC(I) ,DFT,DOFT) 

DO  42  K  -  1,  NC(I) 

OFT(I,K)  -  DOFT(K) 

42  CONTINUE 

SUM  -  0.0 
NF(I)  =  0 
DO  45  K  -  1,  NC(I) 

IF  (0FT(I,K) .LT.TT(I))  THEN 
NF(I)  -  NF(I)  +  1 
END  IF 

45  CONTINUE 

C 

C  GEOMETRIC  COMPONENTS.  DETERMINE  NF(I). 

C 

ELSEIF(TY(I) .EQ.3)  THEN 
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NF(I)  -  0 

DO  50  K  -  1,  TT(I) 

CALL  LRNDPCC SEED. UNI, 1) 

IF  (UNI.GT.PARM(I,1))  THEN 
NF(I)  -  NF(I)  +  I 
ENDIF 

50  CONTINUE 

C 

ENDIF 

60  CONTINUE 
C 

C  Determine  total  number  of  failed  test  components  as  well  as 
C  checking  for  zero  component  failure  or  just  failures  from  a 
C  SINGLE  component  (FLAG  will  be  set  to  .TRUE,  if  so). 

C 

ISUM  -  0 
I COUNT  -  0 
DO  70  I  -  I,  NCOMP 
IF  (NF(I).GT.O)  THEN 
I COUNT  -  I COUNT  +  I 
ENDIF 

ISUM  -  ISUM  +  NF(I) 

70  CONTINUE 
C 

NFC(M)  -  ISUM 
C 

IF  (ICOUNT.LE. I)  THEN 
FU\G  -  .TRUE. 

ELSE 

FLAG  -  .FALSE. 

ENDIF 

C 

Case  A  ;  More  than  ONE  component  type  experienced  failures 
.  in  the  test. 

IF  (.NOT. FLAG)  THEN 

Estimate  the  failure  rate  of  each  component. 

DO  90  I  -  I,  NCOMP 

Exponential  components. 

C 

IF  (TY(I) .EQ.I)  THEN 
ET(I)  -  TT(I) 

ELM(I)  -  FLOAT(NF(I))/TT(I) 

C 

C  WeibuII  components. 

C 

ELSEIF  (TY(I).EQ.2)  THEN 
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DO  75  K  -  1.  NC(I) 

DOFT(K)  -  OFT(I,K) 

75  CONTINUE 

CALL  MLESHAPE(DOFT.NC(I) ,IMAX(1,NF(1)) ,T0L,1.D0,EBETA(I)) 
EBETA(I)  -  BN(NC(1))*EBETA(I) 

SUM  -  0.0 

DO  80  K  -  1,  NF(I) 

SUM  -  SUM  +  OFT(I,K)**EBETA(I) 

80  CONTINUE 

ET(I)  -  FLOAT(NC(I)-NF(I))*TT(I)**EBETA(I)  +  SUM 
ELM(I)  -  FLOAT(NF(I))/ET(I) 

Geometric  components . 

ELSEIF  (TY(I).EQ.3)  THEN 
ET(I)  -  TT(I) 

ELM(I)  -  FLOAT(NF(I))/TT(I) 

ENDIF 

90  CONTINUE 

Determine  the  maximum  failure  rate  estimate 
and  identify  that  component. 

ELMAX(M)  -  0.0 
KEY  -  0 

DO  100  I  -  1,  NCOMP 

IF  (ELM(I) .GT.ELMAX(M))  THEN 
ELMAX(M)  -  ELM(I) 

KEY  -  I 
ENDIF 

100  CONTINUE 

Compute  the  ratios  of  the  failure  rate  estimates  to  their  maximum. 

DO  110  I  -  1,  NCOMP 

ER(I)  -  ELM(I)  /  ELMAX(M) 

110  CONTINUE 

Determination  of  LMU(M) 

SUM  -  0.0 

DO  120  1  -  1,  NCOMP 

SUM  -  SUM  +  (ER(I)*ET(I)) 

120  CONTINUE 

LMU(M)  -  CHISQD(1-ALPHA,NINT(1.3*2*(1+NFC(M))))/(2*SUM) 

C 

C  Compute  estimate  of  overall  reliability  RSL(M)  for  the  system. 
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c 

RSL(M)  -  1.0 
DO  140  J  -  1,  NCS 
PROD  =1.0 

DO  130  I  -  1,  COMP(J,l) 

K  -  COMP(J,I+l) 

IF  (TY(K).EQ.l)  THEN 

PROD  -  PROD*(l  -  SURV(TY(K),LMU(M)*ER(K),EBETA(K),UT(K))) 
ELSEIF  (TY(K).EQ.2)  THEN 

PROD  -  PROD*(l  -  SURV(TY(K) , (Lmj(M)*ER(K) )**(!. /EBETA(K)) , 
*  EBETA(K) ,UT(K))) 

ELSEIF  (TY(K).EQ.3)  THEN 

PROD  -  PROD*(l  -  SURV(TY(K) . l.DO-LMU(M)*ER(K) ,0.D0,UT(K))) 
END  IF 

130  CONTINUE 

REL2(J)  -  1.0  -  PROD 
RSL(M)  -  RSL(M)  *  REL2(J) 

140  CONTINUE 
C 

ENDIF 

C 

C  Case  B  ;  Where  there  are  at  most  ONE  component  which  experienced 
C  ......  failures  during  the  test. 

C 

IF  (FLAG)  THEN 
C 

C  Determine  number  of  complete  systems  (NSYS)  implied  by  the  test 
C  based  on  cut-set  information  after  first  determining  the  number 
C  of  mission  cycles  tested  for  each  component  (to  the  nearest 
C  integer)  (NCYC(I)). 

C 

ISUM  -  0 

DO  150  I  -  1,  NCOMP 

NCYC(I)  -  INT(TT(I)/UT(1)) 

ISUM  -  ISUM  +  NCYC(I) 

ELW(I)  ■=  0 
ER(I)  -  0 
EBETA(I)  -  1.0 
150  CONTINUE 

ELMAX(M)  -  0 
LMU(M)  -  0 
C 

NSYS  -  ISUM 
ISUM  -  0 

DO  170  J  -  1,  NCS 
ISUM  ■>  0 

DO  160  I  -  1,  COMP(J,l) 

ISUM  =  ISUM  +  NCYC(COMP(J,I+l)) 

160  CONTINUE 

ISUM  =  ISUM  /  COMP(J,l) 
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NSYS  -  IMIN(NSYS,ISUM) 

170  CONTINUE 
C 

C  For  zero  failures  in  the  test. 

C 

IF  (ICOUNT.EQ.O)  THEN 

RSL(M)  -  ALPHA** (1.0/FLOAT( NSYS)) 

C 

C  For  failures  experienced  by  a  particular  component  type. 
C 

ELSE 

CALL  GETP( NSYS, ALPHA, RSL(M)) 

ENDIF 

C 

ENDIF 

C 

C  Increment  replication  counter. 

C 

M  -  M  +  1 
C 

ENDDO 

END 


4.  Subroutine  (EVAL) . 

C . 

SUBROUTINE  EVAL 
C 

C  This  subroutine  calls  BUBBLE  to  sort  the  array  RSL(NREP)  in 
C  ascending  order  to  get  an  ordered  array  ORSL(NREP) .  It  also 
C  determine  the  estimate  for  RSLOW  at  the  specified  significance 
C  level  ALPHA  and  the  value  of  LEVEL  in  which  ORSL(LEVEL)  is  closest 
C  to  the  true  reliability  RS . 

C 

C  Include  the  declaration  files 
C  and  declare  the  local  variables. 

C 

INCLUDE  'PARM2  DEF' 

INTEGER  INDEX 
REAL*8  DIFF 
C 

C  Order  the  array  RSL(NREP)  in  ascending  order. 

C 

DO  10  M  =  1,  NREP 
ORSL(M)  -  RSL(M) 

10  CONTINUE 
C 

C  Bubble  Sort.  Sink  the  larger  of  the  pair. 

C 

CALL  BUBBLE(NREP,RSL,ORSL) 

C 
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C  Determine  the  (1 -ALPHA) %  lower  confidence  bound  for  the  system 
C  reliability. 

C 

RSLOW  =  0RSL(NINT(NREP’^'(1 -ALPHA))) 

C 

C  Finding  the  %  confidence  level  for  the  true  reliability  RS . 

C  (ie.  the  proportion  of  RSL(M)  lesser  than  RS) 

C 

DTFF  =  1.0 

INDEX  -  0 

DO  200  M  -  1,  NREP 

IF  (ABS(ORSL(M)-RS) .LT.DIFF)  THEN 
DIFF  -  ABS(ORSL(M)-RS) 

INDEX  -  M 
ENDIF 

200  CONTINUE 
C 

LEVEL  =  FLOAT (INDEX) /NREP 
C 

C  Record  evaluated  parameters  in  RAW2.DAT  (unit  2). 

C 

OPEN(UNIT>=2) 

WRITE(2,300) 

300  F0RMAT(1X,  '  M  LMU(M)  ELMAX(M)  RSL(M) ' , 

*  '  ORSL(M)  NFC(M)’) 

DO  500  M  »  1,  NREP 

WkITE(2,400)  M,1JIU(M) ,ELMAX(M) ,RSL(M) ,0RSL(M) ,NFC(M) 

400  FORMAT ( IX , 16 , 2F12 . 7 , 2F12 .7,110) 

500  CONTINUE 

CLOSE (UNIT=2) 

C 

RETORN 

END 

C . 

C  5.  Report  Generation  Subroutine  (REPORT). 

C . . 

SUBROUTINE  REPORT 
C 

C  This  subroutine  record  the  simulation  results  into  the  ’0UT2.DAT’ 
C  file  as  logic  unit  3. 

C 

C  Include  the  declaration  files 
C  and  declare  local  variables. 

C 

INCLUDE  'PARM2  DEE' 

INTEGER  I.  J,  K.  DUM(IO) 

C 

C  Write  to  output  file  '0UT2.DAT'  designated  as  logic  unit  3. 

C 

0PE:>;a’MT=3) 


c 

WRITE(3,10) 

WRITE(3,20)  NREP 
WRITE(3,25) 

WRITE(3.26) 

WRITE(3,30) 

WRITE(3,40) 

WRITE ( 3 , 50 )  I SEED , NCOMP . ALPHA , TOL . NCS , TCN 
C 

WRITE(3,60) 

DO  200  I  -  1,  NCOMP 

WRITE(3,70)  I,TY(I) ,PARM(1.I) ,PARM(2,I) , UT( I ) , TT( I ) , NC( I ) , NF( I ) 
200  CONTINUE 
WRITE(3,80) 

WRITE(3,90) 

DO  300  I  =  1,  NCOMP 

WRITE(3, 100)  I ,NF(I) , ET ( I ) , ELM( I ) , ER( I ) ,EBETA(I) 

300  CONTINUE 

WRITE(3,110) 

WRITE(3, 120) 

DO  500  J  -  1,  NCS 
DO  400  K  =  1,  10 
DUM(K)  -  CQMP(J.K) 

400  CONTINUE 

WRITE(3.130)  J.DUM.RELKJ)  ,REL2(J) 

500  CONTINUE 

WRITE(3,140) 

WRITE(3, 150)  RS.ELMAX(NREP) ,LMU(NREP) .RSLOW, LEVEL 
C 

10  FORMAT (IX, '0UT2.DAT  :  Output  File  of  the  RETP2  simulation') 

20  FORMAT(lX,'  after  ’,15,'  replications',/) 

25  FORMAT (IX, 'COMMENTS  :  8  COMPONENT  IN  SERIES  ') 

26  FORMATdX,'  DF  -  MINT  (1.3  *  2  *  (1  +  NFC))  ’  ,/) 

30  FORMAT( lx Input  Parameters:',/) 

40  FORMATdX,'  ISEED  NCOMP  ALPHA  TOL  NCS  TCN',/) 

50  FORMATdX.FlO.  1,I8,F8.4,F8.5,2I6,/) 

60  FORMATdX,'  I  ’Ti’d)  PARMl(I)  PARM2  ( I )  UT(I)  TT(I)  NC(I)' 

*  '  NF(I)’,/) 

70  F0RMAT(1X,I3,I6,2F9.5,2F9.2,2I8) 

80  FORMATdX ,/,' Output  Parameters  for  the  LAST  Replication:',/) 

90  FORMATdX,'  I  NF(I)  ET(I)  ELM(I)  ER(I)’ 

*  '  EBETA(I)',/) 

100  FORMAT(lX,I3,I6,E16.7,2F14.7,F14.7) 

110  FORMATdX,/, 'Cut-Set  Data:  '  ./) 

120  FORMATdX,'  J  NUM  Component  List 

*  '  RELl  REL2(M)’,/) 

130  FORMATdX, 13, 15, 913, 2F12. 9) 

140  FORMATdX,/,’  RS  ELMAX(M)  LMU(M)’, 

*  '  RSLOW  LEVEL',/) 

150  FORMATdX,  5F12.  7./) 
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C 

C 


CL0SE(UNIT=3) 

RETURN 

END 


This  portion  of  the  file  contains  functions  and  subroutines 
used  in  the  RETP2  model . 

-  20  Jun  91 

-  by  Yee  Kah-Chee  SMC  2802 


A.  Random  Number  Generating  Subroutine  (LRNDPC) . 

C  (Courtesy  of  Mr.  David  Lim  Hung-Heng) 

C . - . 

SUBROUTINE  LRNDPC  (DSEED.U.N) 

INTEGER  N,  I 

REAL*8  U(N) 

REAL*8  D31M1,  DSEED,  D31 

C  D31M1-2**31  -  1 

C  D31  -2**31 

DATA  D31M1/2147483647 .DO/ 

DATA  D31  /2147483648.D0/ 

DO  5  I-1,N 

C  DSEED  -  DMOD(950706376.DO*DSEED,D31M1) 

DSEED  =  DMOD(16807 .D0*DSEED,D31M1) 

5  U(I)  -  DSEED  /  D31 
RETURN 
END 

C . 

C  B.  Survivability  Function. 

C . 

FUNCTION  SURV(IYTE, PARI ,PAR2, UTIL) 

C 

C  This  function  returns  the  survival  probability  of  the  component  of 
C  different  types  (TYPE)  with  scale  (PARI)  and  shape  (PAR2)  parameters 
C  given  the  specified  utilization  times  or  cycles  (UTIL). 

C 

INTEGER  TYPE,  N 
REAL*8  PARI,  PAR2 ,  UTIL 
C 

IF  (TYPE.EQ.l)  THEN 

SURV  -  EXP(-(PAR1*UTIL)) 

ELSEIF  (TYPE.EQ.2)  THEN 

SURV  -  EXP(-((PAR1*UTIL)**PAR2)) 

ELSE 

N  -  NINT(UTIL) 

SURV  -  PAR1**N 
END  IF 

C 

END 
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C . - . 

C  C.  Bubble  Sort  Routine  in  ASCENDING  Order. 

C . --- . 

SUBROUTINE  BUBBLE(N,LIST,OLIST) 

C 

C  This  subroutine  performs  a  bubble  sort  in  increasing  order  (ie.  sink 
C  the  greater  numeral)  for  the  first  N  terms  in  an  array  LIST  and 
C  returns  the  result  in  OLIST. 

C 

LOGICAL  DONE 
INTEGER  N,  K,  PAIR 
REAL*8  LIST(*),  OLIST(*) 


Sink  the  larger  of  the  pair. 

DO  50  K  -  1,  N 

OLIST(K)  -  LIST(K) 

50  CONTINUE 

PAIR  -  N  -  1 
DONE  -  .FALSE. 

DO  WHILE  (.NOT. DONE) 

DONE  -  .TRUE. 

DO  100  K  -  1,  PAIR 

IF  (OLIST(K) .GT.OLIST(K+l))  THEN 
TEMP  -  OLIST(K) 

OLIST(K)  -  0LIST(K+1) 
OLIST(K+l)  -  TEMP 
DONE  -  .FALSE. 

ENDIF 

100  CONTINUE 

PAIR  -  PAIR  -  1 
ENDDO 
END 


D.  Unbiasing  Factor  for  Biased  MLE  for  Weibull  Shape  Parameter. 


FUNCTION  BN(I) 

This  functon  returns  the  value  of  the  unbiased  factor  for  the  biased 
maximum  likelihood  estimate  of  the  shape  paraiuecer  of  a  Weibull 
distribution  with  a  sample  size  of  N. 

INTEGER  1 
IF  (I.LE.5)  THEN 

BN  -  (I*0.699)/(5.0) 

ELSEIF  (I.EQ.6)  THEN 
BN  -  0.752 

ELSEIF  (I.EQ.7)  THEN 
BN  -  0.786 

ELSEIF  (I.EQ.8)  THEN 
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BN  -  0.82 

ELSEIF  (I.EQ.9)  THEN 
BN  =  0.8395 
ELSEIF  (I.EQ.IO)  THEN 
BN  -  0.859 

ELSEIF  (I.EQ.ll)  THEN 
BN  -  0.871 

ELSEIF  (I.EQ.12)  THEN 
BN  -  0.883 

ELSEIF  (I.EQ.13)  THEN 
BN  -  0.892 

ELSEIF  (I.EQ.14)  THEN 
BN  -  0.901 

ELSEIF  (I.EQ.15)  THEN 
BN  -  0.9075 
ELSEIF  (I.EQ.16)  THEN 
BN  =  0.914 

ELSEIF  (I.EQ.17)  THEN 
BN  -  0.9185 
ELSEIF  (I.EQ.18)  THEN 
BN  -  0.923 

ELSEIF  (I.EQ.19)  THEN 
BN  -  0.927 

ELSEIF  (I.EQ.20)  THEN 
BN  -  0.931 

ELSEIF  (I.LE.25)  THEN 

BN  -  0.931+(I-20)*0. 014/5.0 
ELSEIF  (I.LE.30)  THEN 

BN  -  0.945+(I-25)*0. 01/5.0 
ELSEIF  (I.LE.40)  THEN 

BN  -  0.955+(I-30)*0. 011/10.0 
ELSEIF  (I.LE.60)  THEN 

BN  -  0.966+(I-40)*0. 012/20.0 
ELSEIF  (I.LE.80)  THEN 

BN  -  0.978+(I-60)*0. 006/20.0 
ELSEIF  (I.LE.lOO)  THEN 

BN  -  0.984+(I-80)*0. 003/20.0 
ELSEIF  (I.LE.120)  THEN 

BN  -  0.987+(l-100)*0. 003/20.0 
ELSE 

BN  -  1.0 
ENDIF 
RETURN 
END 

C . - . - . 

C  E.  Biased  MLE  of  Weibull  Shape  Parameter. 

C . . . - . 

SUBROUTINE  MLESHAPE ( T , N , R , DEL. B . BNEW) 

C 

C  This  subroutine  returns  a  biased  estimator  (BNEW)  for  a  Weibull 
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c 

c 

c 

c 

c 

c 


c 


c 

c 


c 


c 


c 


c 

c 

c 

c 

c 

c 

c 


shape  parameter  using  the  Newton-Raphson' s  Method  of  Successive 
Approximation.  The  data  parameters  consist  of  an  ascending  ordered 
list  of  failure  times  (T),  sample  size  (N) ,  number  of  failed  samples 
(R) ,  tolerance  for  convergence  (DEL)  and  an  initial  estimate  of  the 
shape  parameter  (B) . 

LOGICAL  DONE 

INTEGER  N,  R.  I 

REAL*8  GFUNCT,  GPRIME ,  B,  BOLD.  BNEW,  T(*) ,  DEL, 

*  TERMl,  TERM2.  TERM3 .  SUMl ,  SUM2 .  SUM3 ,  SUM4 ,  STEP 

BNEW  -  B 

DONE  -  .FALSE. 


DO  WHILE  ( .NOT. DONE) 

DONE  -  .TRUE. 

TERMl  -  FLOAT(N-R)*(T(R)**BNEW) 

TERM2  -  FLOAT(N-R)*(T(R)**BNEW)*LOG(T(R)) 

TERM3  -  FLOAT(N-R)*(T(R)**BNEW)*LOG(T(R))*LOG(T(R)) 

SUMl  -  0.0 

SUM2  -  0.0 

SUM3  -  0.0 

SUM4  -  0.0 


DO  50  I  -  1,  R 

SUMl  -  SUMl  +  T(I)**BNEW 

SUM2  -  SUM2  +  (T(I)**BNEW)*LOG(T(I)) 

SUM3  -  SUM3  +  (T(I)**BNEW)*LOG(T(I))*LOG(T(I)) 
SUM4  -  SUM4  +  L0G(T(I)) 

50  CONTINUE 


* 


GFUNCT  -  (SUM2+TERM2)/(SUM1+TERM1) 
-  (1.0/ FLOAT (R))*SUM4 


(1.0/BNEW) 


GPRIME  -  (1 .0/(SUMI+TERMl)**2)*(  (SUM1+TERM1)*(SUM3+TERM3) 

*  -  (SUM2+TERM2)**2  ) 

*  +  (1.0/BNEW**2) 


PRINT  60,  GFUNCT, GPRIME, BNEW 

60  FORMAT (IX.  ' GFUNCT  F8 . 3 . '  GPRIME  F8 . 3 , '  BNEW-',F8.3) 


Control  magnitude  of  the  marching  step  towards  convergence 
as  no  more  than  O.I. 


IF  ( (GFUNCT. LT.O)  .AND.  (GPRIME . GT . 0) )  THEN 
STEP  -  VMAX(- .IDO, (GFUNCT/GPRIME)) 

ELSEIF  ( (GFUNCT. GT.O)  .AND.  (GPRIME . LT . 0) )  THEN 
STEP  -  VMAX(- .IDO, (GFUNCT/GPRIME)) 

ELSE 
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STEP  -  VMIN( . IDO, (GFUNCT/GPRIME)) 

ENDIF 

C 

BOLD  -  BNEW 
BNEW  -  BNEW  -  STEP 
C 

C  Check  for  convergence  of  the  MLE  for  the  shape  parameter  B. 
C 

IF  (ABS(BOLD-BNEW) .GT.DEL)  THEN 
DONE  -  .FALSE. 

ENDIF 

C 

C  Avoid  overflow  error  due  to  large  MLE  value  caused  by  small 
C  GPRIME  (slope)  as  GFUNCT  approaches  to  near  zero. 

C  STOP  WHEN  MAGNITUDE  OF  BNEW  EXCEEDS  7. 

C 

IF  (BNEW.GT. 7 .0)  THEN 
BNEW  =  BOLD 
DONE  -  .TRUE. 

ENDIF 

C 

ENDDO 

RETURN 

END 


E.  Chi-Square  Quantile  Function, 

C . 

FUNCTION  CHISQD(P,N) 

C 

C  Modified  version  of  Algorithm  451  from  Comunications  of  the  ACM 
C  Aug  1977  Vol.16  No . 8  . 

C 

C  This  function  evaluates  the  quantile  at  the  probability  level  P 
C  (left  tail  area)  for  the  Chi-square  distribution  with 
C  N  degrees  of  freedom. 

C 

REAL* 8  P 
REAL  X 
INTEGER  IF 

DIMENSION  C(21) ,  A(19) 


DATA 

C/  1.565326E-3, 

★ 

1.060438E-3. 

★ 

-6.950356E-3, 

* 

-1 . 323293E-2 , 

★ 

2.277679E-2, 

★ 

-8.986007E-3, 

* 

-1 . 513904E-2 , 

2.530010E-3, 

-1.450117E-3, 

-k 

5.169654E-3, 
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* 

-1.153761E-2. 

* 

1.128186E-2, 

-k 

2.607083E-2. 

k 

-0.2237368, 

k 

9.780499E-5, 

k 

-8.426812E-4, 

k 

3.125580E-3. 

k 

-8.553069E-3, 

k 

1.348028E-4. 

k 

0.4713941, 

k 

1.0000886  / 

DATA  A/ 

1.264616E-2. 

* 

-1.425296E-2. 

★ 

1.400483E-2, 

k 

-5.886090E-3, 

k 

-1.091214E-2, 

k 

-2.304527E-2, 

k 

3.135411E-3. 

k 

-2.728484E-4, 

k 

-9.699681E-3, 

k 

1.316872E-2. 

k 

2.618914E-2, 

k 

-0.2222222, 

k 

5.406674E-5, 

k 

3.483789E-5. 

k 

-7.274761E-4. 

k 

3.292181E-3. 

k 

-8.729713E-3, 

k 

0.4714045, 

k 

1.  / 

IF  (N-2)  10.  20,  30 
10  CALL  XFROMP( . 5*(1 . -P) ,X,IF) 

CHISQD  -  X 

CHISQD  -  CHISQD*CHISQD 
RETURN 

20  CHISQD  -  -2.*LOG(l. -P) 

RETURN 
30  F  -  N 
FI  -  l./F 

CALL  XFROMP(P,X, IF) 

T  -  X 

F2  -  SQRT(F1)*T 

IF  (N.GE. (2+INT(4,*ABS(T))))  GO  TO  40 

CHISQD  -  (((((((C(1)*F2+C(2))*F2+C(3))*F2+C(4))*F2 

*  +C(5))*F2+C(6))*F2+C(7))*F1+((((((C(8)+C(9)*F2)*F2 

*  +C(10) )*F2+C(11))*F2+C(12))*F2+C(13))*F2+C(14)))*F1+ 

*  (((((C(15)*F2+C(16))*F2+C(17))*F2+C(18))*F2 

*  +C(19))*F2+C(20))*F2+C(21) 

GO  TO  50 

40  CHISQD=( ((A(1)+A(2)*F2)*F1+(((A(3)+A(4)*F2)*F2 


82 


non 


*  +A(5))*F2+A(6)))*F1+(((((A(7)+A(8)*F2)*F2+A(9))*F2 

*  +A(10))*F2+A(11))*F2+A(12)))*F1+(((((A(13)*F2 

*  +A(14))*F2+A(15))*F2+A(16))*F2+A(17))*F2*F2 

*  +A(18))*F2+A(19) 

50  CHISQD  -  CHISQD*CHISQD*CHISQD*F 
RETURN 
END 

C . - . - . 

C  F.  Standard  Normal  Variate  Computation  Subroutine. 

C . - . 

SUBROUTINE  XFROMP(P , X , IFAULT) 

Algorithm  AS  24  J .R. STAT. SOC.  C.  (1969)  Vol.18.  No . 3 . 


non  non  non  n  non 


I FAULT  -  1 

IF  ((P.LE.0.0) .OR. (P.GE.1.0))  GO  TO  100 
I FAULT  -  0 

Get  first  approximation  XO  to  deviate  by  Hastings'  formula 
B  -  P 

IF(B.GT.0.5)  B  -  1.0  -  B 

F  -  -  LOG(B) 

E  -  SQRT(F+F) 

XO  -  -E  +  ((HSTNGS(3)*E+HSTNGS(2))*E+HSTNGS(1))/ 

*  ( ( ( HSTNGS ( 6 ) *E+HSTNGS ( 5 ) ) *E+HSTNGS (4 ) ) *E+1 . 0 ) 

IF  (XO.LT.0.0)  GO  TO  1 
XO  -  0.0 
PO  =  0.5 
XI  -  -RTHFPI 
GO  TO  7 


Find  the  area  PO  corresponding  to  XO 

1  Y  -  XO**2 

IF  (XO.LE. -1.9)  GO  TO  3 
Y  -  -0.5*Y 

(1)  series  approximation 

PO  -  CONNOR(l) 

DO  2  L-2,17 

2  PO  -  P0*Y  +  CONNOR (L) 

PO  -  (PO*Y+1.0)*XO 
XI  -  -  (PO-t-RTHFPI)*EXP( -Y) 

PO  -  P0*RRT2PI  +0.5 
GO  TO  7 

(2)  continued  fraction  approximation 

3  Z  -  1.0 A 
A(2)  -  1.0 
A(3)  -  1.0 
A(4)  -  Z  +  1.0 
A(5)  -  1.0 
W  -  2.0 

C 

4  DO  6  L-1,3.2 
DO  5  J-1,2 

K  -  L  +  J 
KA  -  7  -  K 
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5  A(K)  -  A(KA)  +  A(K)*W*Z 
C 

6  W  -  W  +  1.0 
APPRXU  =  A(2)/A(3) 

APPRXL  -  A(5)/A(4) 

C  -  APPRXU  -  APPRXL 
IF  (C.GE.TERMIN)  GO  TO  4 
XI  -  APPRXL/XO 
PO  -  -X1*RRT2PI*EXP(-0.5*Y) 

C 

C  Get  accurate  value  of  deviate  by  Taylor  Series 
(XI,  X2 ,  X3  are  derivatives  for  the  Taylor  Series 

7  D  -  F  +  LOG(PO) 

X2  -  X0*X1*X1  -XI 
X3  -  Xl**3  +  2.0*XO*X1*X2  -X2 
X  -  ((X3*D/3.0+X2)’^D/2.0+Xl)*D  +  XO 
IF  (P.LE.0.5)  GO  TO  100 
X  -  -X 

100  RETURN 
END 


G.  MAXIMUM  FUNCTIONS. 

C . 

FUNCTION  VMAX(X.Y) 
REAL*8  X.  Y 
IF  (X.GT.Y)  THEN 
VMAX  -  X 
ELSE 

VMAX  -  Y 
ENDIF 
RETURN 
END 

FUNCTION  IMAX(X,Y) 
INTEGER  X.  Y 
IF  (X.GT.Y)  THEN 
IMAX  -  X 
ELSE 

IMAX  =  Y 
ENDIF 
RETURN 
END 

C . 

C  H.  Minimum  Functions. 

C . . 

FUNCTION  VMIN(X,Y) 
REAL*8  X,  Y 
IF  (X.LT.Y)  THEN 
VMIN  -  X 
ELSE 
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VMIN  ■=  Y 
ENDIF 
RETURN 
END 

FUNCTION  IMIN(I.J) 
INTEGER  I ,  J 
IF  (I.LT.J)  THEN 
IMIN  -  I 
ELSE 

IMIN  -  J 
ENDIF 
RETURN 
END 


I.  ROUTINE  TO  EVALUATE  VALUE  OF  P. 

C . - . 

SUBROUTINE  GETP(N . ALPHA , NEWP) 

INTEGER  N 

REAL*8  ALPHA,  OLDP ,  NEWP.  TOL 
LOGICAL  DONE 

OLDP  -  ALPHA** (1.0/FL0AT(N)) 

NEWP  -  1.0 
TOL  -  0.0001 
DONE  -  .FALSE. 

DO  WHILE  ( .NOT. DONE) 

GFUNCT  -  N*OLDP**(N-l)  -  (N- l)*OLDP**N  -  ALPHA 
GPRIME  -  N*(N-l)*OLDP**(N-2)  -  N*(N-1)*0LDP**(N-1) 

NEWP  -  OLDP  -  (GFUNCT/GPRIME) 

IF  ((ABS(NEWP-OLDP) .LE.TOL)  .OR.  (ABS (GFUNCT) . LE . TOL) )  THEN 
DONE  -  .TRUE. 

ENDIF 

OLDP  -  NEWP 
END  DO 
RETURN- 
END 

4.  Program  Output.  (OUT2.DAT) 


The  result  for  the  simulation  run  based  on  the  input  parameters  specified  in 
IN2.DAT  are  computed  and  written  to  the  file  OUT2.DAT.  A  sample  of  this  file  is 
as  follows. 

OUT2.DAT  :  Output  File  of  the  RETP2  simulation 
after  1000  replications 

COMMENTS  :  8  COMPONENT  IN  SERIES 

DF  -  NINT  (1.3  +  2  *  (1  +  NFC)) 
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Input  Parameters; 


ISEED 

NCOMP 

ALPHA 

TOL 

KCS  TCN 

16807.0 

8 

0.2000  0, 

.01000 

8  3 

I 

TY(I) 

PARMl(I) 

PARM2(I) 

UT(I) 

TT(I) 

NC(I) 

NF(I) 

1 

1 

0.00500 

1.00000 

5.00 

5400.00 

20 

36 

2 

1 

0.00500 

1.00000 

5.00 

5400.00 

20 

23 

3 

1 

0.00500 

1.00000 

5.00 

5400.00 

20 

28 

4 

1 

0.00500 

1.00000 

5.00 

5400.00 

20 

25 

5 

2 

0.01000 

2.00000 

15.00 

2700.00 

20 

20 

6 

2 

0.01000 

2.00000 

15.00 

2700.00 

20 

20 

7 

2 

0.01000 

2.00000 

15.00 

2700.00 

20 

20 

8 

2 

0.01000 

2.00000 

15.00 

2700.00 

20 

20 

Output  Parameters  for  the  LAST  Replication: 


I 

NF(I) 

ET(I) 

ELM(I) 

ER(I) 

EBETA(l) 

1 

36 

0.5400000E-(-04 

0.0066667 

1.0000000 

0.0000000 

2 

23 

0. 5400000E+04 

0.0042593 

0.6388889 

0.0000000 

3 

28 

0. 5400000E-H04 

0.0051852 

0.7777778 

0.0000000 

4 

25 

0. 5400000E-t04 

0.0046296 

0.6944444 

0.0000000 

5 

20 

0.2378985E+05 

0.0008407 

0.1261042 

1.5441284 

6 

20 

0.4825114E-t-06 

0.0000414 

0.0062175 

2.1482593 

7 

20 

0.2539978E-f-05 

0.0007874 

0.1181112 

1.5843514 

8 

20 

0.3708000E-^07 

0.0000054 

0.0008091 

2.5186896 

Cut-Set  Data; 

J  NUM  Component  List  RELl  REL?(M) 

1  1100000000  0.975309908  0.955163479 

2  1200000000  0.975309908  0.971117675 

3  1300000000  0.975309908  0.964950144 

4  14C0000000  0.975309908  0.968645990 

5  1500000000  0.977751195  0.927053154 

6  1600000000  0.977751195  0.981007099 

7  1700000000  0.977751195  0.923940539 

8  1800000000  0.977751195  0.993218899 

RS  ELMAX(M)  LMU(M)  RSLOW  LEVEL 

0.8269590  0.0066667  0.0091745  0.7754698  0.9940000 
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APPENDIX  D  :  Evaluation  of  Subroutines  and  Functions 


RANDOM  NUMBER  GENERATOR 
(LRNDPC) 

Evaluation 

One  thousand  unifonn  random  real  numbers  between  0  and  1  are  generated 
using  the  random  number  generating  routine  LRNDPC.  From  these  uniformly 
distributed  numbers,  1000  exponential  (with  scale  parameter  1)  numbers  and  1000 
Weihull  (with  scale  parameter  1  and  shape  parameter  2)  numbers  were  generated. 

Uniform  Random  Variate 


Plotof  U(N+1)  vs  U(N) 
fof  1000  Uniform  RVs  from  LRNDPC 

1  •  ,  .  ... 

o.«  ■ 

OJ  -  .  ■ 

07 

o.« 

^  0-a 

D 

04 

07 

07 

0  1 

0 

0  0.1  07  03  a4  OB  0.0  07  OB  a»  1 

U(N) 

Figure  1  :  Uniform  RVs  generated  by  LRNDPC 


Figure  1  above  shows  a  plot  of  1000  unifonn  real  numbers  against  their 
predecessors.  The  uniformity  of  the  distribution  of  points  over  the  state  space 
confirms  LRNDPC's  adequacy  in  generating  unifonn  random  numbers. 


8<S 


Exponential  and  Weibull  real  numbers  were  generated  using  these  1000 
Uniform(0,l)  random  numbers.  The  cumulative  histograms  of  these  resultant 
random  variates  were  compared  with  their  respective  theoretical  cumulative 
distribution  functions  (cdfs). 


Exponential  Random  Variate 


Distribution  Fitting 

for  1000  Exp(1)  RV«  by  LRNOPC 


Figure  2  :  Exponential  RVs  generated  by  LRNDPC 


Figure  2  shows  the  close  distribution  fit  between  the  theoretical  cdf  (line)  and 
the  cumulative  distribution  of  exponential  RV  generated  using  LRNDPC. 


F{t)  =  exp(->t) 

.-.  t  =  -IlnlFC/)} 
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Weibull  Random  Variate 


Distribution  Fitting 

for  1000  Wal(1,2)  RVt  by  LnND'>C 
(la.  Scala-1,  Shapa-Q 


Value  of  Random  Variable 


Figure  3  :  Weibull  RVs  generated  by  LRNDPC 


Figure  3  shows  the  close  distribution  fit  between  the  theoretical  cdf  (line)  and 
the  cumulative  distribution  of  Weibull  RV  generated  using  LRNDPC. 


F{t)  =  exp{-(Xr)^} 

=  -i[ln{F(0)P 


t 


Plot  of  Unbiasing  Factor  B(N)  vs  N 
for  WeibuII  Shape  Parameter  Estimation 


N  =  Test  Sample  Size 
B(N)  =  Unbiasing  Factor  for  MLE 


Plot  of  BN  vs  N 

(Weibult  Shape  Parameter  Estimation) 


0  20  40  «0  so  1M  120 

Sampie  Size  N 


Figure  1  :  B(N)  vs  N 


The  function  BN(N)  returns  the  linear-interpolated  values  of  the  unbiasing 
factor  for  the  raw  MLE  ^  for  both  RETPl  and  RETP2. 


Evaluation  of  Subroutine  CHISQD  and  XFROMP 

The  x"  statistics  for  1  to  499  degrees  of  freedom  for  a  values  of  0.1  and  0.2  are 
generated  using  the  routines  CHISQD  and  XFROMP.  These  outputs  matched  those 
tabulated  in  the  mathematical  tables  of  any  general  textbook  on  statistics. 
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Prob  [  Xd/  ^  table  value  ]  =  1  -  Of  =  0.9 


df 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

2. 

71 

4. 

61 

6. 

25 

7. 

78 

9. 

2A 

10. 

6A 

12. 

02 

13. 

36 

lA. 

68 

1 

15. 

99 

17. 

27 

18. 

55 

19. 

81 

21. 

06 

22. 

31 

23. 

5A 

2A. 

77 

25. 

99 

27. 

20 

2 

28. 

A1 

29. 

62 

30. 

81 

32. 

01 

33. 

20 

3A. 

36 

35. 

56 

36. 

74 

37. 

92 

39. 

09 

3 

AO. 

26 

Al. 

A2 

A2. 

58 

A3. 

75 

AA. 

90 

A6. 

06 

A7. 

21 

A8. 

36 

Ag. 

51 

50. 

66 

A 

51. 

81 

52. 

95 

5A. 

09 

55. 

23 

56. 

37 

57. 

51 

58. 

6A 

59. 

77 

60. 

91 
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Table  lA  : 


8  Exp  in  Series,  RS  =  0.931  (Hi) 

min  \  =  0.0002  f/hr,  max  X  =  0.0016  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

■  -■ 

Measures  of  Accuracy 

RSLOW 

LEVEL 

1 

Test  5  until 

5  failed. 

NFC=40 

2*NFC 

(80) 

0.1 

0.919 

0.982 

0.2 

0.919 

0.960 

2*(NFC  + 
NCOMP) 

(96) 

0.1 

0.906 

1.000 

0.2 

0.905 

0.999 

2*NFC- 

NCOMP 

(72) 

0.1 

0.927 

0.949 

0.2 

0.927 

0.880 

2*(NFC- 

NCOMP) 

(64) 

0.1 

0.934 

0.821 

0.2 

0.934 

0.702 

2 

Test  15  until 

15  failed. 

NFC  =120 

2*NFC 

(240) 

0.1 

0.928 

0.955 

0.2 

0.927 

0.908 

2*(NFC+ 

NCOMP) 

(256) 

0.1 

0.923 

0.990 

0.2 

0.923 

0.975 

2*NFC- 

NCOMP 

(232) 

0.1 

0.930 

0.916 

0.2 

0.930 

0.833 

2*(NFC- 

NCOMP) 

(224) 

0.1 

0.932 

0.844 

0.2 

0.932 

0.747 

3 

Test  15  until 

11  failed. 

NFC  =  88 

2*NFC 

(176) 

0.1 

0.927 

0.955 

0.2 

0.926 

0.916 

2*(NFC+ 

NCOMP) 

(192) 

0.1 

0.921 

0.996 

0.2 

0.920 

0.988 

2*NFC- 

NCOMP 

(168) 

0.1 

0.930 

0.916 

0.2 

0.929 

0.843 

2*(NFC- 

NCOMP) 

(160) 

0.1 

0.933 

0.843 

0.2 

0.932 

0.735 
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Table  lA  : 


8  Exp  in  Series,  RS  =  0.931  (Hi)  (Cont..,) 

min  A  =  0.0002  f/hr,  max  A  =  0.0016  i/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC=56 

2*NFC 

(112) 

0.1 

0.924 

0.970 

0.2 

0.923 

0.931 

2*(NFC+ 

NCOMP) 

(128) 

0.1 

0.915 

0.998 

0.2 

0.913 

0.994 

2*NFC- 

NCOMP 

(104) 

0.1 

0.929 

0.919 

0.2 

0.928 

0.853 

2*(NFC- 

NCOMP) 

(%) 

0.1 

0.934 

0.835 

0.2 

0.933 

0.720 

5 

Test  15  until 

3  failed. 

NFC =24 

2*NFC 

(48) 

0.1 

0.915 

0.986 

0.2 

0.912 

0.975 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.891 

1.000 

0.2 

0.888 

1.000 

2*NFC- 

NCOMP 

(40) 

0.1 

0.927 

0.944 

0.2 

0.926 

0.860 

WmmMm 

0.1 

0.939 

0.753 

0.2 

0.939 

0.634 
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Table  IB  : 


8  Exp  in  Series,  RS  =  0.803  (Lo) 

min  \  =  0.001  f/hr,  max  \  =  0.0045  f/hr,  UT  =  10  hrs 


Test 

Deg  of 

S/N 

Plan 

Freedom 

Measures  of  Accuracy 


Test  5  until 
5  failed. 

NFC=40 


Test  15  until 
15  failed. 

NFC  =  120 


Test  15  until 
11  failed. 

NFC=88 


2*NFC 

(80) 


2*(NFC  + 
NCOMP) 
(96) 


2*NFC- 

NCOMP 

(72) 


2*(NFC- 

NCOMP) 

(64) 


2*NFC 

(240) 


2*(NFC+ 

NCOMP) 

(256) 


2*NFC- 

NCOMP 

(232) 


2*(NFC- 

NCOMP) 

(224) 


2*NFC 

(176) 


2*(NFC+ 

NCOMP) 

(192) 


2*NFC- 

NCOMP 

(168) 


a 

RSLOW 

LEVEL 

0.1 

0.773 

0.986 

0.2 

0.773 

0.963 

0.1 

0.738 

1.000 

0.2 

0.736 

0.999 

0.1 

0.792 

0.953 

0.2 

0.792 

0.892 

0.1 

0.811 

0.831 

0.2 

0.812 

0.713 

0.1 

0.794 

0.%2 

0.2 

0.793 

0.916 

0.1 

0.783 

0.993 

0.2 

0.782 

0.981 

0.1 

0.800 

0.924 

0.2 

0.799 

0.840 

0.1 

0.806 

0.858 

0.2 

0.805 

0.755 

0.1 

0.792 

0.966 

0.2 

0.790 

0.921 

0.1 

0.776 

0.996 

0.2 

0.774 

0.989 

0.1 

0.800 

0.922 

0.2 

0.798 

0.855 

0.1 

0.808 

0.855 

0.2 

0.806 

0.746 
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Table  IB  : 


8  Exp  in  Series,  RS  =  0.803  (Lo)  (Cont...) 

min  X  =  0.001  f/hr,  max  X  =  0.0045  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy  || 

— 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC  =  56 

2*NFC 

(112) 

0.1 

0.785 

0.974 

0.2 

0.782 

0.938 

2*(NFC+ 

NCOMP) 

(128) 

0.1 

0.760 

0.998 

0.2 

0.756 

0.996 

2*NFC- 

NCOMP 

(104) 

0.1 

0.798 

0.925 

0.2 

0.795 

0.860 

2*(NFC- 

NCOMP) 

(96) 

0.1 

0.811 

0.841 

0.2 

0.808 

0.727 

5 

Test  15  until 

3  failed. 

NFC  =  24 

2*NFC 

(48) 

0.1 

0.759 

0.989 

0.2 

0.755 

0.977 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.700 

1.000 

0.2 

0.793 

1.000 

2*NFC. 

NCOMP 

(40) 

0.1 

0.791 

0.949 

0.2 

0.789 

0.872 

0.1 

0.825 

0.763 

0.2 

0.825 

0.642 
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Table  2A  : 


8  Wei  in  Series,  RS  =  0.980  (Hi) 

min  X  =  0.001  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


Test 

Deg  of 

Measures  of  Accuracy  | 

S/N 

Plan 

Freedom 

a 

RSLOW 

LEVEL 

1 

Test  5  until 

2*NFC 

0.1 

0.947 

0.992 

5  failed. 

(80) 

0.2 

0.930 

0.989 

NFC=40 

2*(NFC+ 

0.1 

0.937 

0.994 

NCOMP) 

(96) 

0.2 

0.918 

0.993 

2*NFC- 

0.1 

0.951 

0.989 

NCOMP 

(72) 

0.2 

0.937 

0.986 

2*(NFC- 

0.1 

0.956 

0.985 

NCOMP) 

(64) 

C.2 

0.943 

0.981 

2 

Test  15  until 

2*NFC 

0.1 

0.978 

0.918 

15  failed. 

(240) 

0.2 

0.974 

0.913 

NFC =120 

2*(NFC  + 

0.1 

0.977 

0.931 

NCOMP) 

(256) 

0.2 

0.972 

0.924 

2*NFC- 

0.1 

0.979 

0.914 

NCOMP 

(232) 

0.2 

0.975 

0.901 

2*(NFC- 

0.1 

0.980 

0.904 

NCOMP) 

(224) 

0.2 

0.975 

0.889 

3 

Test  15  until 

2*NFC 

0.1 

0.982 

0.876 

11  failed. 

(176) 

0.2 

0.977 

0.860 

NFC  =  88 

2*(NFC+ 

0.1 

0.980 

0.894 

NCOMP) 

(192) 

0.2 

0.975 

0.882 

2*NFC- 

0.1 

0.983 

0.861 

NCOMP 

(168) 

0.2 

0.978 

0.839 

2*(NFC- 

0.1 

0.983 

0.840 

NCOMP) 

(160) 

0.2 

0.979 

0.819 
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Table  2A  :  8  Wei  in  Series,  RS  =  0.980  (Hi)  (Cont...) 

min  \  =  0.001  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC  =  56 

2*NFC 

(112) 

0.1 

0.987 

0.800 

0.2 

0.981 

0.779 

2*(NFC  + 
NCOMP) 
(128) 

0.1 

0.985 

0.839 

0.2 

0.978 

0.824 

2*NFC- 

NCOMP 

(104) 

0.1 

0.988 

0.776 

0.2 

0.982 

0.753 

2*(NFC- 

NCOMP) 

(%) 

0.1 

0.989 

0.746 

0.2 

0.983 

0.732 

5 

Test  15  until 

3  failed. 

NFC  =  24 

2*NFC 

(48) 

0.1 

0.994 

0.621 

0.2 

0.991 

0.584 

2*(NFC  + 
NCOMP) 

(64) 

0.1 

0.993 

0.705 

0.2 

0.988 

0.685 

2*NFC- 

NCOMP 

(40) 

0.1 

0.995 

0.548 

0.2 

0.992 

0.514 

mmsm 

0.1 

0.996 

0.468 

0.2 

0.993 

0.417 
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Table  2B  : 


8  Wei  in  Series,  RS  =  0.832  (Lo) 

min  X  =  0.003  f/hr,  max  X  =  0.024  f/hr,  UT  =  10  hrs 


Deg  of 
Freedom 


Measures  of  Accuracy 


Test  5  until 
5  failed. 

NFC=40 


Test  15  until 
IS  failed. 

NFC  =  120 


Test  15  until 
11  failed. 

NFC=88 


2*NFC 

(80) 


2*(NFC+ 

NCOMP) 

(96) 


2*NFC- 

NCOMP 

(72) 


2*(NFC- 

NCOMP) 

(64) 


2*NFC 

(240) 


2*(NFC+ 

NCOMP) 

(256) 


2*NFC- 

NCOMP 

(232) 


2»(NFC- 

NCOMP) 

(224) 


2*NFC 

(176) 


2*(NFC  + 
NCOMP) 
(192) 


2*NFC- 

NCOMP 

(168) 


a 

RSLOW 

LEVEL 

0.1 

0.736 

0.985 

0.2 

0.688 

0.983 

0.1 

0.696 

0.992 

0.2 

0.641 

0.991 

0.1 

0.757 

0.980 

0.2 

0.713 

0.973 

0.1 

0.778 

0.968 

0.2 

0.739 

0.954 

0.1 

0.834 

0.895 

0.2 

0.812 

0.876 

0.1 

0.825 

0.920 

0.2 

0.801 

0.904 

0.1 

0.839 

0.880 

0.2 

0.817 

0.858 

0.1 

0.844 

0.866 

0.2 

0.823 

0.838 

0.1 

0.854 

0.838 

0.2 

0.831 

0.808 

0.1 

0.842 

0.882 

0.2 

0.817 

0.861 

0.1 

0.860 

0.809 

0.2 

0.837 

0.776 

0.1 

0.865 

0.777 

0.844 

0.734 
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Table  2B  : 


8  Wei  in  Series,  RS  =  0.832  (Lo)  (Cont„.) 

min  \  =  0.003  f/hr,  max  X  =  0.024  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC=56 

2*NFC 

(112) 

01 

0.878 

0.738 

0.2 

0.857 

0.708 

2*(NFC+ 

NCOMP) 

(128) 

0.1 

0.864 

0.804 

0.2 

0.839 

0.780 

2*NFC- 

NCOMP 

(104) 

0.1 

0.886 

0.702 

0.2 

0.866 

0.655 

2*(NFC- 

NCOMP) 

(96) 

0.1 

0.894 

0.642 

0.2 

0.875 

0.587 

5 

Test  15  until 

3  failed. 

NFC =24 

2*NFC 

(48) 

0.1 

0.910 

0.611 

0.2 

0.888 

0.560 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.885 

0.752 

0.2 

0.856 

0.719 

2*NFC- 

NCOMP 

(40) 

0.1 

0.923 

0.515 

0.2 

0.905 

0.462 

■wM 

0.1 

0.936 

0390 

0.2 

0.922 

0.306 
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Table  3A  :  4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.980  (Hi) 

min  X  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


Test 

Deg  of 

S/N 

Plan 

Freedom 

Measures  of  Accuracy 


Test  5  until 
5  failed. 

NFC  =  40 


Test  15  until 
15  failed. 

NFC =120 


Test  15  until 
11  failed. 

NFC  =  88 


2*NFC 

(80) 


2*(NFC+ 

NCOMP) 

(96) 


2*NFC- 

NCOMP 

(72) 


2*(NFC- 

NCOMP) 

(64) 


2*NFC 

(240) 


2*(NFC+ 

NCOMP) 

(256) 


2*NFC- 

NCOMP 

(232) 


2*(NFC- 

NCOMP) 

(224) 


2*NFC 

(176) 


2*(NFC+ 

NCOMP) 

(192) 


2*NFC- 

NCOMP 

(168) 


a 

RSLOW 

LEVEL 

0.1 

0.979 

0.942 

0.2 

0.978 

0.905 

0.1 

0.975 

0.987 

0.2 

0.974 

0.976 

0.1 

0.981 

0.881 

0.2 

0.980 

0.805 

0.1 

0.983 

0.771 

0.2 

0.982 

0.684 

0.1 

0.981 

0.863 

0.2 

0.980 

0.800 

0.1 

0.979 

0.941 

0.2 

.  . 

0.979 

0.898 

0.1 

0.981 

0.881 

0.2 

0.980 

0.805 

0.1 

0.982 

0.725 

0.2 

0.981 

0.631 

0.1 

0.981 

0.864 

0.2 

0.980 

0.801 

0.1 

0.979 

0.951 

0.2 

0.978 

0.907 

0.1 

0.982 

0.802 

0.2 

0.981 

0.698 

0.1 

0.982 

0.702 

0.982 

0.591 
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Table  3A  : 


4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.980  (Hi)  (Com..,) 
min  X  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

Q 

Measures  of  Accuracy 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC=56 

2*NFC 

(112) 

0.1 

0.981 

0.865 

0.2 

0.980 

0.787 

2*(NFC+ 

NCOMP) 

(128) 

0.1 

0.978 

0.952 

0.2 

0.978 

0^20 

2*NFC- 

NCOMP 

(104) 

0.1 

0.982 

0.769 

0.2 

0.982 

0.676 

2*(NFC- 

NCOMP) 

(96) 

0.1 

0.983 

0.644 

0.2 

0.983 

0.523 

5 

Test  15  until 

3  failed. 

NFC =24 

2*NFC 

(48) 

0.1 

0.982 

0.843 

0.2 

0.981 

0.762 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.976 

0.970 

0.2 

0.975 

0.941 

2*NFC- 

NCOMP 

(40) 

0.1 

0.984 

0.684 

0.2 

0.984 

0.580 

0.1 

0.987 

0.459 

0.2 

0.987 

0.356 
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Table  3B  : 


4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.809  (Lo) 
min  X  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


Test  15  until 
15  failed. 


Test  15  until 
11  failed. 

NFC =88 


Deg  of 

Measures  of  Accuracy 

Freedom 

RSLOW 

LEVEL 

2*NFC 

0.1 

0.788 

0.%1 

(80) 

0.2 

0.777 

0.930 

2*(NFC+ 

0.1 

0.754 

0.996 

NCOMP) 

(96) 

0.2 

0.741 

0.987 

2*NFC- 

0.1 

0.805 

0.910 

NCOMP 

(72) 

0.2 

0.796 

0.874 

2*(NFC- 

0.1 

0.823 

0.838 

NCOMP) 

(64) 

0.2 

0.815 

0.741 

2*NFC 

0.1 

0.808 

0.909 

(240) 

0.2 

0.805 

0.842 

2*(NFC+ 

0.1 

0.797 

0.%2 

NCOMP) 

(256) 

0.2 

0.794 

0.929 

2*NFC- 

0.1 

0.813 

0.854 

NCOMP 

(232) 

0.2 

0.811 

0.778 

2*(NFC- 

0.1 

0.819 

0.787 

NCOMP) 

(224) 

0.2 

0.817 

0.717 

2*NFC 

0.1 

0.811 

0.885 

(176) 

0.2 

0.807 

0.820 

2*(NFC  + 

0.1 

0.797 

0.%2 

NCOMP) 

(192) 

0.2 

0.792 

0.925 

2*NFC- 

0.1 

0.818 

0.821 

NCOMP 

(168) 

0.2 

0.814 

0.737 

2*(NFC- 

0.1 

0.826 

0.741 

NCOMP) 

(160) 

0.2 

0.822 

0.647 

Table  3B  : 


4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.809  (Lo)  (Cont...) 
min  X  =  0.002  f/hr,  max  X  =  0.008  f/hr,  UT  =  10  hrs 


S/N 

Test 

Plan 

Deg  of 
Freedom 

a 

Measures  of  Accuracy  | 

RSLOW 

LEVEL 

4 

Test  15  until 

7  failed. 

NFC =56 

2*NFC 

(112) 

0.1 

0.814 

0.872 

0.2 

0.810 

0.792 

2*(NFC  + 
NCOMP) 
(128) 

0.1 

0.792 

0.%3 

0.2 

0.787 

0.931 

2*NFC- 

NCOMP 

(104) 

0.1 

0.825 

0.775 

0.822 

0.685 

2*(NFC- 

NCOMP) 

(96) 

0.1 

0.836 

0.663 

0.2 

0.834 

0.550 

5 

Test  15  until 

3  failed. 

NFC =24 

2*NFC 

(48) 

0.1 

0.825 

0.836 

0.2 

0.815 

0.755 

2*(NFC+ 

NCOMP) 

(64) 

0.1 

0.780 

0.970 

0.2 

0.766 

0.940 

2*NFC- 

NCOMP 

(40) 

0.1 

0.849 

0.679 

0.2 

0.842 

0.553 

0.1 

0.874 

0.437 

0.2 

0.869 

0.338 
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APPENDIX  F  :  Taburated  Run  Results  for  RETP2 
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Table  4 A  :  8  Exp  in  Series,  RS  =  0.961  (Hi) 

A  =  0.001  f/hr,  UT  =  5  hrs 


S/N 

Degrees 

of 

Freedom 

K  /  EfNFCJ 

Measures  of  Accuracy 

(TT) 

a 

RSLOW 

LEVEL 

1 

2*(1  +  NFC) 

0.25  /  1.2 

0.1 

0.950 

0.851 

(225) 

0.2 

0.935 

0.851 

0.1 

0.957 

0.857 

0.2 

0.954 

0.857 

1.0  /  4.8 

0.1 

0.957 

0.941 

(900) 

0.2 

0.957 

0.850 

2.0  /  9.6 

0.1 

0.958 

0.916 

(1800) 

0.2 

0.960 

0.850 

3.0  /  14.4 
(2700) 

0.1 

0.959 

0.916 

0.2 

0.959 

0.809 

4.0  /  19.2 
(3600) 

0.1 

0.959 

0.937 

0.2 

0.960 

0.843 

5.0  /  24 

0.1 

0.960 

0.926 

(4500) 

0.2 

0.960 

0.814 

10.0  /  48 

0.1 

0.960 

0.924 

(9000) 

0.2 

0.960 

0.809 

20.0  /  96 

0.1 

0.960 

0.914 

(18000) 

0.2 

0.961 

0.820 

30.0  /  144 

0.1 

0.961 

0.906 

(27000) 

0.2 

0.961 

0.804 
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Table  4 A  ; 


8  Exp  in  Series,  RS  =  0.961  (Hi)  (Cont...) 
X  =  0.001  f/hr,  UT  =  5  hrs 


S/N 

Degrees 

of 

Freedom 

K  /  E[NFC1 
(TT) 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

2 

1.3* 

2*(1+NFC) 

0.25  /  1.2 
(225) 

0.1 

0.950 

0.851 

0.2 

0.935 

0.851 

0.1 

0.957 

0.857 

0.2 

0.941 

0.857 

1.0  /  4.8 
(900) 

0.1 

0.946 

0.981 

0.2 

0.945 

0.941 

2.0  /  9.6 
(1800) 

0.1 

0.947 

0.989 

0.2 

0.950 

0.948 

3.0  /  14.4 
(2700) 

0.1 

0.949 

0.997 

0.2 

0.949 

0.969 

4.0  /  19.2 
(3600) 

0.1 

0.948 

0.998 

0.2 

0.950 

0.987 

5.0  /  24 
(4500) 

0.1 

0.949 

0.998 

0.2 

0.949 

0.985 

10.0  /  48 
(9000) 

0.1 

0.949 

1.000 

0.2 

0.949 

0.997 

20.0  /  96 
(18000) 

0.1 

0.949 

1.000 

0.2 

0.950 

1.000 

30.0  /  144 
(27000) 

0.1 

0.950 

1.000 

0.2 

0.949 

1.000 
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Table  4B  :  8  Exp  in  Series,  RS  =  0.819  (Lo) 

X  =  0.005  f/hr,  UT  =  5  hrs 


S/N 

Degrees 

of 

Freedom 

K  /  E(NFC] 

Measures  of  Accuracy 

(TT) 

a 

RSLOW 

LEVEL 

1 

2*(1  +  NFC) 

0.25  /  1.2 

0.1 

0.774 

0.851 

(45) 

0.2 

0.702 

0.851 

0.1 

0.801 

0.857 

Wmm 

0.2 

0.788 

0.857 

1.0  /  4.8 

0.1 

0.801 

0.941 

(180) 

0.2 

0.803 

0.850 

2.0  /  9.6 

0.1 

0.807 

0.916 

(360) 

0.2 

0.814 

0.847 

3.0  /  14.4 
(540) 

0.1 

0.812 

0.916 

0.2 

0.812 

0.809 

4.0  /  19.2 
(720) 

0.1 

0.809 

0.923 

0.2 

0.817 

0.840 

5.0  /  24 

0.1 

0.813 

0.925 

0.2 

0.816 

0.814 

10.0  /  48 

0.1 

0.816 

0.919 

(1800) 

0.2 

0.817 

0.809 

20.0/% 

0.1 

0.817 

0.814 

(3600) 

0.2 

0.818 

0.820 

30.0  /  144 

0.1 

0.819 

0.907 

- 1 

(5400) 

0.2 

-  -  - 

0.818 

0.804 

no 


Table  4B  ; 


8  Exp  in  Series,  RS  =  0.819  (Lo)  (Cont...) 
A  =  0.005  f/hr,  UT  =  5  hrs 


S/N 

Degrees 

of 

Freedom 

K  /  E[NFC1 
(IT) 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

2 

1.3* 

2*(1  +  NFC) 

0.25  /  1.2 
(45) 

0.1 

0.774 

0.851 

0.2 

0.702 

0.851 

0.5  /  2.4 
(90) 

0.1 

0.801 

0.857 

0.2 

0.736 

0.857 

1.0  /  4.8 
(180) 

0.1 

0.759 

0.981 

0.2 

0.753 

0.941 

2.0  /  9.6 
(360) 

0.1 

0.762 

0.989 

0.2 

0.771 

0.948 

3.0  /  14.4 
(540) 

0.1 

0.770 

0.991 

0.2 

0.768 

0.969 

4.0  /  19.2 
(720) 

0.1 

0.766 

0.997 

0.2 

0.772 

0.987 

5.0  /  24 
(900) 

0.1 

0.769 

0.998 

0.2 

0.771 

0.985 

10.0  /  48 
(1800) 

0.1 

0.771 

1.000 

0.2 

0.771 

0.997 

20.0  /  96 
(3600) 

0.1 

0.771 

1.000 

0.2 

0.772 

1.000 

30.0  /  144 
(5400) 

0.1 

0.773 

1.000 

0.2 

0.772 

1.000 
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Table  5A  : 


8  Wei  in  Series,  RS  =  0.956  (Hi)  (*) 
A  =  0.005  f/hr,  UT  =  15  hrs 


S/N 

Degrees 

of 

Freedom 

K  /  EfNFCJ 
(TT) 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

1 

2*(1  +  NFC) 

0.25  /  1.2 
(45) 

0.1 

1.000 

0.186 

0.2 

1.000 

0.158 

0.5  /  2.4 
(90) 

0.1 

0.986 

0301 

0.2 

0.979 

0.458 

1.0  /  4.8 
(180) 

0.1 

0.%7 

0.767 

0.2 

0.960 

0.732 

2.0  /  9.6 
(360) 

0.1 

0.957 

179 

0.2 

0.952 

0.854 

3.0  /  14.4 
(540) 

0.1 

0.952 

0.934 

0.2 

0.946 

0.922 

4.0  /  19.2 
(720) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

5.0  /  24 
(900) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

10.0  /  48 
(1800) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

20.0  /  % 
(3600) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

30.0  /  144 
(5400) 

0.1 

0.952 

0.940 

0.2 

0.946 

0.928 

(•)  20  test  items  for  each  Weibull  component. 
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Table  5A  : 


8  Wei  in  Series,  RS  =  0.956  (Hi)  (*)  (Cont...) 
X  =  0.005  f/hr,  UT  =  15  hrs 


Degrees 

of 

Freedom 

K  /  ElNFCl 
(TT) 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

1.3* 

2*(1  +  NFC) 

0.25  /  1.2 
(45) 

0.1 

1.000 

0.258 

0.2 

0.999 

0.224 

0.5  /  2.4 
(90) 

0.1 

0.983 

0.635 

0.2 

0.973 

0.593 

190  /  4.8 
(180) 

0.1 

0.958 

0.884 

0.2 

0.949 

0.866 

2.0  /  9.6 
(360) 

0.1 

0.946 

0.963 

0.2 

' 

0.939 

0.956 

3.0  /  14.4 
(540) 

0.1 

0.939 

0.984 

0.2 

0.930 

0.976 

4.0  /  19.2 
(720) 

0.1 

0.938 

0.987 

0.2 

0.930 

0.981 

5.0  /  24 
(900) 

0.1 

0.938 

0.987 

0.2 

0.930 

0.981 

10.0  /  48 
(1800) 

0.1 

0.938 

0.987 

0.2 

0.930 

0.981 

20.0  /  % 
(3600) 

0.1 

0.938 

0.987 

0.2 

0.930 

0.981 

30.0  /  144 
(5400) 

0.1 

0.938 

0.987 

0.2 

0.930 

0.981 
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Table  5B  : 


8  Wei  in  Series,  RS  =  0.835  (Lo)  (*) 
X  =  0.01  f/hr,  UT  =  15  hrs 


K  /  E[NFC1 
(TT) 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

0.25  /  1.2 
(22.5) 

0.1 

0.966 

0.222 

0.2 

0.963 

0.161 

0.5  /  2.4 
(45) 

0.1 

0.932 

0.367 

0.2 

0.913 

0.319 

1.0  /  4.8 
(90) 

0.1 

0.874 

0.704 

0.2 

0.858 

0.658 

2.0  /  9.6 
(180) 

0.1 

0.842 

0.851 

0.2 

0.832 

0.819 

3.0  /  14.4 
(270) 

0.1 

0.829 

0.924 

0.2 

0.814 

0.902 

4,0  /  19.2 
(360) 

0.1 

0.827 

0.928 

0.2 

0.813 

0.908 

5.0  /  24 
(450) 

0.1 

0.827 

0.928 

0.2 

0.813 

0.908 

10.0  /  48 
(900) 

0.1 

0.827 

0.928 

0.2 

0.813 

0.908 

20.0  /  96 
(1800) 

0.1 

0.827 

0.928 

0.2 

0.813 

0.908 

30.0  /  144 

0.1 

0.827 

0.928 

(2700) 


Table  5B  : 


8  Wei  in  Series,  RS  =  0.835  (Lo)  (*) 
X  =  0.01  f/hr,  UT  =  15  hrs 


Degrees 

of 

Freedom 

K  /  EINFCJ 
(IT) 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

1.3* 

2*(1  +  NFC) 

0.25  /  1.2 
(22.5) 

0.1 

0.957 

0.326 

0.2 

0.953 

0.268 

0.5  /  2.4 
(45) 

0.1 

0.914 

0.572 

0.2 

0.890 

0.522 

1.0  /  4.8 
(90) 

0.1 

0.842 

0.881 

0.2 

0.821 

0.862 

— 

0.1 

0.802 

0.968 

0.2 

0.788 

0.962 

3.0  /  14.4 
(270) 

0.1 

0.786 

0.986 

0.2 

0.766 

0.983 

4.0  /  19.2 
(360) 

0.1 

0.784 

0.989 

0.2 

0.766 

0.986 

0.1 

0.784 

0.989 

0.2 

0.766 

0.986 

10.0  /  48 
(900) 

0.1 

0.784 

0.989 

0.2 

0.766 

0.986 

20.0/96 

(1800) 

0.1 

0.784 

0.989 

0.2 

0.766 

0.986 

30.0  /  144 
(2700) 

0.1 

0.784 

0.989 

0.2 

0.766 

0.986 
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Table  6A  :  4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.958  (Hi)  (*) 

A(exp)  =  0.001  f/hr,  UT(exp)  =  5  hrs 
X(wei)  =  0.005  f/hr,  UT(wei)  =  15  hrs 


S/N 

Degrees 

of 

Freedom 

Measures  of  Accuracy 

a 

RSLOW 

LEVEL 

1 

2*(1  +  NFC) 

0.25  /  1.2 
(225) 

(45) 

0.1 

1.000 

0.620 

0.2 

0.995 

0.451 

0.5  /  2.4 
(450) 

(90) 

0.1 

0.982 

0.623 

0.2 

0.975 

0.573 

1.0  /  4.8 
(900) 

(180) 

0.1 

0.971 

0.736 

0.2 

0.965 

0.684 

2.0  /  9.6 
(1800) 

(360) 

0.1 

0.964 

0.803 

0.2 

0.960 

0.765 

3.0  /  14.4 
(2700) 

(540) 

0.1 

0.960 

0.874 

0.2 

0.956 

0.841 

4.0  /  19.2 
(3600) 

(720) 

0.1 

0.960 

0.873 

0.2 

0.957 

0.839 

5.0  /  24 
(4500) 

(900) 

0.1 

0.959 

0.887 

0.2 

0.956 

0.861 

10.0  /  48 
(9000) 

(1800) 

0.1 

0.959 

0.891 

0.2 

0.956 

0.862 

20.0  /  96 
(18000) 
(3600) 

0.1 

0.959 

0.892 

0.2 

0.955 

0.867 

0.1 

0.960 

0.877 

0.2 

0.956 

0.862 

Table  6A  ;  4  Exp  and  4  Wei  (Mixed)  in  Series.  RS  =  0.958  (Hi)  (•)  (Cont...) 

A(exp)  =  0.001  f/hr,  UT(exp)  =  5  hrs 
X(wei)  =  0.005  f/hr,  UT(wei)  =  15  hrs 


Degrees 

of 

Freedom 

a 

Measures  of  Accuracy 

RSLOW 

LEVEL 

1.3* 

2*(1  +  NFC) 

0.25  /  1.2 
(225) 

(45) 

0.1 

1.000 

0.723 

0.2 

0.994 

0.697 

0.5  /  2.4 
(450) 

(90) 

0.1 

0.977 

0.730 

0.2 

0.969 

0.681 

1.0  /  4.8 
(900) 

(180) 

0.1 

0.964 

0.853 

0.2 

0.955 

0.837 

2.0  /  9.6 
(1800) 

(360) 

0.1 

0.954 

0.934 

0.2 

0.949 

0.921 

3.0  /  14.4 
(2700) 

(540) 

0.1 

0.949 

0.961 

0.2 

0544 

0.949 

4.0  /  19.2 
(3600) 

(720) 

0.1 

0.948 

0.975 

0.2 

0.945 

0.966 

5.0  /  24 
(4500) 

(900) 

0.1 

0.948 

0.984 

0.2 

0.944 

0.978 

10.0  /  48 
(9000) 

(1800) 

0.1 

0.948 

0.971 

0.2 

0.943 

0.965 

20.0  /  % 
(18000) 
(3600) 

0.1 

0.948 

0.983 

0.2 

0542 

0.979 

IHKSSSMH 

0.1 

0.949 

0.978 

0.2 

0.944 

0.971 
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Table  6B  :  4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.827  (Lo)  (*) 

A(exp)  =  0.005  f/hr,  UT(exp)  =  5  hrs 
X(wei)  =  0.010  f/hr,  UT(wei)  =  15  hrs 


S/N 

Degrees 

of 

1 

Freedom 

0.25  / 1.2 
(45) 


0.5  /  2.4 
(90) 
(45) 


1.0  /  4.8 
(180) 
(90) 


3.0  /  14.4 
(540) 
(270) 


4.0  /  19.2 
(720) 
(360) 


10.0  /  48 
(1800) 
(900) 


20.0  /% 
(3600) 
(1800) 


Measures  of  Accuracy  | 

a 

RSLOW 

LEVEL 

0.1 

0.971 

0.691 

0.2 

0.938 

0,685 

0.1 

0.916 

0.600 

0.2 

0.893 

0.532 

0.1 

0.876 

0.684 

0.2 

0.855 

0.638 

0.1 

0.849 

0.787 

0.2 

0.836 

0.750 

0.1 

0.838 

0.855 

0.2 

0.824 

0.817 

0.1 

0.834 

0.870 

0.2 

0.824 

0.819 

0.1 

0.829 

0.890 

0.2 

0.822 

0.838 

0.1 

0.827 

0.900 

0.2 

0.820 

0.863 

0.1 

0.826 

0.902 

0.2 

0.819 

0.871 

0.1 

0.831 

0.888 

0.2 

0.821 

0.854 

Table  6B  ; 


4  Exp  and  4  Wei  (Mixed)  in  Series,  RS  =  0.827  (Lo)  (•)  (Cont...) 
>(exp)  =  0.005  f/hr,  UT(exp)  =  5  hrs 
X(wei)  =  0.010  f/hr,  UT(wei)  =  15  hrs 


S/N 

Degrees 

of 

Freedom 

Measures  of  Accuracy 

Q 

RSLOW 

LEVEL 

2 

1.3* 

2*(1  +  NFC) 

mmam 

0.1 

0.964 

0.718 

0.2 

0.922 

0.700 

0.5  /  2.4 
(90) 

(45) 

0.1 

0.895 

0.728 

0.2 

0.866 

0.688 

1.0  /  4.8 
(180) 

(90) 

0.1 

0.846 

0.864 

0.2 

0.818 

0.838 

2.0  /  9.6 
(360) 

(180) 

0.1 

0.811 

0.939 

0.2 

0.795 

0.921 

3.0  /  14.4 
(540) 

(270) 

0.1 

0.798 

0.964 

0.2 

0.780 

0.952 

4.0  /  19.2 
(720) 

(360) 

0.1 

0.792 

0.981 

0.2 

0.780 

0.972 

5.0  /  24 
(900) 

(450) 

0.1 

0.787 

0.986 

0.2 

0.777 

0.984 

10.0  /  48 
(1800) 

(900) 

0.1 

0.784 

0.988 

0.2 

0.774 

0.983 

20.0  /  96 
(3600) 

(1800) 

0.1 

0.783 

0.993 

0.2 

0.772 

0.988 

30.0  /  144 
(5400) 

(2700) 

0.1 

0.788 

0.996 

0.2 

0.775 

0.994 
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